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ABSTRACT 

Using the photometric parallax method we estimate the distances to ~48 million stars detected by 
the Sloan Digital Sky Survey (SDSS) and map their three-dimensional number density distribution in 
the Galaxy. The currently available data sample the distance range from 100 pc to 20 kpc and cover 
6,500 deg^ of sky, mostly at high galactic latitudes (|&| > 25). These stellar number density maps 
allow an investigation of the Galactic structure with no a priori assumptions about the functional 
form of its components. The data show strong evidence for a Galaxy consisting of an oblate halo, a 
disk component, and a number of localized overdensities. The number density distribution of stars 
as traced by M dwarfs in the Solar neighborhood [D < 2 kpc) is well fit by two exponential disks 
(the thin and thick disk) with scale heights and lengths, bias-corrected for an assumed 35% binary 
fraction, of Hi = 300 pc and Li = 2600 pc, and H2 = 900 pc and L2 = 3600 pc, and local thick-to- 
thin disk density normalization pthick^Ro) / pthiniRo)) — 12%. We use the stars near main-sequence 
turnoff to measure the shape of the Galactic halo. We find a strong preference for oblate halo models, 
with best-fit axis ratio c/a = 0.64, pH oc r^^-^ power-law profile, and the local halo-to-thin disk 
normalization of 0.5%. Based on a series of Monte-Carlo simulations, we estimate the errors of 
derived model parameters not to be larger than ~ 20% for the disk scales and ~ 10% for the density 
normalization, with largest contributions to error coming from the uncertainty in calibration of the 
photometric parallax relation and poorly constrained binary fraction. While generally consistent with 
the above model, the measured density distribution shows a number of statistically significant localized 
deviations. In addition to known features, such as the Monoceros stream, we detect two overdensities 
in the thick disk region at cylindrical galactocentric radii and heights {R,Z) ^ (6.5,1.5) kpc and 
{R, Z) ~ (9.5,0.8) kpc, and a remarkable density enhancement in the halo covering over a thousand 
square degrees of sky towards the constellation of Virgo, at distances of ~6-20 kpc. Compared to 
counts in a region symmetric with respect to the I = 0° line and with the same Galactic latitude, 
the Virgo overdensity is responsible for a factor of 2 number density excess, and may be a nearby 
tidal stream or a low-surface brightness dwarf galaxy merging with the Milky Way. The u — g color 
distribution of stars associated with it implies metallicity lower than that of thick disk stars, and 
consistent with the halo metallicity distribution. After removal of the resolved overdensities, the 
remaining data are consistent with a smooth density distribution; we detect no evidence of further 
unresolved clumpy substructure at scales ranging from ~ 50 pc in the disk, to '--^ 1 — 2 kpc in the halo. 
Subject headings: Galaxy: disk. Galaxy: halo, Galaxy: structure, Galaxy: fundamental parameters 
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In the canonical model of Milky Way formation 
(|Eggen. Lvnden-Bell. fc Sandaed 11962 ) the Galaxy be- 
gan with a relatively rapid (~ lO^yr) radial collapse of 
the initial protogalactic cloud, followed by an equally 
rapid settling of gas into a rotating disk. This model 
readily explained the origin and general structural, 
kinematic and metallicity correlations of observation- 
ally identified p opulations of field stars (jBaadd 119441 : 
lO'Connell |1958( ): low metallicity Population II stars 
formed during the initial collapse and populate the ex- 
tended stellar halo; younger Population I and Intermedi- 
ate Population II stars formed after the gas has settled 
into the Galactic plane and constitute the disk. 

The observationally determined distribution of disk 
stars is commonly desc r ibed by exponential density 
laws (IBahcall fc Soneiral fT980t IGilmore fc ReidI 119831 : 
IGilmore et al.lll989D . while power-laws or flattened de 
Vaucoule urs spheroids are usuall y used to describe the 
halo fe.g.- lWvse fc Gilmord[l989l : lLarsenlfT99l see also a 
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review bv lMaiewskilll99^ . In both disk and the halo, the 
distribution of stars is expected to be a smooth function 
of position, perturbed only slightly by localized bursts of 
star formation or spiral structure induced shocks. 

However, fo r some time, starting with the pioneer- 
ing work of ISearle fc ZinnI (|l978) . continuing with 
the studies of stellar counts and count asymmetries 
from Palomar Observatory Sky S urvey (e.g. iLarsenl 
[l99l iLarsen fc Humphrevd [ToM iParker et all I2003D . 
and most recently with the data from modern large- 
scale sky surveys (e .g., the Sloan Digital Sky Sur- 
vey, lYork et al.ri2000l: The Two M icron All Sky Sur- 
vey, 2 MASS. iMaiewski et all [20mL and the QUEST 
survey iVivas et all I2001D evidence has been mounting 
for a more complex picture of the Galaxy and its for- 
mation. Unlike the smooth distribution easily cap- 
tured by analytic density laws, new data argue for 
much more irregular substructure, especially in the 
stellar halo. Exa mples include the Sgr dwarf tidal 
stream in the halo (Ivezic et al. 2000: Ya nnv et al.ll2000l : 
IVivas et all 120011: [Majewski et al. 2003) or the Mono- 
ceros stream closer to the Gala ctic plane ()Newberg et al.l 
I2002t iRocha-Pinto et al.ll2003h . The existence of ongo- 
ing merging points to a likely significant role of accretion 
events in the early formation of the Milky Way's compo- 
nents, making the understanding of both the distribution 
of merger remnants, and of overall Milky Way's stellar 
content, of considerable theoretical interest. 



The majority (> 90%) of Galactic stellar content re- 
sides in the form of main-sequence (MS) stars. How- 
ever, a direct measurement of their spatial distribution 
requires accurate estimates of stellar distances to faint 
flux levels, as even the most luminous main sequence 
stars have V ^ 15 — 18 for the 1-10 kpc distance range. 
This requirement, combined with the need to cover a 
large sky area to sample a representative portion of the 
Galaxy, have historically made this type of measurement 
a formidable task. 

A common workaround to the first part of the problem 
is to use bright tracers for which reasonably accurate dis- 
tance estimates are possible (e.g. RR Lyrae stars, A-type 
stars, M giants), and which arc thought to correlate with 
the overall stellar number density distribution. These 
tracers, however, represent only a tiny fraction of stars 
on the sky, and their low number density prevents tight 
constraints on th e Galactic structure model parameters 
(jReid et al.ll 19961 ). For the same reason, such tracers are 
unsuitable tools for finding localized overdensities with 
small contrast ratios over their surroundings. 

Preferably, given precise enough multiband photome- 
try, one would avoid the use of tracers and estimate the 
distances to MS stars directly using a color- absolute mag- 
nitude, or "photometric parallax", relation. However, 
until now the lack of deep, large-area optical^^ surveys 
with sufficiently accurate multiband photometry has pre- 
vented an efficient use of this method. 

For example, near-IR colors measured by the all-sky 2MASS 
survey are not well suited for this purpose, because they only probe 
the Rayleigh- Jeans tail of the stellar spectral energy distribution 
and thus are not very sensitive to the effective temperature. 



Surveying a wide area is of particular importance. For 
example, even the largest Galactic structure oriented 
data set to date to use accurate optical CCD photom- 
etry (jSiegel et al.l 1200^ covered only ~ 15 deg^, with 
~ 10^ stars. To recover the overall Galactic density field 
their study, as others before it, it has had to resort to 
model fitting and assume a high degree of regularity in 
the density distribution and its functional form. This 
however, given that typical disk-|-halo models can have 
up to 10 free parameters, makes parameter estimation 
vulnerable to bias by unrecognized clumpy substructure. 

Indeed, a significant spread in results coming from 
different studies has existe d for quite some time (e.g., 
ISiegel et al.l |2002| . Table 1; iBilir et al.l [20061 ) indicating 
that either the unidentified substructures are confusing 
the model fits, that there are a multitude of degenerate 
models that are impossible to differentiate from using a 
limited number of lines of sight, or that the usual mod- 
els provide an inappropriate description of the large-scale 
distribution of stars in the Galaxy. A direct model-free 
determination of the stellar number density distribution 
in a large volume of the Galaxy would shed light on, and 
possibly resolve, all these issues. 

The large area covered by the SDSS, with accurate 
photometric measurements (~0.02 mag) and faint flux 
limits (r < 22), allows for a novel approach to studies of 
the stellar distribution in the Galaxy: using a photomet- 
ric parallax relation appropriate for main sequence stars, 
we estimate distances for a large number of stars and di- 
rectly map the Galactic stellar number density without 
the need for an a-priori model assumption^^. In this 
paper, we describe a study based on ^--^48 million stars 
detected by the SDSS in - 6500 deg^ of sky. An ad- 
vantage of this approach is that the number density of 
stars as a function of color and position in the Galaxy, 
p{X, Y, Z, r—i) can be measured without assuming a par- 
ticular Galactic model (e.g. the luminosity function and 
functional forms that describe the density laws for disks 
and halo). Rather, with minimal assumptions about the 
nature of the observed stellar population (that the large 
majority of the observed stars are on the main sequence) 
and by using an adequate photometric parallax relation, 
the computed stellar number density maps can be used to 
get an overall picture about the distribution of stars first, 
and a-posteriori constrain the density laws of Galactic 
components and look for deviations from them. 

This is the first paper, in a series of three^^, that em- 
ploys SDSS data and a photometric parallax relation to 
map the Galaxy. Here, we focus on the stellar number 
density distribution. In Ivezic et al. (2007, in prep., 
hereafter Paper II) we discuss the distribution of photo- 
metric metallicity (calibrated using SDSS spectroscopic 
data), and in Bond et al. (2007, in prep., hereafter Pa- 
per HI) we analyze the stellar kinematics using radial 
velocity and proper motion measurements. 

We begin by describing the SDSS data, the photo- 

The use of photometric parallax to determine Galactic model 
parameters is n ot particularl y novel, having a long history going 
back to at least iGilmore fcR cid (1983). The novelty in our ap- 
proach is to use the photometric parallax and wide area of SDSS 
to construct stellar density distribution maps first, and look for 
structure in the maps and fit analytic Galactic models second. 

We credit the late J.R.R. Tolkien for demonstrating the 
virtues of this approach. 
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metric parallax relations, and the construction of stellar 
number density maps in the following Section. Analysis 
of overall trends and identification of localized density 
features (substructure) is described in Section [3l In Sec- 
tion |4] we use the maps to derive best-fit parameters of 
density model for the Galactic disk and stellar halo. Sec- 
tionOdiscusses the details of a remarkably large overden- 
sity of stars identified in Section [3l Our results and their 
theoretical implications are summarized and discussed in 
Section [HI 

2. DATA AND METHODOLOGY 

In this Section we list the basic characteristics of the 
SDSS imaging survey, discuss the adopted photomet- 
ric parallax relation used to estimate the distance to 
each star, and describe a method for determining three- 
dimensional number density distribution as a function of 
Galactic coordinates. 

2.1. The Basic Characteristics of the SDSS Imaging 

Survey 

The SDSS is a digital photometric and spectroscopic 
survey which will cover up to one quarter of the celes- 
tial sphere in the North Galactic cap, and produce a 
smaller area (~ 225 deg^) but much deeper survey in 
the Southern Galactic he misphere^^ ( | York et all 120001: 
^oughton ct al. 2002; Abaz aiian et al.ll2003l:lGunn et al l 
2006t iTucker et al.. 20061. The flux densities of detected 
objects are measured almost simultaneously in five bands 
(u, g, r, i, and z) with effective wavelengths of 3540 A, 
4760 A, 6280 A, 7690 A, and 9250 A (iFukugit a et al" 



1996t iGunn et all fl998l : iSmith et all [2002; Hogg^eLal 



20011) . The completeness of SDSS catalogs for point 
sources is ~99.3% at t he bright end (r ~ 14, where the 
SDSS CCDs saturate. Ilvezi6 et al.ll200l . and drops to 
95% at magnitudes!'^ of 22.1, 22.4, 22.1, 21.2, and 20.3 
in u, g, r, i and z, res pectively. All magni tudes are given 
on the ABy system (Okc fc GunnI Il983l for additional 
discussion regarding the SDSS pho tometric system see 
IFukugita et all Il996l and 119991 ). The final survey 
sky coverage of about 8,000 deg^ will result in photo- 
metric measurements to the above detection limits for 
about 80 million stars and a similar number of galax- 
ies. Astrometric positions are accurate to about 0.1 arc- 
sec per coordina te for sources brighter than r ^20.5™ 
(jPier et al.l l2003l ). and the morphological information 
from the images allows robust point sourc e-galaxy sep- 
aration to r ~ 21.5™ (|Lupton et al.ll200l ). The SDSS 
photometric accuracy is 0.02 mag (root-mean-square, at 
the bright e nd), with well co ntrolled tails of the error 
distribution (jlvezic et al.ll2003| ). The absolute zero point 
calibration o f the SDSS photorn etry is accurate to within 
'--^ 0.02 mag ([Ivezic et al.ll2004bl) . A compendium o f tech- 
nical d etails about SDSS can be found in .Stoughton et all 
(|2002t l. and on the SDSS web site ( |http:// www.sdss.org 



2.2. The Photometric Parallax Method 

SDSS is superior to previous optical sky surveys be- 
cause of its high catalog completeness and accurate 

See also http;/ /www. astro. princeton.edu/PBOOK/welcome. htm' 
These values are determined by comparing multiple scans of 

the same area obtained during the commissioning year. Typical 

seeing in these observations was l.SitO.l arcsec. 
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Fig. 1. — A comparison of photometric parallax relations, ex- 
pressed in the Johnson system, from the literature. The relation 
from Henry et al. (1999) is valid for stars closer than 10 pc, while 
other relations correspond to the Hyades main sequence. Note that 
the latter differ by a few tenths of a magnitude. The relation from 
Laird, Carney & Latham (1988) is also shown when corrected for 
two different metallicity values, as marked in the legend. The gra- 
dient dMy /d[Fe/ H] given by their prescription is about 1 mag/dex 
at the blue end, and about half this value at the red end. 
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Fig. 2. — A comparison of photometric parallax relations in 
the SDSS ugriz system from the literature and adopted in this 
work. The two relations adopted here are shown by the dashed 
("bright" normalization) and solid ("faint" normalization) linos. 
Other lines show photometric parallax relations from the literature, 
as marked. The lower (thin) curve from Siegel et al. corresponds to 
low metallicity stars. The large symbols show SDSS observations 
of globular cluster M13. 

multi-band CCD photometry to faint flux limits over a 
large sky area. The majority of stars de tected by SDSS 
are main-sequence stars (~98%, iFinlator et al.l [2000') . 
which have a fairly well-deflned color-luminosity rela- 
tion^''. Thus, accurate SDSS colors can be used to esti- 
mate luminosity, and hence, distance, for each individual 
star. While these estimates are incorrect for a fraction 
of stars such as multiple systems and non-main sequence 
stars, the overall contamination is small or controllable. 
There are a number of proposed photometric parallax 

The uniqueness of color-luminosity relation breaks down for 
stars at main sequence tu rn-off (r — t ~ .11 mag for disk, and 
r - i ~ 0.06 for halo stars. [Chen et al.ll200lh . Those are outside of 
all but the bluest bin of the r — i range studied here. 
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relations in the literature. They differ in the methodol- 
ogy used to derive them (e.g., geometric parallax mea- 
surements, fits to globular cluster color-magnitude se- 
quences), photometric systems, and the absolute magni- 
tude and metallicity range for which they are applicable. 
Not all of them are mutually consistent, and most exhibit 
significant intrinsic scatter of order a half a magnitude 
or more. Even the relations corresponding to the same 
cluster, such as the Hyades, can differ by a few tenths of 
a magnitude (see Fig.[T]). 

In Fig. [2] we compare several recent photometric paral- 
lax relations found in the literature. They are all based 
on geometric parallax measurements, but the stellar col- 
ors are measured in different photometric systems. In or- 
der to facilitate comparison, we use photometric transfor- 
mations between the Johns on and SPSS system s derived 
for main-sequence stars by llvezic et"aLl (|2Q07cD , and fits 
to the stellar locus in SDSS color-color diagrams from 
llvezic et al.l (|2004b[ ). As evident, different photometric 
parallax relations from the literature are discrepant at 
the level of several tent hs to a magnitude. Fu rthermore, 
the relation proposed bv I Williams et all (|2002r) is a piece- 
wise fit to restricted color ranges, and results in a discon- 
tinuous relation. The behavior of Kurucz model atmo- 
spheres suggests that these discontinuities are probably 
unphysical. 

We constructed a fit, shown in Figure [2l that attempts 
to reconcile the differences between these relations. We 
require a low-order polynomial fit that is roughly consis- 
tent with the three relations at the red end, and properly 
reproduces the SDSS observations of the position of the 
turn-off (median Mr = 5 at r — z = 0.10 ) for globula r 
cluster M13 (using a distance of 7.1 kpc, Illarrisl[l996l ). 
The adopted relation 



Mr 



4.0 -f 11.86 {r-i)- 10.74 (r - if 
+bM {r-if - 1.20 {r-if 



(1) 



is very similar to the IWilliams et al (l2002f) rela t ion at 
the red end, and agrees well with the lSiegel et al.l ()2002f ) 
relation at the blue end. 

In order to keep track of uncertainties in our results 
due to systematic errors in photometric parallax rela- 
tion, we adopt another relation. The absolute magnitude 
difference between the two relations covers the plausible 
uncertainty range, and hence the final results are also ex- 
pected to bracket the truth. While we could arbitrarily 
shift the normalization of eq. [T] for this purpose, we in- 
stead use a relation that has an independent motivation. 

In Paper III, we propose a novel method to constrain 
the photometric parallax relation using kinematic data. 
The method relies on the large sky coverage by SDSS 
and simultaneous availability of both radial velocities and 
proper motion data for a large number of stars. These 
data can be successfully modeled using simple models 
such as a non-rotating halo and a disk rotational lag that 
is dependent only on the height above the Galactic plane. 
The best-fit models that are independently constrained 
using radial velocity and proper motion measurements 
agree only if the photometric parallax relation is correct. 
That is, the tangential velocity components, that are pro- 
portional to distance and measured proper motions, are 
tied to the radial velocity scale by adopting an appropri- 
ate distance scale. As a result of such kinematic analysis. 



we adopt a second photometric parallax relation 

Mr = 3.2 + 13.30 {r-i) - 11.50 [r - if (2) 
+5.40 {r-if {r-if. 

This relation is 0.66 mag brighter at the blue end {r — i^ 
0.1), and matches eq. [l]at r — z = 1.38 (see Fig. [2] for a 
comparison). The normalization differences between the 
two relations at the blue end correspond to a systematic 
distance scale change of ±18%, relative to their mean. 

To distinguish the two relations, we refer to the rela- 
tion from eq. [T] as the "faint" normalization, and to the 
relation from eq. [2] as the "bright" normalization. We 
note that, encouraging ly, the Hipparc os-hased Mn vs. 
R — I relation from iReid et al.l ()2001[ ) falls in between 
these two relations. 

In sections to follow, we perform all the analysis sep- 
arately for each relation, and discuss the differences in 
results when they are noticeable. For all figures, we use 
the bright normalization, unless noted otherwise. 

Equations [1] and [2] are quite steep, for example, 
AMr/A{r — i) ~ 10 mag/mag at the blue end (r — i ~ 
0.1). Because of this steepness^*, very accurate photom- 
etry (0.01-0.02 mag) is required to reach the intrinsic 
accuracy of the photometric relation (about 0.2 mag or 
better for individual globular clusters; for metallicity ef- 
fects see below). Older photographic surveys have pho- 
tometric errors of ^ 0.1-0.2 mag (Sesar et al. 2006), and 
inaccurate color measurements result in Mr errors ex- 
ceeding ~1 mag. Hence, with the SDSS, the intrinsic 
accuracy of the photometric parallax method can be ap- 
proached to a faint flux limit and over a large sky area 
for the first time. 

2.2.1. Effects of Metallicity on the Photometric Parallax 

Relation 

The main source of systematic errors in photometric 
parallax relation is it s depen dence on metallicity. For 
example, ISiegel et al.l (|2002f ) address this problem by 
adopting different relations for low- and high-metallicity 
stars (c.f. Fig. Another approach is to estimate 
metallicity, either from a spectrum or using photomet- 
ric methods such as a UV excess based S method (e.g. 
ICarnevll979( l. and then apply a correction to the adopted 
photometric parallax relation that depends both on color 
and metallicity (e.g. iLaird et al. 19 8^), as illustrated in 
Fig. [H We have tested the iLaird. Carney. &: LathamI 
metallicity correct ion by scaling the Hyades main se- 
quence, as given bv lKaraali et al.l (|2003D . using [Fe/H] — 
— 1.5 appropriate for M13, and comparing it to SDSS ob- 
servations of that cluster. As shown in Fig.|3l the agree- 
ment is very good (^0.1 mag). 

An application of 6 method to SDSS photome tric sys- 
tem was recently attempted by iKaraali et al.l (^005). 
However, as they pointed out, their study was not based 
on SDSS data, and thus even small differences between 
different photometric systems may have a significant ef- 
fect on derived metallicities (especially since the SDSS u 
band photometry is n ot precisely on the AB system, see 
lEisenstein et al.H2006D . 

This is not an artifact of the SDSS photometric system, or the 
adopted photometric parallax relation. For example, even for the 



linear My vs. B-V relation from lLaird et al.l JTOSl) dMv/d(B - 
V) = 5.6 mag/mag 
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Fig. 3. — A comparison of photometric parallax relations from 
the literature and adopted in this work, shown using blue SDSS 
bands (stars with spectral types later than ~M0 have g — r ~ 1.4). 
The two relations adopted here are shown by the dotted ("bright" 
normalization) and solid ("faint" normalization) lines. These are 
the same relations as shown in Fig. [2] translated here into the 
Mg{g — r) form using the r — i = f{g — r) relation appropriate for 
main sequence stars on the main stellar locus. Other lines show 
photometric parallax relations from the literature, as marked. The 
line marked Girardi et al. shows the range of model colors for 
Mg = 6. The lower (thin) curve from Siegel et al. corresponds 
to low metallicity stars. The triangle, circle and square show the 
SDSS observations of globular clusters Pal 5 ({Fe/ H\ = —1.4), and 
the Hyades [[Fe/H] = 0.1) and M48 ([Fe/H\ = -0.2) open clus- 
ters, respectively. The three large dots show the SDSS observations 
of globular cluster M13 {[Fe/H\ = —1.5). Note the good agree- 
ment between these observations and the Hyades sequence scaled 
to MlS's metallicity using the prescription from Laird, Carney & 
Latham (1988). For reference, B-V = 0.95 (g-r)-f 0.20 to within 
0.05 mag. 

The expected correlation of metallicity and the 
SD SS u — q an d g — r colors was recently discussed 
by llvezic et all ()2007bf ). Using SDSS photometry 
and metallicity estimates derived from SDSS spectra 
(|Allende Prieto et al.ll2006( ). they demonstrated a very 
strong dependence of the median metallicity on the po- 
sition in the g — r vs. u — g color-color diagram. For 
example, for stars at the blue tip of the stellar locus 
[u — g < 1, mostly F stars), the expression 

[Fe/i/] = 5.11 - 6.33 (3) 

reproduces the spectroscopic metallicity with an rms of 
only 0.3 dex. This relation shows that even in this favor- 
able case (it is much harder to estimate metallicity for 
red stars), a 0.1 mag error of the u — g color would in- 
troduce an error of [Fe/ H] as large as 0.5 dex, resulting 
in an error in the absolute magnitude of ^--^0.5 mag. 

We aim here to study the Galaxy to as large a distance 
limit as the SDSS photometric sample of stars allows. 
While metallicity could be estimated for bright blue stars 
using the above expression, for most stars in the sample 
the SDSS u band photometry is not sufficiently accurate 
to do so reliably. For example, the ra ndom error of u — q 
color becomes 0.1 mag at m '--^ 20.5 (jlvezic et al.ll2003l ). 
which corresponds to ~ 19.5 or brighter even for the 
bluest stars. Therefore, metallicity estimates based on 
the u — g color would come at the expense of a more than 
2 mag shallower sample. Hence, we choose not to correct 
the adopted photometric parallax relation for metallicity 
effects, and only utilize the correlation between metal- 



licity and u — g color when constraining the metallicity 
distribution of a large halo overdensity discussed in Sec- 
tion [5] 

We point out that the adopted relations do account for 
metallicity effects to some extent. The metallicity distri- 
bution shows a much larger gradient perpendicular to the 
Galactic plane tha n in the radial direction (see Fig. 3 in 
llvezic et all 120079 ) ■ As we only consider high Galactic 
latitude data, the height above the plane is roughly pro- 
portional to distance. At the red end, the adopted rela- 
tions are tied via geometric parallax to nearby metal-rich 
stars, and even the faintest M dwarfs in SDSS sample are 
only ~1 kpc away. At the blue end, the adopted relations 
are tied to globular clusters and halo kinematics, which 
is appropriate for the bluest stars in the sample, that 
are detected at distances from several kpc to ~10 kpc. 
Thus, in some loose "mean" sense, the adopted relation 
smoothly varies from a relation appropriate for nearby, 
red, high-metallicity stars to a relation appropriate for 
more distant, blu e, low-metallicity stars^^. Furthermore, 
iReid et all ()2001[ ) show that photometric parallax rela- 
tions constructed using red photometric bands, such as 
our Mr vs. r ■— i relation, are much less sensitive to 
metallicity than the traditional My vs. B — V relation 
(compare the top left and bottom right panel in their 
Fig. 15). 

Nevertheless, to further control metallicity and other 
systematic effects, we perform analysis in narrow color 
bins, as described in more detail in Section l2?3l 

2.2.2. A Test of the Photometric Parallax Relation using 
Resolved Binary Stars 

The number of close stellar pairs in the SDSS survey 
with distances in the 2-5 arcsec range shows an excess 
relative to the extrapolation from larger distances (Sesar 
et al. 2007, accepted to ApJ). Statistically, they find 
that ^70% of such pairs are physically associated bi- 
nary systems. Since they typically have different col- 
ors, they also have different absolute magnitudes. The 
difference in absolute magnitudes, AM, can be com- 
puted from an adopted photometric parallax relation 
without the knowledge of the system's distance, and 
should agree with the measured difference of their appar- 
ent magnitudes, Am. The distribution of the difference, 
S — Am — AM should be centered on zero and should 
not be correlated with color if the shape of photometric 
parallax relation is correct (the overall normalization is 
not constrained, but this is not an issue since the rela- 
tion can be anchored at the red end using nearby stars 
with geometric parallaxes) The width of the S dis- 
tribution provides an upper limit for the intrinsic error 
of the photometric parallax method (note, however, that 
S is not sensitive to systematic errors due to metallicity 
since the binary components presumably have the same 
metallicity) . 

When the adopted photometric parallax relation is applied 
to the Sun {r — i = 0.10), the resulting absolute magnitude is too 
faint by about 0.5 mag. This is an expected result, because the 
relation is anchored to a low-metallicity globular cluster at the blue 
end. For both relations, the predicted absolute magnitudes of low- 
temperature, low-metallicity stars are systematically too bright. 
However, the majority of such stars (e.g., distant halo M-dwarfs) 
are faint, and well beyond the flux limit of the survey. 

Note the similar ities of this m ethod, and the method of re- 
duced proper motions ILuvtenI HimB ). 
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Fig. 4. — The distribution of the median <5 for a sam- 
ple of ~ 17,000 candidate wide-angle binaries in the (r — i)i 
(color of brighter pair member; the primary) vs. (r — 1)2 (color 
of fainter member; the secondary) color-color diagram. Here, 
5 = {Mr,2 — Mr,i) — (r2 — Ti), IS the difference of two estimates 
(one from the absolute, and the other from the apparent magni- 
tudes) of brightness difference between the two components. In 
the top panel, the absolute magnitudes were estimated using eq. [2] 
(the "bright" paralax relation; the dotted line in Figure [2]l, and 
in the bottom panel using eq. [T] (the "faint" paralax relation; the 
solid line in Figure [2]l Inset histograms show the distribution of 
the median S evaluated for each color-color pixel. The distribution 
medians are 0.07 (top panel) and -0.004 (bottom panel), and the 
dispersions (determined from the interquartile range) are 0.13 and 
0.10 mag, respectively. 

We have performed such a test of adopted parallax re- 
lations using a sample of 17,000 candidate binaries from 
SDSS Data Release 5. Pairs of stars with 14 < r < 20 
are selected as candidate wide binaries if their angular 
separation is in the 3-4 arcsec range. The brighter star 
(in the r band) is designated as the primary (subscript 
1), and the fainter one as the secondary (subscript 2). 
For each pair, we calculated S twice - once assuming the 
bright photometric parallax relation (eq.[2]), and once as- 



suming the faint relation (eq.[T]). We further remove from 
the sample all pairs with \S\ > 0.5, those likely being the 
interlopers and not actual physical pairs. 

The results of this analysis are summarized in Fig- 
ure m The color-coded diagrams show the dependence 
of S on the r — i colors of the primary and the secondary 
components. The median S value in each {rii,ri2) pixel 
measures whether the absolute magnitude difference ob- 
tained using the parallax relation for stars of colors rii 
and ri2 is consistent with the difference of their appar- 
ent magnitudes (in each bin, the d distribution is much 
more peaked than for a random sample of stars, and is 
not much affected by the \S\ < 0.5 cut). If the shape of 
the photometric parallax relation is correct, the median 
d should be close to zero for all combinations of rii and 
ri2- 

The distributions of the median 5 for each pixel are 
fairly narrow (~ 0.1 mag), and centered close to zero 
(the medians are 0.07 mag for the bright relation and 
—0.004 mag for the faint relation). Irrespective of color 
and the choice of photometric parallax relation, the devi- 
ations are confined to the ~ ±0.25mag range, thus plac- 
ing a stringent upper limit on the errors in the shape of 
the adopted relations. The S distributions root-mean- 
square width of ^ 0.1 mag implies average distance error 
of about 5%. Nevertheless, the binary stars in a candi- 
date pair are of presumably identical metallicities. As a 
large fraction of the intrinsic scatter of Mr{r — i) comes 
the dependence of absolute magnitude on metallicity, we 
adopt a conservative value of ctm,. = 0.3. 

The coherent deviations seen in Figure 3] (e.g. around 
rii ~ 0.3 and ri2 ~ 0.5) indicate that the adopted par- 
allax relations could be improved. Given already satis- 
factory accuracy of the adopted relations, such a study 
is presented separately (Sesar et al. 2007, accepted to 
ApJ). 

2.2.3. Contamination by Giants 

The photometric parallax method is not without pit- 
falls, even when applied to the SDSS data. Tradition- 
ally, the application of this method was prone to sig- 
nificant errors due to sample contamination by evolved 
stars (subgiants and giants, hereafter giants for simplic- 
ity), and their underestimated distances. This effect is 
also present in this study, but at a much less significant 
level because of the faint magnitudes probed by SDSS. 
At these flux levels, the distances corresponding to giants 
are large and sometimes even beyond the presumed edge 
of the Galaxy (up to ~100 kpc). The stellar density at 
these distances is significantly smaller than at distances 
corresponding to main sequence stars with the same ap- 
parent magnitude. The contamination with evolved stars 
rapidly asymptotes (e.g., assuming a halo profile) 

and may decline when the edge of the halo is reached. 

A quantitative illustration of this effect is shown in 
Fig. \5\ for a fiducial Galaxy model. The worst case sce- 
nario corresponds to G giants with g — r ^ 0.4 — 0.5 and 
r — i ~ 0.15 — 0.20, and their most probable fraction is 
about 5%. This color range and the fraction of giants 
was determined using the SDSS data for the globular 
cluster M13 (the data for the globular cluster Pal 5 im- 
ply similar behavior). To be conservative, we have also 
tested a model with a twice as large fraction of giants. 
This analysis (see bottom panel) shows that the effect of 
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Fig. 5. — An illustration of the effects of misidentifying giants 
as main sequence stars. The top panel shows the Z dependence of 
stellar density at R=8 kpc for a fiducial model consisting of two 
disks with scale heights of 300 pc and 1200 pc. The contribution 
of the disks is shown by the short-dashed line, and the long-dashed 
line shows the contribution of a power-law spherical halo with the 
power-law index of 3. The middle panel shows the contribution 
of misidentified giants from disks (short-dashed) and halo (long- 
dashed) for an assumed giant fraction of 5%, and underestimated 
distances by a factor of 3. The "contaminated" model is shown by 
dotted line, just above the solid line, which is the same as the solid 
line in the top panel. The ratio of the "contaminated" and true 
density is shown in the bottom panel (note the different horizontal 
scale). 

misidentifying giants as main sequence stars is an overall 
bias in estimated number density of {^8 % when 

the fraction of giants is 10%), with little dependence on 
distance from the Galactic plane beyond 500 pc. This is 
the distance range probed by stars this blue, and thus the 
worst effect of contamination by giants is a small overall 
overestimate of the density normalization. Shorter dis- 
tances are probed by redder stars, M dwarfs, for which 
the contamination by M giants is negligible because the 
luminosity difference between red giants and dwarfs is 



very large (e.g. there are tens of millions of M dwarfs 
in our sample, while the 2MASS survey revealed only a 
few thousand M gian ts in the same region of the sky, 
iMaiewski et al.|[2003f) . Hence, the misidentified giants 
are not expected to significantly impact our analysis. 

2.2.4. Unrecognized Multiplicity 

Multiplicity may play a significant role by systemati- 
cally making unresolved multiple systems, when misiden- 
tified as a single star, appear closer then they truly are. 
The net effect of unrecognized multiplicity on derived 
distance scales, such as scale height and scale length, is 
to underestimate them by up to ~35% (see Section l4.3.4l 
here and Siegel et al. 2002). The magnitude of this bias 
is weakly dependent on the actual composition of the bi- 
naries (e.g. their color difference and luminosity ratio), 
but it is dependent on the fraction of multiple systems in 
the Galaxy. Since this fraction is not well constrained, 
for the purpose of constructing the number density maps 
(Section 12. 4p we assume all observed objects are single 
stars. This biases the distance scales measured off the 
maps, making them effectively lower limits, and wc a- 
posteriori correct for it, after making the Galactic model 
fits (Sections 14.3.41 and 14.3.11"]) . Note that this bias can- 
not affect the shapes of various density features seen in 
the maps, unless the properties of multiple systems varies 
greatly with the position in the Galaxy. 

2.2.5. Distance Range Accessible to SDSS Observations of 
Main-Sequence Stars 

A disadvantage of this method is its inability, when 
applied to main sequence stars, to probe distances as 
large as those probed by RR Lyrae and M giants (20 kpc 
vs. 100 kpc). However, a significant advantage of using 
main sequence stars is the vastly larger number of stars 
(the number ratio of main sequence to RR Lyrae stars 
in the SPSS sample is ^10,000, and even l arger for M 
giant s (jlvezic et al.ll2003l : llvezic et al.ll2003bt llvezic et al.l 
l2005f ). This large number of main-sequence stars allows 
us to study their number density distribution with a high 
spatial resolution, and without being limited by Poisson 
noise in a large fraction of the observed volume. 

2.3. The SDSS Stellar Sample 

In this Section we describe the stellar sample utilized in 
this work, and the methods used to construct the three- 
dimensional number density maps. 

2.3.1. The Observations 



We utilize observations from 248 SDSS 
runs obtained in a 5 year period through 
ber 



imagmg 

mat) year period tlirougli Septem- 
2003, which cover 6,538 deg^ of the sky This 
is a superset of imaging runs desc ribed in SDSS 
Data Release 3 (jAbazaiian et al.l 120051 ) , complemented 
by a number of runs from S DSS Data Release 4 
(|Adelman-McCarthv et all l2006l) and the so called 
"Orion" runs ( Finkbeiner et al.l l2004l ). The sky cov- 
erage of these 248 runs is shown in figure |6l They 
cover 5450 deg^ in the northern Galactic hemisphere, and 
1088 deg2 in the south. 

We start the sample selection with 122 million detec- 
tions classified as point s ources (stars) by the SDSS pho- 
tometric pipeline. Photo (|Lupton et al.ll200^ . For a star 
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Fig. 6. — Footprint on the sky of SDSS observations used in this work shown in Lambert equal area projection (hatched region). The 
circles represent contours of constant Galactic latitude, with the straight lines showing the location of constant Galactic longitude. For this 
study, observations from 248 SDSS imaging runs were used, obtained over the course of 5 years. The data cover 5450 deg^ of the north 
Galactic hemisphere, and a smaller but more frequently sampled area of 1088 deg^ in the southern Galactic hemisphere. 



to be included in the starting sample, we require that 
r < 22, and that it is also detected (above 5(t) in at least 
the g or i band. The latter requirement is necessary to be 
able to compute either the g — rovr — i color. The two re- 
quirements reduce the sample to 87 million observations. 
For each magnitude measurement. Photo also provides a 
fairly reliable estimate of its accuracy (Ivczic et al. 2003,) , 
hereafter Ug, Gr and ai. We correct all me asurements for 
the in terstellar dust extinction using the ISchlegel et al.l 
(IT998h (hereafter SFD) maps. 

2.3.2. T/ie Effects of Errors in Interstellar Extinction 
Corrections 

The SFD maps are believed to be correct within 10%, 
or better. This uncertainty plays only a minor role in this 
work because the interstellar extinction is fairly small at 
the high galactic latitudes analyzed here (|&| > 25): the 
median value of the extinction in the r band, A^, is 0.08, 
with 95% of the sample with Ar < 0.23 and 99% of the 
sample with Ar < 0.38. Thus, only about 5% of stars 
could have extinction correction uncertain by more than 
the photometric accuracy of SDSS data (~0.02 mag). 
The SFD maps do not provide the wavelength depen- 
dence of the interstellar correction, only its magnitude. 
The extinction corrections in the five SDSS photometric 
bands are computed from the SFD maps using conver- 
sion coefficients derived from an Ry = 3.1 dust model. 
Analysis of the position of the stellar locus in the SDSS 
color-color diagrams suggests that these coefficients are 
satisfactory at th e level of accuracy a nd galactic latitudes 
considered here (|Ivezic et al.ll2004b[ ). 

We apply full SFD extinction correction to all stars in 
the sample. This is inappropriate for the nearest stars 
because they are not beyond all the dust. Distances to 
the nearest stars in our sample, those with r — i ^ 1.5 
(the red limit) and r ^ 14 (approximately the SDSS r 
band saturation limit), are ^30 pc (distance determina- 



tion is described in the next two sections). Even when 
these stars are observed at high galactic latitudes, it is 
likely that they are over-corrected for the effects of in- 
terstellar extinction. To estimate at what distances this 
effect becomes important, we have examined the depen- 
dence of the g — r color on apparent magnitude for red 
stars, selected by the condition r — i > 0.9, in the re- 
gion defined by 210 </< 240 and 25 < 6 < 30. The 
distribution of the intrinsic g — r color for these stars is 
practically independent of their r — i color (see Fig. [8]) , 
with a median of 1.40 and a standard deviation of only 
0.06 mag (Ivezic et al. 2004b). This independence allows 
us to test at what magnitude (i.e. distance) the applied 
SFD extinction corrections become an overestimate be- 
cause, in such a case, they result in g — r colors that are 
bluer than the expected value of 1.40. We find that 
for r > 15 the median g — r color is nearly constant - it 
varies by less than 0.02 mag over the 15 < r < 20 range. 
On the other hand, for stars with r < 15 the median g — r 
color becomes much bluer - at r = 14.5 the median value 
is 1.35. This demonstrates that stars at r > 15 are al- 
ready behind most of the dust column. With the median 
r — i color of 1.17, the implied distance corresponding to 
r = 15 is ~80 pc. For the probed galactic latitude range, 
this indicates that practically all the dust is confined to a 
region within ~70 pc from the galactic midplane (here we 
define midplane as a plane parallel to the galactic plane 
that has Z — —25 pc, because the Sun is offset from 
the midplane towards the NGP by ~25 pc; for more de- 
tails see below). We arrive to the same conclusion about 
the dust distribution when using an analogous sample in 
the south galactic plane with \b\ ~ 12 (in this case the 
median g — r color is systematically bluer for r < 19, 
due to different projection effects and the Sun's offset 
from the midplane). Hence, in order to avoid the effects 
of overestimated interstellar extinction correction for the 
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Fig. 7. — The top panel shows the mean fractional distance error 
as a function of the r — i color and r band magnitude, assuming the 
intrinsic photometric parallax relation scatter of om^ = 0.3 mag. 
The solid lines are contours of constant fractional distance error, 
starting with ao/D = 15% (lower right) and increasing in incre- 
ments of 5% towards the top left corner. The dotted lines are 
contours of constant distance, and can be used to quickly estimate 
the distance errors for an arbitrary combination of color and mag- 
nitude/distance. Fractional distance errors are typically smaller 
than ~ 20%. Note that the distance errors act as a o-jvf^ wide 
convolution kernel in magnitude space, and leave intact structures 
larger than the kernel scale. In particular, they have little effect 
on the slowly varying Galactic density field and the determination 
of Galactic model parameters. 

nearest stars, we exclude stars that are within 100 pc 
from the galactic plane when fitting galaxy models (de- 
scribed below). Only 0.05% of stars in the sample are at 
such distances. In summary, the effects of overestimated 
interstellar extinction correction, just as the effects of 
sample contamination by giants, are not very important 
due to the faint magnitude range probed by SDSS. 

2.3.3. The Treatment of Repeated Observations 

SDSS imaging data are obtained by tracking the sky 
in six parallel scanlines, each 13.5 arcmin wide. The six 
scanlines from two runs are then interleaved to make a 
filled stripe. Because of the scan overlaps, and because of 
the convergence of the scans near the survey poles, about 
40% of the northern survey is surveyed at least twice. 
Additionally, the southern survey areas will be observed 
dozens of times to search for variable objects and, by 
stacking the frames, to push the fiux limit deeper. For 
these reasons, a significant fraction of measurements are 
repeated observations of the same stars. 

We positionally identify observations as corresponding 
to the same object if they are within 1 arcsec of each 
other (the median SDSS seeing in the r band is 1.4 arc- 
sec). Out of the initial ^^122 million observations, the 
magnitude cuts and positional matching produce a cat- 
alog of 47.7 million unique stars (the "star catalog". Ta- 
ble [T|). They span the MK spectral types from ^F9 to 
~M3 (Table [5]). There are two or more observations for 
about 36% (17.2 million) of observed stars. For stars with 
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Note. — Repeat observations in the stellar sample: Because 
of partial imaging scan overlaps and the convergence of scans 
near survey poles, a significant fraction of observations are re- 
peated observations of the same stars. In columns 7V(r < 22) 
and N{r < 21.5) we show the number of stars observed Napp 
times for stars with average magnitudes less than r = 22 and 
r = 21.5, respectively. The final two rows list the total number 
of stars in the samples, and the total number of observations. 

multiple observations we take the catalog magnitude of 
the star to be equal to the weighted mean of all obser- 
vations. In this step there is a tacit assumption that 
the variability is not important, justified by the main- 
sequence nature of the stellar sample under consideration 
(for the variabilit y analysis of the SDSS stellar sample see 
iSesar et al.l[2006l ). 

As discussed in Section [2. 21 an accurate determination 
of stellar distances by photometric parallax hinges on a 
good estimate of the stellar color and magnitude. In the 
bottom panel of Fig. [7] we show the mean r magnitude 
error of stars in the catalog as a function of the r band 
magnitude. The photometric errors are ^0.02 mag for 
bright objects (limited by errors in modeling the point 
spread function) , and steadily increase towards the faint 
end due to the photon noise. At the adopted sample 
limit, r — 22, the r band photometric errors are ~0.15 
mag. The g and i band magnitude errors display similar 
behavior as for the r band. 

2.3.4. Maximum Likelihood Estimates of True Stellar 

Colors 

The photometric parallax relation (eq. [2]) requires only 
the knowledge of r — i color to estimate the absolute mag- 
nitude. The accuracy of this estimate deteriorates at the 
faint end due to increased r — i measurement error. It 
also suffers for blue stars (r — i < 0.2) of all magnitudes 
because the slope of the photometric parallax relation, 
AA/r/A(r — i), is quite large at the blue end - for these 
stars it would be better to use the g — r {or u — g) color 
to parametrize the photometric parallax relation. On the 
other hand, the g — r color is constant for stars later than 
~M0 [g ~r ~ 1.4), and cannot be used for this purpose. 
These problems can be alleviated to some extent by uti- 
lizing the fact that colors of main sequence stars form a 
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Note. — The number density map parameters. Each of the 19 maps is a volume hmited three-dimensional density map of 
stars with rio < r — 2 < ri\, corresponding to MK spectral types and mean g — r column listed in columns SpT and < gr > , 
respectively. Median absolute magnitude Mr, distance limits Dq — Di (in parsecs) and binning pixel scale dx (also in parsecs) 
are given in columns labeled "Bright" and "Faint", for the bright (Equation [2j and faint (Equation photometric parallax 
relation. The number of stars in each r — i bin is given in Ngtars column (in millions). 
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blue stars and 0.06 mag for red stars flvezic "ei~aIll2004H ). 
which is significantly smaller than the measurement er- 
ror at the faint end. To a very good approximation, any 
deviation of observed colors from the locus can be at- 
tributed to photometric errors. We use this assumption 
to improve estimates of true stellar colors and apparent 
magnitudes at the faint end, and thus to increase the 
sample effective distance limit by nearly a factor of two. 

As illustrated in Fig. [51 for each point and a given 
error probability ellipse, we find a point on the locus with 
the highest probability^^, and adopt the corresponding 
(g— r)e and {r — i)e colors. The error ellipse is not aligned 
with the g — r and r — i axes because the g — r and r — i 
errors are correlated {u'^-r.r-i 



g,r 



0.5 



1.5 



Fig. 8. — The distribution of ~48 million stars analyzed in 
this work in the r — i vs. g — r color-color diagram, shown by 
isodensity contours. Most stars lie on a narrow locus, shown by 
the dashed line, whose width at the bright end is 0.02 mag for blue 
stars {g — r ^ 1) and 0.06 mag for red stars (3 — r ^ 1.4). The 
inserts illustrate the maximum likelihood method used to improve 
color estimates: the ellipses show measurement errors, and the 
crosses are the color estimates obtained by requiring that a star 
lies exactly on the stellar locus. Note that the principal axes of 
the error ellipses are not aligned with the axes of the color-color 
diagram because both colors include the r band magnitude. 

very narrow, nearly one-dimensional locus. 

The r — t vs. g — r color-color diagram of stars used 
in this work is shown in Fig. [8] We find that the stellar 
locus is well described by the following relation: 



We exclude all points further than 0.3 mag from the 
locus, as such large deviations are inconsistent with 
measurement errors, and in most cases indicate the 
source is not a main-sequence star. This requir ement 
effect ively removes hot white dwarfs (IKlcinma n et al 
2004), low-redshift quasars (z < 2.2, [Richards et al 



2002 , and wh ite dwarf/red dwarf unresolved binaries 
([Smolcic et ali l2004). 

Using the maximum likelihood colors, we estimate the 
magnitudes (g^, r^, ie) by minimizing: 

+ -5 + -2 , (5) 



X 



Or 



which results in 

Wrr -f w, 



+ "Wgia - (.9 - r)e) + Wjji + {r - i)e) 

Wr + Wo + Wi 



(6) 



g~r= 1.39(1 
-2Ah{r-if - 1. 



exp[-4.9(r - if 
8(r-i) - 0.050]) 



(4) 



which is shown by the solid line in the figure. 
The intrinsic width of the stellar locus is 0.02 mag for 



This is effectively a Bayesian maximum likelihood (ML) pro- 
cedure with the assumption of a uniform prior along the one- 
dimensional locus. As seen from from Fig [8] the real prior is not 
uniform. We have tested the effects of non-uniform priors. Adopt- 
ing an observationally determined (from FiglSJl non- uniform prior 
would change the loci of posterior maxima by only ^ 0.005 mag 
(worst case), while further complicating the ML procedure. We 
therefore retain the assumption of uniform prior. 
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9e {g - r)e + re (7) 
ie^{r~i)e-re (8) 

where Wj = V"'! ^'^^ J ~ 5: '"i 

The adopted (r — i)e color and rg magnitude uniquely 
determine (through eqs.[T]and[2]) the absolute magnitude 
Mr for each star in the catalog. We dub this procedure 
a "locus projection" method, and refer to the derived 
colors as "locus-projected colors" . In all subsequent cal- 
culations we use these "locus-projected" colors, unless 
explicitly stated otherwise. This method is the most nat- 
ural way to make use of all available color information, 
and performs well in cases where the measurement of 
one color is substantially worse than the other (or even 
nonexistent). It not only improves the color estimates 
at the faint end, but also helps with debiasing the esti- 
mate of density normalization in regions of high gradients 
in {g — r,r — i) color-color diagram (e.g., near turnoff). 
This and other aspects of locus projection are further 
discussed in Appendix A. 

2.3.5. The Contamination of Stellar Counts by Quasars 

The stellar samples selected using the g — r and r — i 
colors, as described above, are contaminated by low- 
redshift quasars. While easily recognizable with the aid 
oi u — g color, a significant fraction of quasars detected 
by SDSS have the g — r and r — i colors similar to 
those of turn-off stars. The SDSS sample of spectro- 
scopically confirmed quasars is flux- limited at i = 19.1 
(|Richards et al.l [20021 . and references therein) and thus 
it is not deep enough to assess the contamination level 
at the fai nt end relevant her e. Instead, we follow analy- 
sis from ( Ivezic et al.ll2004al) . who were interested in the 
contamination of quasar samples by stars, and obtain 
an approximate contamination level by comparing the 
counts of faint blue stars and photometrically selected 
quasar candidates. We use a catalog of coadded pho- 
tometry based on about ten repeated SDSS observations 
recently constructed by ([Ivezic et al.ll2007a[ ). The cata- 
log covers a 300 deg^ large sky region at high galactic 
latitudes (|&| ^ 60°) and thus the estimated contamina- 
tion fraction represents an upper limit. With its signifi- 
cantly improved u — g color measurements relative to sin- 
gle SDSS scans, this catalog allows efhcient photometric 
selection of low-redshift quasar candidates to flux levels 
below r = 21. 

As summarized in Fig. [HI the largest contamination 
of stellar sample by quasars is expected in blue bins. 
The bluest bin (0.10< r — i < 0.15) includes stars with 
0.2< g ~ r <0.5, and ~5% of sources in the r < 21.5 
subsample have u — r < 0.8, consistent with quasars. 
Even if we restrict the sample to 0.2< g — r <0.3, and 
thus maximize the sample contamination by quasars, the 
estimated fraction of quasars does not exceed 10% for 
r < 21.5 (see the bottom right panel). 

2.3.6. Estimation of Distances 

Given the photometric parallax relation (eq[l]), the 
locus-projected maximum likelihood r band magnitude, 
and r — i color, it is straightforward to determine the 
distance D to each star in the catalog using 

D = 10"^+^ pc, (9) 



Depending on color and the chosen photometric parallax 
relation, for the magnitude range probed by our sample 
(r=15-21.5) the distance varies from ^100 pc to ~20 
kpc. 

Due to photometric errors in color, magnitude, and the 
intrinsic scatter of the photometric parallax relation, the 
distance estimate has an uncertainty, cr^i, given by: 

2 / \2 2 l^^\2 2 m\ 

'^D = f^M.(r-^) + (^) '^r (Hj 

where ctm^ is the intrinsic scatter in the photometric par- 
allax relation. With an assumption of (T,^_j « 2(t^, this 
reduces to a simpler form: 

^ = 0-46^(l + 2(^|^)>,2 + .,2,^ (12) 

The fractional distance error, cro/D, is a function of 
color, apparent magnitude and magnitude error (which 
itself is a function of apparent magnitude). In the top 
panel in figure [7| we show the expected ao I D as a func- 
tion of r and r — i with an assumed intrinsic photomet- 
ric relation scatter of ctm^ = 0.3 mag. This figure is a 
handy reference for estimating the distance accuracy at 
any location in the density maps we shall introduce in 
Section 12.41 For example, a star with r — i = 0.5 and 
r = 20 (or, using eq. [2 at a distance of D = 3 kpc) has a 
^--^18% distance uncertainty. Equivalently, when the stars 
are binned to three-dimensional grids to produce density 
maps (Section l2.4p . this uncertainty gives rise to a nearly 
Gaussian kernel smoothing the maps in radial direction, 
with color and distance dependent variance a^. Note 
that this convolution leaves intact structures larger than 
the kernel scale and, in particular, has little effect on the 
slowly varying Galactic density field and determination 
of Galactic model parameters fSection l4.3.3p . 

To summarize, due to measurement errors, and uncer- 
tainty in the absolute calibration of the adopted photo- 
metric parallax relations, the derived density maps, de- 
scribed below, will differ from the true stellar distribu- 
tion. First, in the radial direction the spatial resolution 
is degraded due to the smoothing described above. A 
similar effect is produced by misidentification of binaries 
and multiple systems as single stars. Second, the dis- 
tance scale may have systematic errors, probably color 
and metallicity dependent, that "stretch or shrink" the 
density maps. Third, for a small fraction of stars, the dis- 
tance estimates may be grossly incorrect due to contam- 
ination by giants and multiple unresolved systems. Fi- 
nally, stars with metallicities significantly different than 
assumed at a particular r—i int the parallax relation may 
be systematically placed closer or farther away from the 
origin (the Sun). 

However, all of these are either small (e.g., contami- 
nation by giants), have a small total effect on the un- 
derlying Galactic density field (radial smearing due to 
dispersion in distance estimates), or cause relative radial 
displacements of entire clumps of stars with metallicities 
different than that of the background while not affecting 
their relative parallaxes, and thus allowing the discrim- 
ination of finer structure. Altogether, the maps fidelity 
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Fig. 9. — The dots in the top left panel shows point sources from the SDSS Stripe 82 catalog of coadded observations in the r vs. u — r 
color-magnitude diagram. Sources with 0.3 < g — r < 0.5 are marked blue. Sources with u — r< 0.8 are dominated by low-redshift quasars, 
those with m — r ~ 1.3 by low-metallicity halo stars, and the bright stars (r < 18) with ti — r ~ 1.6 are dominated by thick disk stars. Note 
the remarkable separation of halo and disk stars both in magnitude (a distance effect) and color (a metallicity effect) directions. The top 
right panel shows a subset of sources with r < 21 in the g — r vs. u — g color-color diagram. Cumulative counts of sources from several 
regions of this diagram (blue: hot stars, dominated by white dwarfs; red: quasars; magenta: blue horizontal branch stars; cyan: halo stars; 
green: thick disk stars) are shown in the lower left panel, with the same color coding. The solid lines have slopes of 0.13 (blue) and 0.34 
(red) for thick disk and halo stars, while the quasar counts change slope at r ~ 20 from ~0.7 to ~0.4, as indicated by the dashed lines. 
The bottom right panel compares cumulative counts of two subsets of sources with 0.2 < g — r < 0.3 that are separated by the u — r = 0.8 
condition. The fraction of it - r < 0.8 sources is ~10% for r < 21.5 and ~34% for 21.5 < r < 22. 

type (Covey et al. 2005), or stellar mass. The range of 
spectral types corresponding to each r — i bin is given in 
SpT column of table H 

For each color bin we select a volume limited sample 
given by: 

^0 = 10"^^""^ -+'pc, (14) 

Di^lO s 



will be fairly well preserved, making them a powerful 
tool for studying the Milky Way's stellar number density 
distribution. 

2.4. The Construction of the Density Maps 

The distance, estimated as described above, and the 
Galactic longitude and latitude, (I, b), fully determine the 
three-dimensional coordinates of each star in the sample. 
To better control the systematics, and study the depen- 
dence of density field on spectral type, we divide and 
map the sample in 19 bins in r — i color^^: 



'-+1 



pc. 



Here r„ 



riQ < r — I < rii 



(13) 



Typically, the width of the color bins, Ari = rii — rio, is 
Ari = 0.1 for bins redder than r — i = 0.7 and Ari = 0.05 
otherwise. The bin limits rig and rii for each color bin 
are given in the second column of tabled This color bin- 
ning is roughly equivalent to a selection by MK spectral 

To avoid excessive usage of parenthesis, we sometimes drop 
the minus sign when referring to the colors (e.g. g — r = gr or 
(r - i)i = rii). 



= 15 and rmax — 21.5 are adopted as bright 
and faint magnitude limits (SDSS detectors saturate at 
r ~ 14). In each color bin Di/Dq ~15, and for the full 
sample Di/Dq -300. 

We define the "Cartesian Galactocentric coordinate 
system" by the following set of coordinate transforma- 
tions: 



X = Rq - D cos{l) cos{b) 
Y = -Dsm{l)cos{b) 
Z^Dsm{b) 



(15) 
(16) 



where R(7) — 8 kpc is the adopted distance to the Galactic 
center ()Reidlll993D . 
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The choice of the coordinate system is motivated by 
the expectation of cylindrical symmetry around the axis 
of Galactic rotation Z, and mirror symmetry of Galac- 
tic properties with respect to the Galactic plane. Its 
{X, Y) origin is at the Galactic center, the X axis points 
towards the Earth, and the Z axis points towards the 
north Galactic pole. The Y = Z x X axis is defined so 
as to keep the system right handed. The X — Y plane is 
parallel to the plane of the Galaxy, and the Z — plane 
contains the Sun. The Galaxy rotates clockwise around 
the Z axis (the rotational velocity of the Sun is in the 
direction of the —Y axis). 

We bin the stars onto a three dimensional rectangular 
grid in these coordinates. The choice of grid pixel size 
is driven by compromise between two competing require- 
ments: keeping the Poisson noise in each pixel at a rea- 
sonable level, while simultaneously avoiding over-binning 
(and related information loss) in high-density regions of 
the maps. By manual trial-and-error of a few different 
pixel sizes, we come to a size (for each color bin) which 
satisfies both requirements. The adopted pixel sizes are 
listed in table [5] For bins with r — i > 0.3 the median 
number of stars per pixel is ^ 10, growing to ^ 30 for 
the bluest r ^ i bin. 

For each volume limited (rio,rii) color bin sample, 
this binning procedure results in a three-dimensional 
data cube, a map, of observed stars with each {X, Y, Z) 
pixel value equal to the number of stars observed in 
{X-dx/2,X+dx/2), {Y-dx/2,Y+dx/2), {Z-dx/2,Z+ 
dx/2) interval. 

Not all of the pixels in the maps have had their volume 
fully sampled by the SDSS survey. This is especially true 
near the edges of the survey volume, and at places where 
there are holes in the footprint of the survey (cf. figure 
[6|) . In order to convert the number of stars observed in a 
particular pixel (X, F, Z) to density, we must know the 
fraction of pixel volume that was actually sampled by the 
survey. Although simple in principle, the problem of ac- 
curately binning the surveyed volume becomes nontrivial 
due to overlap of observing runs, complicated geometry 
of the survey, and the large survey area. We solve it by 
shooting a dense, horizontal, rectangular grid of verti- 
cal {Xr = const, Yr = const) rays through the observed 
volume, with horizontal spacing of rays dxr being much 
smaller than the pixel size dx (typically, dxr/dx = 0.1). 
For each ray, we calculate the intervals in Z coordinate 
in which it intersects each imaging run ( "ray-run inter- 
sections"). Since imaging runs are bounded by simple 
geometric shapes (cones, spheres and planes), the ray- 
run intersection calculation can be done almost entirely 
analytically, with the only numerical part being the com- 
putation of roots of a 4'^ order polynomial. For each ray, 
the union of all ray-run intersections is the set of Z in- 
tervals {[Zo, Zi),[Z2, Z3),[Z4, Zi^), ...) at a given column 
{Xr, Yr) which were sampled by the survey. It is then a 
simple matter to bin such interval sets in Z direction, and 
assign their parts to pixels through which they passed. 
Then, by approximating that the ray sweeps a small but 
finite area dx^ , the survey volume swept by the ray con- 
tributing to pixel {X, Y, Z) is simply dx^ times the length 
of the ray interval(s) within the pixel. By densely cover- 
ing all of the {X, Y) plane with rays, we eventually sweep 



the complete volume of the survey and partition between 
all of the {X, Y, Z) pixels. This ray-tracing method is 
very general and can handle any survey geometry in any 
orientation, as long as the survey geometry can be rep- 
resented by a set of runs along great circles. Using this 
approach, we compute the volume observed within each 
pixel with an accuracy of one part in 10^. 

In summary, for each of the 19 r—i color bins, we finish 
with a three-dimensional map in which each (X, Y, Z) 
pixel holds the number of observed stars (N) and the 
observed volume {V). We estimate the number density 
in the pixel by simply dividing the two: 



p{X,Y,Z)^ 



N{X,Y,Z) 



(17) 



ViX,Y,Z)- 

with the error in density estimate due to shot noise being 



apiX,Y,Z) = 



^N{X,Y,Z) 



(18) 



ViX,Y,Z) 

For each pixel we also track additional auxiliary informa- 
tion (e.g. a list of all contributing SDSS runs), mainly 
for quality assurance and detailed a posteriori analysis. 

3. STELLAR NUMBER DENSITY MAPS 

In this Section we analyze the 19 stellar number density 
maps constructed as described above. The 0.10 < r — 
i < 1.40 color range spanned by our sample probes a 
large distance range - as the bin color is varied from the 
reddest to the bluest, the maps cover distances from as 
close as 100 pc traced by M dwarfs {r—i ^ 1.3), to 20 kpc 
traced by stars near the main sequence turnoff (r — i ~ 
0.1). We begin the analysis with a qualitative survey 
of the various map cross-sections, and then proceed to 
a quantitative description within the context of analytic 
models. 

3.1. The Number Density Maps in the R — Z Plane 
We first analyze the behavior of two-dimensional maps 
in the R — Z plane, where R = \fX^^^VY^ and Z are 
the galactocentric cylindrical coordinates. Assuming the 
Galaxy is circularly symmetric (we critically examine this 
assumption below), we construct these maps from the 
three-dimensional maps by taking a weighted mean of 
all the values for a given Z — R pixel (i.e. we average 
over the galactocentric polar angle (j) — arctanX). 

We show a subset of these maps in Fig. [TOl They 
bracket the analyzed r — i range; the remaining maps 
represent smooth interpolations of the displayed behav- 
ior. 

The bottom two panels in Fig. (TU] correspond to the 
reddest bins, and thus to the Solar neighborhood within 
~2 kpc. They show a striking simplicity in good agree- 
ment with a double exponential disk model: 



p(i?,Z) = p(i?o,0)e- 



R Z + Zp, 



expi-- 



H 



(19) 



. Here p is the number density of disk stars, Rq and Zq 
are the cylindrical coordinates of the Sun, and L and H 
are the exponential scale length and scale height, respec- 
tively. This model predicts that the isodensity contours 
have the linear form 



\Z + Za 



c 



H 



R, 



(20) 
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Fig. 10. — The stellar number density as a function of Galactic cylindrical coordinates R (distance from the axis of symmetry) and Z 
(distance from the plane of the Sun), for different r — i color bins, as marked in each panel. Each pixel value is the mean for all polar 
angles (f>. The density is shown on a natural log scale, and coded from blue to red (black pixels are regions without the data). Note that 
the distance scale greatly varies from the top left to the bottom right panel - the size of the the bottom right panel is roughly equal to 
the size of four pixels in the top left panel. Each white dotted rectangle denotes the bounding box of region containing the data on the 
subsequent panel. 
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behavior. 

We test whether the number density maps are circu- 
larly symmetric by examining isodensity contours on a 
cylindrical surface at i? = Rq kpc. Fig. [Tl] shows such 
projections for two color bins, where we plot the de- 
pendence of isodensity contours on galactocentric po- 
lar angle 0, and distance from the plane Z. In case 
of cylindrical symmetry, the contours would be horizon- 
tal. The top panel shows the isodensity contours for the 
1.0 < r — i < 1.1 color bin and is representative of all 
bins redder than r — i > 0.35 mag. The contours are 
horizontal, and the number density maps are indeed ap- 
proximately cylindrically symmetric. However, for bins 
r ~ i < 0.35 mag, detectable deviations from cylindri- 
cal symmetry do exist, especially at large distances from 
the Galactic plane (a few kpc and beyond). We show 
an example of this in the bottom panel, where there is 
a slight upturn of the isodensity contour at Z ~10,000 
and (f) 40°, indicating the presence of an overdensity. 
We will discuss such overdensities in more detail in the 
following section. 
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Fig. 11. — The azimuthal dependence of the number density for 
-R = Rq cyhnder around the Galactic center. The shaded region 
is the area covered by the SDSS survey, and the lines show con- 
stant density contours for two color bins (1.0 < r — i < 1.1 in the 
top panel and 0.10 < r — i < 0.15 in the bottom panel). The fact 
that isodensity contours are approximately horizontal supports the 
assumption that the stellar number density distribution is cylindri- 
cally symmetric around the Galactic center, and at the same time 
indicates that the assumed photometric parallax distribution is not 
grossly incorrect. Nevertheless, note that deviations from cylindri- 
cal symmetry do exist, e.g. at Z ~ 10 kpc and <f) 40° in the 
bottom panel. 

where C is an arbitrary constant, a behavior that is in 
good agreement with the data. 

As the bin color becomes bluer (the middle and top 
panels), and probed distances larger, the agreement with 
this simple model worsens. First, the isodensity contours 
become curved and it appears that the disk flares for 
i? > 14 kpc. Further, as we discuss below, the Z depen- 
dence deviates significantly from the single exponential 
given by eq. I19[ and additional components or a differ- 
ent functional form, are required to explain the observed 



3.2. The X ~Y Slices of the 3-dimensional Number 
Density Maps 

Instead of contracting the three-dimensional maps by 
taking the mean of all values for a given Z — R pixel, 
two-dimensional analysis can be based on simple cross- 
sections parallel to an appropriately chosen plane. A con- 
venient choice is to study the X — Y cross-sections that 
are parallel to the Galactic plane. A series of such pro- 
jections for the bluest color bin is shown in Figs. [T^HHl 
Their outlines are determined by the data availability. In 
particular, the gap between the two largest data regions 
will be eventually filled in as more SDSS imaging data 
becomes available^^. 

An unexpected large overdensity feature is easily dis- 
cernible in five of the six panels in Fig. [T^l In all stan- 
dard Galaxy models, the stellar density in the upper half 
{Y > 0) should mirror the bottom half {Y < 0), and in 
most models density depends only on the distance from 
the center of the Galaxy (each annulus enclosed by two 
successive circles should have roughly the same color). In 
contrast, the observed density map, with a strong local 
maximum offset from the center, is markedly different 
from these model predictions. This is the same feature 
that is responsible for the structure visible at Z ~10 kpc 
and R ^5 kpc in the top left panel in Fig. [TOl and for 
the upturn of the isodensity contour at Z ^10,000 and 
(j) ^ 40° in the bottom panel in Fig. [TT] We discuss this 
remarkable feature in more detail in Section [5l 

The top three panels (Z=3-5 kpc) in Fig. [13] clearly 
show another local overdensity at R ^ 16 kpc and 
Y ~ 0. This i s the "Monoceros Stream" discovered by 
iNewberg et al.l ()2002l ) using a subset of the data analyzed 
here (this overdensity is also discernible in the top left 
panel in Fig. [10] at i? ~ 16 kpc and Z - 3 kpc). The 
maps discussed here suggest that the stream is well lo- 
calized in the radial direction with a width of '--^ 3 kpc. 
This well-defined width rules out the hypothesis that this 
overdensity is due to disk flaring. 

An alternative hypothesis, that of a "ring" around the 

^■^ This region of the sky has already been imaged, and will be a 
part of SDSS Data Release 6 projected to be released in July 2007. 
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Fig. 12. — The stellar number density for the same color bin as in the top left panel in Fig. 1101 (0.10 < r — i < 0.15), shown here in slices 
parallel to the Galactic plane, as a function of the distance from the plane. The distance from the plane varies from 17.5 kpc (top left) to 
6 kpc (bottom right), in steps of 2 and 2.5 kpc. The circles visualize presumed axial symmetry of the Galaxy, and the origin marks the 
location of the Galactic center (the Sun is at X = 8, y = kpc). Note the strong asymmetry with respect to the Y = line. 
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Fig. 13. — Analogous to Fig. 1121 except that three symmetric slices at Z=3, 4 and 5 kpc above and below the plane are shown. The 
color stretch in panels for Z=3, 4 and 5 kpc is optimized to bring out the Monoceros overdensity at i? ~ 16 kpc and y ~ 0. 
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Fig. 14. — Analogous to Fig. 1121 except that here three symmetric sHces at Z=300, 600 and 900 pc above and below the plane are 
shown, for the 1.00 < r — i < 1.10 color bin. Note that at these distance scales there is no obvious discernible substructure in the density 
distribution. 
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Galaxy, was proposed by llbata et al.l (l2003l) , but put 
question b y observations oflRqcha-Pinto et aLl ()2003D . In 
particular, iRocha-Pinto et al.l analyzed the distribution 
of 2MASS M giants in the Monoceros feature and con- 
cluded its morphology was inconsistent with a homoge- 
neously dense ring surrounding the Milky Way. Instead, 
a more likely explanation is a merging dwarf galaxy 
with tidal arms. The inhomogeneity of the stream ap- 
parent in top three panels of Fig. [131 as well as i? = 
co nst, projections of these maps and a theoretical study 
by iPenarrubia et al.l (|2005. ) , support this conclusions as 
well. 

Closer to the plane, at distances of less than about 1 
kpc, the number density maps become smoother and less 
asymmetric, with deyiations from a simple exponential 
model giyen by eq. [H] not exceeding 30-40% (measured 
upper limit). This is true of all color bins for which 
region closer than 2 kpc is well sampled, and is shown 
in Fig. [T3] for 1.0<r — i<l.l color bin. 

3.3. Overall Distribution of Stellar Number Density 

Traditionally, the stellar distribution of the Milky Way 
has been decomposed into several components: the thin 
and thick disks, the central bulge, and a much more ex- 
tended and tenuous halo. While it is clear from the pre- 
ceding discussion that there are a number of overdensities 
that complicate this simple model, the dynamic range of 
the number density variation in the Galaxy (orders of 
magnitude) is large compared to the local density excess 
due to those features (a factor of few). Hence, it should 
still be possible to capture the overall density variation 
using analytic models. 

Before attempting a full complex multi-parameter fits 
to the overall number density distribution, we first per- 
form a simple qualitative exploration of density variation 
in the radial (R) and vertical (Z) directions. This type of 
analysis serves as a starting point to understand of what 
types of models are at all compatible with the data, and 
to obtain reasonable initial values of model parameters 
for global multi-parameter fits (Section [ 



3.3.1. The Z-dependence of the Number Density 

Fig. [Tni shows the stellar number density for several 
color bins as a function of the distance Z from the plane 
of the Sun at i? = Rq. The behavior for red bins, which 
probe the heights from 50 pc to ~2 kpc, is shown in the 
top panel. They all appear to be well fit by an exponen- 
tial profile^"' with a scale height of ^ 270 pc^^. While 
the best-fit value of this scale height is uncertain up to 
10-20%, it is encouraging that the same value applies to 
all the bins. This indicates that the slope of the adopted 
photometric parallax relation is not greatly incorrect at 
the red end. 



Motivated by theoretical reasoning (e.g.. iBinnev &: Tremain^ 
119871) . sometimes the sech^ function is used instead of exponen- 
tial dependence. However, the exponential provides a significantly 
better description of the data than sech-^. For example, the ex- 
ponential distribution is a good fit all the way towards the plane 
to 1/6 or so of the scale height, where the sech'^ function would 
exhibit significant curvature in the ln(p) vs. Z plot. 

Note that this is just an initial estimate for the scale height, 
based on a single effecti ve li ne of sight (SGP - NGP) and limited 
Z coverage. In Section 14.21 we will derive the values of Galactic 
model parameters using the entire dataset. 
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Fig. 15. — The vertical (Z) distribution of SDSS stellar counts 
for R = 8 kpc, and different r — i color bins, as marked. The lines 
are exponential models fitted to the points. The dashed lines in 
the top panel correspond to a fit with a single, exponential disk 
having a 270 pc scale height. The vertical dot-dashed line marks 
the position of the density maximum, and implies a Solar offset 
from the Galactic plane of ~ 20 pc. The dashed line in the middle 
panel correspond to a sum of two disks with scale heights of 270 
pc and 1200 pc, and a relative normalization of 0.04 (the "thin" 
and the "thick" disks). The dot-dashed line is the contribution of 
the 1200 pc disk. Note that its contribution becomes important 
for \Z\ > 1000 pc. The dashed line in the bottom panel (closely 
following the data points) corresponds to a sum of two disks (with 
scale heights of 260 pc and 1000 pc, and the relative normalization 
of 0.06), and a power-law spherical halo with power-law index of 
2, and a relative normalization with respect to the 260 pc disk of 
4.5x10"^. The dashed line is the contribution of the 260 pc disk, 
the dot-dashed line is the contribution of the 1000 pc disk, and 
the halo contribution is shown by the dotted line. Note that both 
the disk and halo models shown here are just the initial estimates 
of model parameters, based solely on this Z crossection. As we 
discuss in Section 14.3.91 these are not the only combinations of 
model parameters fitting the data, and the true model param eter s 
fitting all of the data are in fact substantially different (Table [TOt . 
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The extrapolations of the best exponential fits for 
Z < and Z > to small values of \Z\ cross at 
Z ~ —25 pc. This is the well-known Solar offset from 
the Galact ic plane towards the north Galactic pole (e.g. 
iReidl llOOBV which is here determined essentially directly 
using a few orders of magnitude greater number of stars 
(several hundred thousand) than in previous work. 

By selecting bluer bins, the Z dependence of the num- 
ber density can be studied beyond 1 kpc, as illustrated 
in the middle panel. At these distances the number den- 
sity clearly deviates from a single exponential disk model. 
The excess of stars at distances beyond 1 kpc, compared 
to this model, is usually interpreted as evidence of an- 
other disk component, the thick disk. Indeed, the data 
shown in the middle panel in Fig. [15] can be modelled 
using a double-exponential profile. 

The need for yet another, presumably halo, compo- 
nent, is discernible in the bottom panel in Fig. I15[ which 
shows the number density for the bluest color bin. The 
data show that beyond 3-4 kpc even the thick disk com- 
ponent underpredicts the observed counts. The observa- 
tions can be explained by adding a power-law halo com- 
ponent, such as described by eg. [Ml 

3.3.2. The R- dependence of the Number Density 

We examine the dependence of number density on 
the (cylindrical) distance from the Galactic center in 
Figs. [TBI [T7I and [T5] Each figure shows the number den- 
sity as a function of R for a given r — i color bin at differ- 
ent heights above the Galactic plane. For red bins, which 
probe the Solar neighborhood within ^2 kpc, the den- 
sity profiles are approximately exponential (i.e. straight 
lines in \n{p) vs. R plot, see Fig. [TB]) . The exponential 
scale length seems to increase with the distance from the 
Galactic plane, or alternatively, requires the introduction 
of an additional exponential dependence with a different 
scale. Due to the small baseline this variation or the 
scale lengths are not strongly constrained with plausible 
values around L 3.5 kpc and an uncertainty of at least 
30%. 

At distances from the Galactic plane exceeding 1-2 
kpc, the exponential radial dependence becomes a fairly 
poor fit to the observed density distribution (Fig. [T7]) . 
The main source of discrepancy are several overdensi- 
ties noted in Section 13.21 In particular, the Monoceros 
stream is prominent at Z ~2-8 kpc, especially when the 
density profiles are extracted only for |F| < 1 kpc slice 
(Fig.[T8|). 

4. GALACTIC MODEL 

The qualitative exploration of the number density 
maps in the preceding section, as well as the analysis 
of the density variation in the radial R and vertical Z di- 
rections, suggest that the gross behavior can be captured 
by analytic models. These typically model the number 
density distribution with two exponential disks, and a 
power-law (or de Vaucouleurs spheroid ) elliptical halo. 

Following earli er work (e.g. iMajewskH Il993l . 
ISiegel et al.ll2002l [Chen et all [200 If) , we decompose the 
overall number density into the sum of disk and halo 
contributions 

piR,Z)^ PDiR,Z) + pH{R,Z). (21) 
We ignore the bulge contribution because the maps an- 
alyzed here only cover regions more than 3-4 kpc from 
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Fig. 16. — The radial distribution of SDSS stellar counts for 
different r — i color bins, and at different heights above the plane, 
as marked in each panel (pc). The two dashed lines show the 
exponential radial dependence of density for scale lengths of 3000 
and 5000 pc (with arbitrary normalization). 

the Galactic center, where the bulge contribution is neg- 
ligible compared to the disk and halo contributions (for 
plausible bulge parameters determined using IRA S dat a 
for asymptotic g iant stars, see e.g. Jackson et al.ll2002[ ). 
Following iBahcall fc Soneiral ( 1980f) and 



iGilmore fc ReidI ( 1983f ). we further decompose the 
disk into a sum of two exponential components (the 
"thin" and the "thick" disk), allowing for different scale 
lengths and heights of each component: 



Pd{R,Z) 
where 



Pd{R, Z- ii, i/i) + !pd{R. Z- L2,H2) (22) 



Pd{R.Z;L,H)= pD{RQ,0)e- 



exp 1 



R Z + Zr. 



H 



(23) 

Here Hi, H-i and Li and L2 are the scale heights and 
lengths for the thin and thick disk, respectively, / is the 
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Fig. 17. — Analogous to Fig. 1161 except for bluer color bins, 
which probe larger distances. 

thick disk normalization relative to the thin disk at (R = 
Rq, Z = 0), and Zq is the Solar offset from the Galactic 
plane. From previous work, typical best-fit values are 
Hi ^300 PC H2 ~l-2 k pc, / -1-10% and Zq ~10-50 pc 
(e.g. iSiegel et al.ll2002l . table 1). 

We also briefly explored models where thin and thick 
disk had the same scale length that was allowed to vary 
linearly with distance distance from the Galactic plane 
(L = Lo + kZ), but found these to be unnecessary as the 
two-disk formalism was able to adequately capture the 
behavior of the data. 

We model the halo as a two-axial power-law ellipsoid^^ 



PHiR,Z)^pDiRQ,0)fH 



Rr. 



yi?2 + (Z/g^)2 



(24) 

The parameter q}j controls the halo ellipticity, with the 

^® For the halo component, Z + Zq Z is a very good approx- 
imation. 
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Fig. 18. — The radial distribution of SDSS stellar counts for 
0.10 < r — i < 0.15 color bin, with the data restricted to \y\ < 1 
kpc. The selected heights are, from top to bottom, (2,3,4), (4,5,6) 
and (6,8,10) kpc. The Monoceros stream at is easily visible as local 
maxima at R = 16 — 17 kpc, and the Virgo overdensity as the wide 
bump at i? ~ 6 kpc. 

ellipsoid described by axes a — b and c ~ qn a- For 
qn < ^ the halo is oblate, that is, "squashed" in the 
same sense as the disk. The halo normalization relative 
to the thin disk at [R = Rq , Z = 0) is specified by fu ■ 
From previous work, typical best-fit values are uh ~2.5- 
3.0, fH - 10"^ and qn 0.5 - 1. 

4.1. Dataset Preparation 

The fitting of models described by eqs. [2TH24l will be af- 
fected by overdensities identified in Sect ion [5T^ and other, 
smaller overdensities that may be harder to see at first. 
If unaccounted for, such overdensities will almost cer- 
tainly bias the best-fit model parameters. In general, as 
we discuss later in Section 14.3.71 their effect is to artifi- 
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Fig. 19. — The regions with large overdensities excluded from 
Galactic model fits. The pixels within the rectangle in the top 
panel are excluded to avoid contamination by the Virgo ovcrden- 
sity (Section [5J. The pixels enclosed by the two rectangles in the 
bottom panel, centered at i? ~ 18 kpc, exclude the Monoceros 
stream. 



daily increase the scale heights of the disks, in order to 
compensate for the localized density excesses away from 
the plane. We therefore exclude from the dataset the re- 
gions where there are obvious localized deviations from 
smooth background profile^^. The excluded regions are 
shown in Figure [T9l 

We exclude the newly found large overdensity dis- 
cernible in Fig. [T2] (the "Virgo overdensity" ) by masking 
the pixels that simultaneously satisfy: 

- 5 < X7kpc<25 

y'>-4kpc 
(X - 8 kpc)2 + y2 + > (2.5 kpc)2 

Note that we are excluding overdensities, but not underden- 
sities, as there are physical reasons to expect the Galaxy to have a 
smooth distribution with overdense regions (e.g., due to mergers, 
clusters, etc.). 



X' 



cos 30° -sin 30° 
sin 30° cos 30° 



The third condition excludes from the cut pixels closer 
than 2.5 kpc to the Sun, which are uncontaminated by 
the overdensity. The excluded region is is shown on Fig- 
ure [19] bounded by the rectangle in the top panel. 

The Monoceros stream is located at an approximately 
constant galactocentric radius. We exclude it by masking 
out all pixels that satisfy either of the two conditions: 

14kpc <R< 24kpcA0 <Z < 7kpc 
16kpc <R< 24kpcA7 <Z < lOkpc 

. These correspond to the region bounded by two white 
rectangles in the bottom panel of Fig. [191 

After the removal of Virgo and Monoceros regions, the 
initial fit for bins redder than r — i = 1.0 resulted in 
measured thin and thick scale heights of Hi ^ 280 and 
H2 ^ 1200. The residuals of this fit showed clear signa- 
tures of at least two more major overdensities (~ 40% 
above background), one near (i?, Z) ^ (6.5, 1.5) kpc and 
the other near (i?, Z) ^ (9, 1) kpc. We therefore went 
back and further excluded the pixels satisfying: 

- 90° < arctan(^-=-^5^5^) < 18° A Z > 
it — 8.6kpc 

R < 7.5kpc /\ Z >0 

The remaining pixels are averaged over the galactocen- 
tric polar angle 0, to produce the equivalent of (i?, Z) 
maps shown in fig. 1101 We additionally imposed a cut 
on Galactic latitude, excluding all pixels with b < 20° 
to remove the stars observed close to the Galactic disk. 
This excludes stars that may have been overcorrected for 
extinction (Section I2.3.2p . and stars detected in imag- 
ing runs crossing the Galactic plane where the efficiency 
of SDSS photometric pipeline drops due to extremely 
crowded fields. Other, less significant, r — i bin-specific 
cuts have also been applied, for example the exclusion of 
\Z\ > 2500 pc stars in r — « > 1.0 bins to avoid contami- 
nation by halo stars. 

We show all 19 "cleaned up" maps in figure [20l The 
contours denote the locations of constant density. The 
gray areas show the regions with available SDSS data. 
Compared to Fig. [TOl the constant density contours are 
much more regular, and the effect of the Virgo overden- 
sity is largely suppressed. The regularity of the density 
distribution is particularly striking for redder bins (e.g., 
for r - i > 0.7). In the bluest bin (0.10 < r - i < 0.15), 
there is a detectable departure from a smooth profile in 
the top left part of the sampled region. This is the area 
of the (i?, Z) plane where the pixels that are sampled 
far apart in {X, Y, Z) space map onto adjacent pixels in 
{R, Z) space. Either deviations from axial symmetry or 
small errors in photometric parallax relation (perhaps 
due to localized metallicity variations) can lead to devi- 
ations of this kind. Unfortunately, which one of the two 
it is, is impossible to disentangle with the data at hand. 

4.2. Model Fit 
4.2.1. Fitting Algorithm 
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Fig. 20. — "Cleaned up" (R, Z) maps of the Galaxy, analogous to figure [TOl but with pixels in obvious overdensities (fig. I19I I excluded 
from azimuthal averaging. We show the maps for all 19 color bins, with the bluest bin in the top left corner and the reddest bin in the 
bottom right. The contours are the lines of constant density, spaced at constant logarithmic intervals. 



The model fitting algorithm is based on the standard 
Leven berg-Marquardt n onlinear minimization algo- 
rithm (jPress et al.lll992^ . with a multi-step iterative out- 
lier rejection. 

The goal of the iterative outlier rejection procedure is 
to automatically and gradually remove pixels contami- 
nated by unidentified overdensities, single pixels or small 
groups of pixels with large deviations (such as those due 
to localized star clusters, or simply due to instrumental 
errors, usually near the edges of the volume), and allow 



the fitter to "settle" towards the true model even if the 
initial fit is extremely biased by a few high-tr outliers. 

The outlier rejection works as follows: after initial fit 
is made, the residuals are examined for outliers from the 
model higher than a given number of standard devia- 
tions, CTi. Outlying data points are excluded, the model 
is refitted, and all data points are retested with the new 
fit for deviations greater than (72, where (72 < fi- The 
procedure is repeated with (73 < a2, etc. The removal 
of outliers continues until the last step, where outliers 
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higher than ctjv are excluded, and the final model refit- 
ted. The parameters obtained in the last step are the 
best fit model parameters. 

The <7i sequence used for outlier rejection is ct; = 
{50, 40, 30, 20, 10, 5}. This slowly decreasing sequence al- 
lows the fitter to start with rejecting the extreme outliers 
(which themselves bias the initial fit), and then (with the 
model now refitted without these outliers, and therefore 
closer to the true solution) gradually remove outliers of 
smaller and smaller significance and converge towards a 
solution which best describes the smooth background. 

4.2.2. A Measurement of Solar Offset 

We begin the modeling by fitting a single exponential 
disk to the three reddest color bins to find the value 
of the Solar offset Zq. To avoid contamination by the 
thick disk, we only use pixels with \Z\ < 300 pc, and 
to avoid effects of overestimated interstellar extinction 
correction for the nearest stars (Section l2.3.2p . we further 
exclude pixels with \Z\ < 100. We further exclude all 
pixels outside of 7600 < R < 8400 pc range, to avoid 
contamination by clumpy substructure. 

We obtain 

,bright — 

(25±5)pc (25) 
^o,fai„t = (24±5)pc (26) 

for the Solar offset, where the 2'Q^bright is the offset ob- 
tained using the bright photometric parallax relation, 
and ^Q^faint using the faint. The quoted uncertainty is 
determined by simply assuming a 20% systematic uncer- 
tainty in the adopted distance scale, and does not imply 
a Gaussian error distribution (the formal random fitting 
error is smaller than 1 pc). 

Our value of the Solar offset agrees favorably with re- 
cent independent measurements {Zq = (27.5 ± 6) pc , 
iChen et all [19991 : Zq = (27 ± 4) pc, iChen et all [200ll : 
(24.2 ± 1.7) pc obtained from trigonometric Hipparcos 
data bv lMafz-Apellani9l2001f ). We keep the value of the 
Solar offset fixed in all subsequent model fits. 

4.2.3. Disk Fits 

We utilize the R—Z density maps of the four r—i > 1.0 
bins to fit the double-exponential disk model. These 
color bins sample the thin and thick disk, with a neg- 
ligible halo contribution (less than ~ 1% for plausible 
halo models). Furthermore, the photometric relations in 
this range of colors are calibrated to metallicities of disk 
dwarfs, thus making these bins optimal for the measure- 
ment of disk model parameters. 

We simultaneously fit all double-exponential disk 
model parameters (p. Hi, Li, /, H2, L2) to the data, for 
both bright and faint photometric parallax relations. To 
avoid contamination by the halo, we only use the pixels 
with \Z\ < 2500 pc. To avoid effects of overestimated in- 
terstellar extinction correction for the nearest stars (Sec- 
tion [2?3l2]) , we further exclude pixels with \Z\ < 100. 

We jointly fit the data from all four color bins, and sep- 
arately for each bin. In the former, "joint fit" case, only 
the densities p{Rq , 0) are allowed to vary between the 
bins, while the scale lengths, heights and thick-to-thin 
disk normalization / are constrained to be the same for 
stars in each bin. As the color bins under consideration 
sample stars of very similar mass, age and metallicity, we 



expect the same density profile in all bins^*. The best 
fit parameters for the joint fit to r — z > 1.0 bins are 
given in top row of Tables [3] and [H calculated assuming 
the bright (Equation [2]) and faint (Equation [T]) photo- 
metric parallax relation, respectively. Two-dimensional 
cross sections of reduced hyper-surface around best- 
fit values are shown in Figure [2T] (for bright relation only 
- analogous crossections obtained with the faint relation 
look qualitatively the same). 

In case of separate fits, all parameters are fitted inde- 
pendently for each color bin. Their variation between 
color bins serves as a consistency check and a way to as- 
sess the degeneracies, significance and uniqueness of the 
best fit values. The best-fit values are shown in the top 
four rows of Table [6] (bright photometric parallax rela- 
tion) and the top five^^ rows of Table [7] (faint relation) . 

In all cases we are able to obtain good model fits, with 
reduced in the range from 1.3 to 1.7. The best-fit solu- 
tions are mutually consistent. In particular, the thin disk 
scale height is well constrained to Hi = 250 pc (bright) 
and Hi — 230 — 240 (faint), as are the values of p(Rq, 0) 
which give the same results in individual and joint fits at 
the 5% level. 

The thick-to-thin disk density normalization is 10%, 
with / = 0.10 - 0.13 (bright) and / = 0.10 - 0.14 
(faint). The thick disk scale length solutions are in 
H2 = 750-900 pc (bright) and H2 = 660-900 pc (faint) 
range. Thick disk normalization and scale heights ap- 
pear less well constrained; however, note that the two are 
fairly strongly correlated (/ vs H2 panel in Figure [21]). 
As an increase in density normalization leads to a de- 
crease in disk scale height and vice versa with no appre- 
ciable effect on , any two models with so correlated dif- 
ferences of scale height and normalization of up to 20% 
to 30% are practically indistiguishable. This interplay 
between p and H2 is seen in Tables |6| and [3 most ex- 
tremely for 1.1 < r — i < 1.2 bin (Table [71 third row). 
With this in mind, the fits are still consistent with a sin- 
gle thick disk scale height H2 and density normalization 
/ describing the stellar number density distribution in 
all r — i > 1.0 color bins. 

Constraints on disk scale lengths are weaker, with the 
goodness of fit and the values of other parameters be- 
ing relatively insensitive on the exact values of Li and 
L2 (Figure [5T1 first two columns). This is mostly due 
to a short observation baseline in the radial (i?) di- 
rection. The best fit parameters lie in the range of 
Li = 1600 - 2400 pc, L2 = 3200 - 6000 pc (bright) 
and Li = 1600 - 3000 pc, L2 = 3000 - 6000 pc (faint 
parallax relation). Note that the two are anticorrelated 
(Figure [21] top left panel), and combinations of low Li 
and high L2, or vice versa can easily describe the same 
density field with similar values of reduced (the be- 
havior seen in Tables [6| and [7|) . The disk scale length 
fits in individual color bins are also consistent with there 
being a single pair of scale lengths Li and L2 applicable 
to all color bins. 

Note also that being O.lmag wide, with a typical magnitude 
errors of ct,. > 0.02m ag th e adjacent bins are not independent. The 
histograms in Figure l46l illustrate this well. 

The fit for 0.9<r — i<1.0 bin when using the faint photo- 
metric relation and including a halo component (see Section l4.2.5ll . 
failed to converge to physically reasonable value. We have therefore 
fitted this bin with disk components only. 
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Fig. 21. — Two-dimensional cross sections of of reduced hyper-surface around best-fit values for 1.0 < r — i < 1.4 data (TableO first 
row). The fit was obtained assuming the "bright" photometric paralax relation (Equation [2J. Analogous cross sections for fits obtained 
assuming Equation [T] (Table |4] first row) show qualitatively same features. The innermost contour is at 1.1 X Xmin l^vel, while the rest are 
logarithmically spaced in steps of 0.5 dex, starting at logx^ = 0.5. 



TABLE 3 

Best Fit Values (.Joint Fits, Bright Paralax Relation) 





Bin 




p(i?Q,0) 


Li 


Hi 


/ 


L2 


H2 


fH 


1.61 


1.3 < r - 
1.2 < r - 
1.1 < r - 
1.0 < r - 


i < 1.4 
i < 1.3 
i < 1.2 
i < 1.1 


0.0058 
0.0054 
0.0046 
0.0038 


2150 


245 


0.13 


3261 


743 




1.70 


0.9 < r - 
0.8 < r - 
0.7 < r - 
0.65 < r - 


i < 1.0 
i < 0.9 
i < 0.8 
i < 0.7 


0.0032 
0.0027 
0.0024 
0.0011 


2862 


251 


0.12 


3939 


647 


0.00507 



Note. — Best fit values of Galactic Model parameters derived assuming the "bright" 
photometric paralax relation (Equation[2j. The fit to 0.65 < r — i < 1.0 bins (bottom 
row) includes the halo component. Its shape was kept fixed (Table [5] top row) and 
only the normalization fn was allowed to vary. 
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TABLE 4 

Best Fit Values (Joint Fits, Faint Paralax Relation) 



Bin 



p{Rq,0) 



Hi f 



L2 H2 



1.59 



2.04 



1.3 < r - 


i < 1.4 


0.0064 


1.2 < r - 


i < 1.3 


0.0063 


1.1 < r - 


i < 1.2 


0.0056 


1.0 < r - 


i < 1.1 


0.0047 


0.9 < r - 


i < 1.0 


0.0043 


0.8 < r - 


i < 0.9 


0.0036 


0.7 < r - 


i < 0.8 


0.0032 


0.65 < r - 


i < 0.7 


0.0015 



2037 229 0.14 3011 662 



2620 225 0.12 3342 583 0.00474 



Note. — Best fit values of Galactic Model parameters derived assuming the "faint" 
photometric paralax relation (Equation[^. The fit to 0.65 < r — i < 1.0 bins (bottom 
row) includes the halo component. Its shape was kept fixed (Table[5] bottom row) and 
only the normalization fn was allowed to vary. 
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Reduced ^ = 5.048 
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Fig. 22. — Reduced X'^ surface of halo parameters nj] and 
qn around the best-fit values (Table [S] first row). The innermost 
contour is at 1.1 X Xmin l^vel, while the rest are logarithmically 
spaced in steps of 0.5 dex, starting at logx'^ = 0.5. 



TABLE 5 
Halo Shape and Profile Fit 



Paralax Relation 








Bright 


3.05 


0.64 ± 0.01 


2.77 ± 0.03 


Faint 


2.48 


0.62 ±0.01 


2.78 ±0.03 


Note. — Best fit 


values 


of halo power 


law index nn 



and axis ratio qn — c/a, assuming the "bright" (top) and 
"faint" (bottom row) photometric paralax relation. 

4.2.4. Halo Fits 

For bluer color bins {r — i < 1.0) the probed distance 
range is larger, and the stellar halo component starts to 
appreciably contribute to the total density near the far 
edge of survey volume. As seen in the middle and bottom 
panel of FigurefTSl the disk-only solution becomes visibly 
unsatisfactory at .Z > 4000 kpc. Also, the reduced 
values of disk-only models begin to climb to higher than 
a few once we attempt to fit them to data in r — i < 1.0 
bins. 

Before we move on to adding and fitting the halo com- 
ponent, there are a few significant caveats that must be 
discussed, understood and taken into account. Firstly, 
the presence of dumpiness and merger debris in the halo, 
if unaccounted for, will almost certainly bias (an make 
difficult, or even impossible to determine) the model pa- 
rameters. An initial survey of the density field (Sec- 
tion l2.3p . the identification, and careful removal of identi- 
fied overdensities fSection l4.1[) are essential for obtaining 
a reasonable fit. 

Secondly, the photometric parallax relations 
(eqs. [T] and [2]) do not explicitly depend on stellar 
metallicity. Implicitly, as discussed in Section 12. 2. H 
they take metallicity into account by virtue of being 
calibrated to disk M-dwarfs on the red, and metal-poor 
halo stars at the blue end. This makes them correct for 
low metallicity stars ([Fe/H] < -1.5) near r — i ^ 0.15, 
and high metalhcity ([Fe/H] > -0.5) at r - i > 0.7. 
They are therefore appropriate for the study of halo 



shape and parameters only at the blue end, and disk 
shape and parameters only on the red end. Conversely, 
they are inappropriate for the study of disk shape and 
parameters at the blue, or halo shape and parameters 
at the red end. For the same reason, it is difficult 
to simultaneously fit the halo and the disk in the 
intermediate r — i bins, as the application of photometric 
parallax relation inappropriate for the low metallicity 
halo induces distortions of halo shape in the density 
maps. 

Therefore, to measure the shape of the halo, we only se- 
lect the data points from the three bluest, 0.1 < r — i < 
0.25 bins, and only in regions of (i?, Z) plane where a 
fiducial qn — 0.5, uh — 2.5, fn = 0.001 halo model 
predicts the fraction of disk stars to be less than 5%. 
This allows us to fit for the power law index , and the 
axis ratio qn of the halo. Because we explicitly excluded 
the disk, we cannot fit for the halo-to-thin-disk normal- 
ization Jh (but see Section 14.2.51 later in the text for a 
workaround) . 

The best fit parameters obtained for the halo are shown 
in Table [5] for both the bright and faint photometric re- 
lation, and the reduced surface for the fit is shown in 
Figure [H] (bright relation only - the surface looks qual- 
itatively the same for the fit made assuming the faint 
relation). 

The fits are significantly poorer than for the disks, with 
reduced = 2 — 3. Formal best-fit halo parameters 
are nn = 2.8, qn = 0.64, but given the relatively high 
and shallow minimum, and the shape of the surfaces 
in Figure [511 it is better to think of the fit results as 
constraining the parameters to a range of values - the 
power law index to uh = 2.5 — 3, and the oblateness 
parameter qn = 0.5 — 0.8. 

Fig. [Ml shows residual maps for the bluest color bin 
and for four different halo models, with the thin and 
thick disk parameters kept fixed at values determined 
using redder bins (Table [3]). Individual panels illustrate 
the changes in residuals when the halo power law index 
is varied while keeping the axis ratio fixed (top row), and 
when the ellipticity of the halo is changed, from oblate 
to spherical while keeping the power law index hh fixed 
(bottom row). The Monoceros and Virgo overdensities, 
and the overdensity at R ^--^6.5 kpc and Z ~ 1.5 kpc, are 
clearly evident, but their detailed properties depend sig- 
nificantly on the particular halo model subtracted from 
the data. 

We further find that a power-law halo model always 
over- or underestimates the stellar counts in the far 
outer halo (Figure [53]), suggesting the use of a dif- 
ferent profile may be more appropriate and consistent 
with "dual-halo" p rofiles favo red by (among o t hers) 
ISommer-Larsen fc Zhen (1990l):lAnen et al.l (fT99ll): IZinnl 
(199^ : ICarnev et all (119961); IChiba fc BeersI j2000l ) and 
more recently discussed by iSiegel et al.l (|2002[ ). 

However, no matter what the exact shape of the profile 
or the power law index is, only significantly oblate halos 
provide good fits to the data (compare the bottom right 
to other panels in Fig. [53]). Specifically, given the reduced 
surface in Figure [22l spherical or prolate halo can be 
ruled out, and this remains to be the case irrespective of 
the details of the photometric parallax relation'^°. 

^'^ Aspherical halos could be artificially favored by the anal- 
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Fig. 23. — Data-model residuals, normalized to the model, for color bin 0.10 < r — i < 0.15, using for four different halo models. All 
four models have identical thin and thick disk parameters, and only the halo parameters are varied. Panels in the top row illustrate the 
changes in residuals when the halo power law index nu is varied while keeping the axis ratio fixed. Panels of the bottom row illustrate the 
effects of axis ratio qn change, while keeping the power law index constant. While njy is not strongly constrained, the data strongly favor 
an oblate halo. 



4.2.5. Simultaneous Disk and Halo Fits 

Keeping the best fit values of halo shape parameters 
qn and uh constant, we next attempt to simultaneously 
fit the thin and thick disk parameters and the halo nor- 
malization, fn, in four 0.65 < r — i < 1.0 bins. 

These bins encompass regions of (i?, Z) space where 
the stellar number density due to the halo is not neg- 
ligible and has to be taken into account. Simultaneous 
fits of both the disk and all halo parameters are still 
unfeasible, both because halo stars still make up only 
a small fraction of the total number density, and due to 
poor applicability of the disk-calibrated photometric par- 
allax relations in this r — i range to low-metallicity halo 
stars. However, knowing the halo shape from the blue, 

ysis, as a way to parametrize away any existing halo inhomogene- 
ity. However, given the analysis of residuals in Section 14.3.81 we 
consider this to be a very unlikely explanation of the measured 
oblateness. 



low-metallicity calibrated bins, we may keep qn and nn 
fixed and fit for the halo-to-thin-disk normalization, fn- 
Given the uncertainty in its current knowledge, thusly 
obtained value of Jh is still of considerable interest de- 
spite the likely biases. 

We follow the procedure outlined in Section [4.2.3i and 
fit the data in 0.65 < r — i < 1.0 bins jointly for all and 
separately in each color bin, with both bright and faint 
photometric parallax relations. 

The results of the joint fits are given in bottom rows 
of Tables [3] and [H Results for individual bins are given 
in bottom rows of Tables [H] and [7] for the bright and faint 
photometric relation, respectively. 

We obtain satisfactory model fits, with reduced in 
1.4 to 2.0 range. As was the case for fits to r — z > 1.0 
bins, the best fit disk parameter values are consistent be- 
tween bins, and with the joint fit. The reduced surface 
cross-sections, shown in Figure [24l are qualitatively the 
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Fig. 24. — Two-dimensional cross sections of of reduced hyper-surface around best-fit values for 0.65 < r — i < 1.0 data (Table [3] 
second row). The fit was obtained assuming the "bright" photometric paralax relation (Equation [2J and includes the contribution of the 
halo. Analogous cross sections for fits obtained assuming Equation[l](Table|4l second row) show qualitatively same features. The innermost 



contour is at 1.1 X X^in while the rest are logarithmically spaced in steps of 0.5 dex, starting at 



: 0.5. 



TABLE 6 

Best Fit Values (Individual Fits, Bright Paralax Relation) 



Color bin 






Li 


m 


/ 


L2 


H2 


fH 


1.3 < r - 


i < 1.4 


1.34 


0.0062 


1590 


247 


0.09 


5989 


909 




1.2 < r - 


i < 1.3 


1.31 


0.0055 


1941 


252 


0.11 


5277 


796 




1.1 < r - 


i < 1.2 


1.58 


0.0049 


2220 


250 


0.09 


3571 


910 




1 < r - 


i < 1.1 


1.64 


0.0039 


2376 


250 


0.10 


3515 


828 




0.9 < r 


~i<l 


1.38 


0.0030 


3431 


248 


0.14 


2753 


602 


0.0063 


0.8 < r - 


i < 0.9 


1.48 


0.0028 


3100 


252 


0.10 


3382 


715 


0.0039 


0.7 < r - 


i < 0.8 


1.83 


0.0024 


3130 


255 


0.09 


3649 


747 


0.0037 


0.65 <r- 


i < 0.7 


1.69 


0.0011 


2566 


273 


0.05 


8565 


861 


0.0043 



Note. — Best fit values of Galactic Model parameters, fitted separately for each 
r — i bin assuming the "bright" photometric paralax relation (Equation [2j. In fits 
which include the halo component, the shape of the halo was kept fixed (Table[5] top 
row), and only the normalization fn was allowed to vary. 
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TABLE 7 

Best Fit Values (Individual Fits, Faint Paralax Relation) 



Color bin 










/ 


L2 


H2 


fH 


1.3 < r - i < 1.4 


1.32 


0.0064 


1599 


246 


0.09 


5800 


893 




1.2 < r - i < 1.3 


1.40 


0.0064 


1925 


242 


0.10 


4404 


799 




1.1 < r - i < 1.2 


1.56 


0.0056 


2397 


221 


0.17 


2707 


606 




1 < r - i < 1.1 


1.71 


0.0049 


2931 


236 


0.10 


2390 


760 




0.9 < r - i < 1 


1.62 


0.0043 


3290 


239 


0.07 


2385 


895 




0.8 < r - i < 0.9 


1.69 


0.0038 


2899 


231 


0.08 


2932 


759 


0.0021 


0.7 < r - i < 0.8 


2.59 


0.0034 


2536 


227 


0.09 


3345 


671 


0.0033 


0.65 < r - i < 0.7 


1.92 


0.0016 


2486 


241 


0.05 


6331 


768 


0.0039 



Note. — Best fit values of Galactic Model parameters, fitted separately for each 
r — i bin assuming the faint" photometric paralax relation (Equation^. In fits which 
include the halo component, the shape of the halo was kept fixed (Table [5] bottom 
row), and only the normalization fn was allowed to vary. 
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same as those in Figure [21] and the entire discussion of 
Section [4.2.31 about fit parameters and their interdepen- 
dencies applies here as well. 

Comparison of top and bottom rows in Tables [3| and [4| 
shows consistent results between r ~ i > 1.0 and 0.65 < 
r — i < 1.0 bins. In particular, the scale heights of the 
thin disk are the same, and the thick-to-thin disk nor- 
malization is the same to within 8 — 15%, still within 
fit uncertancies. The scale lengths are still poorly con- 
strained, and on average 10—30% larger than in disk-only 
fits. Given the poor constraint on scale lengths, it is dif- 
ficult to asses whether this effect is physical, or is it a 
fitting artifact due to the addition of stellar halo compo- 
nent. The scale height of the thick disk, H2 is ~ 14% 
smaller than in disk-only fits. This is likely due to the 
reassignment to the halo of a fraction of stellar number 
density previously assigned to the thick disk. 

For fn, the halo-to-thin-disk normalization at {R = 
8 kpc, Z = 0), the best fit values are in 0.3 — 0.6% range, 
with the best fit value for joint fits being fn = 0.5% both 
for the bright and faint parallax relation. In particular, 
note how insensitive fn is on the choice of photomet- 
ric parallax relation. In this region of r — i colors, the 
average difference between the bright and faint parallax 
relations is AMr — 0.25 mag; therefore even in case of 
uncertainties of ^ half a magnitude, the change in fn 
will be no greater than ~ 10 — 20%. 

4.3. Analysis 

The Galactic model parameters as fitted in the pre- 
ceding Section are biased'^^ by unrecognized stellar mul- 
tiplicity, finite dispersion of the photometric parallax re- 
lation and photometric errors. They are further made 
uncertain by possible systematics in calibration of the 
photometric parallax relations, and a simplified treat- 
ment of stellar metallicities. 

In this section, we analyze all of these (and a number of 
other) effects on a series of Monte Carlo generated mock 
catalogs, and derive the corrections for each of them. We 
also look at the resolved local overdensities found in the 
data, discuss the question of possible statistical signa- 
tures of further unresolved overdensities, and questions 
of uniqueness and degeneracy of our best-fit model. 

After deriving the bias correction factors, we close the 
Section by summarizing and writing out the final de- 
biased set of best fit SDSS Galactic model parameters, 
together with their assumed uncertainties. 

4.3.1. Monte-Carlo Generated Mock Catalogs 

To test the correctness of the data processing and 
fitting procedure and derive the correction factors for 
Malmquist bias, stellar multiplicity and uncertainties due 
to photometric parallax systematics, we developed a soft- 
ware package for generating realistic mock star cata- 
logs. These catalogs are fed to the same data-processing 
pipeline and fit in the same manner as the real data. 

The mock catalog generator, given an arbitrary Galac- 
tic model (which in our case is defined by eas. [2TH24l a lo- 
cal position-independent luminosity function, and binary 
fraction), generates a star catalog within an arbitrarily 
complex footprint on the sky. The code can also include 

Or "apparent" , in the terminology of lKroupa et al.1ll993l 



realistic magnitude-dependent photometric errors (Fig- 
ure [71 bottom panel) and the errors due to Gaussian 
dispersion (Tj\/^ around the photometric parallax mean, 
Mr{r-i). 

Using this code, we generate a series of mock cata- 
logs within the footprint of the SDSS data used in this 
study (Figure [Q using a fiducial model with parame- 
ters listed in the top r ow of Table [8} For t he luminos- 
ity function, we use the iKroupa et al.l (|1993f l luminosity 
function, transformed from (j){Mv) to (f>(Mr) and renor- 
malized to p{R = 8000, Z = 0) = 0.04 stars pc"^ mag~i 
in the 1.0 < r — i < 1.1 bin. As we will be making com- 
parisons between the simulation and r — i > 1.0 bins, we 
do not include the halo component {Jh ~ 0). 

For all tests described in the text to follow, we generate 
the stars in 0.7 < r — i < 1.6 color and 10 < r < 25 mag- 
nitude range, which is sufficient to include all stars that 
may possibly scatter into the survey flux (15 < r < 21.5) 
and disk color bins (1.0 < r — i < 1.4) limits, either 
due to photometric errors, uncertainty in the photomet- 
ric parallax relation, or an added binary companion. To 
transform from distance to magnitude, we use the bright 
photometric parallax relation (eq. [2|). 

4.3.2. Correctness of Data Processing and Fitting Pipeline 

We first test for the correctness of the data process- 
ing and fitting pipeline, by generating a "perfect" cat- 
alog. Stars in this catalog have no photometric errors 
added, and their magnitudes and colors are generated 
using eqs. [2l and [5l 

We fit this sample in the same manner as the real data 
in Section [4.2.31 The results are given in the second row 
of Table \8\ The fit recovers the original model parame- 
ters, with the primary source of error being the "cosmic 
variance" due to the finite number of stars in the catalog. 

This test confirms that fitting and data processing 
pipelines introduce no additional uncertainty to best-fit 
model parameters. It also illustrates the limits to which 
one can, in principle, determine the model parameters 
from our sample assuming a) that stars are distributed in 
a double-exponential disk and b) the three-dimensional 
location of each star is perfectly known. These limits are 
about 1 — 2%, significantly smaller than all other sources 
of error. 

4.3.3. Effects of Malmquist Bias 

We next test for the effects of photometric errors, and 
the errors due to the finite width of the photometric par- 
allax relation. We model the photometric errors as Gaus- 
sian, with a magnitude-dependent dispersion ar mea- 
sured from the data (Figure [3 bottom panel) . Median 
photometric errors range from cr,. = 0.02 on the bright 
to cr,. = 0.12 on the faint end. We assume the same 
dependence holds for g and i band as well. We model 
the finite width of the photometric parallax relation as 
a Gaussian = 0-3 dispersion around the mean of 
Mr{r — i). The two sources of effective photometric error 
add up in quadrature and act as a source of a Malmquist 
bias, with the photometric parallax relation dispersion 
giving the dominant effect (eq. [T2ll . 

The best fit parameters obtained from this sample are 
given in the third row of Table[5| The thin and thick disk 
scale heights are underestimated by ~ 5%. The density 
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TABLE 8 

Monte Carlo Generated Catalog Fits 



Simulation ^1 Hi f L-z H2 

True Model ■■■ 2500 240 0.10 3500 770 

Perfect Catalog 1.03 2581 ± 44 244 ± 2 0.094 ± 0.009 3543 ± 69 791 ± 18 

Photometric and Paralax Errors 0.95 2403 ± 40 230 ± 2 0.111 ± 0.010 3441 ± 57 725 ± 13 
25% binary fraction 0.97 2164 ± 39 206 ± 1 0.119 ±0.011 3199 ± 47 643 ± 9 

50% binary fraction 0.97 1986 ± 34 193 ± 1 0.115 ± 0.011 2991 ± 41 611 ± 7 

100% binary fraction 1.02 1889 ± 31 178 ± 1 0.104 ± 0.010 2641 ± 31 570 ± 6 

Note. — The true model parameters (top row), and the best fit values of model parameters recovered from 
a scries of Monte-Carlo generated catalogs. These test the correctness of data processing pipeline and the 
effects of cosmic variance ("Perfect Catalog"), the effects of photometric paralax dispersion and photometric 
errors ("Photometric and Paralax Errors"), and the effects of varying fraction of unresolved binary stars in 
the data (last three row^s) . 



normalization / is overestimated by ~ 10% (note how- 
ever that this is still within the statistical uncertainty). 
The scale lengths are also slightly underestimated, with 
the effect less pronounced for the thick disk. 

We conclude that the Malmquist bias due to photo- 
metric errors and the dispersion around the photometric 
parallax relation has a relatively small effect on the de- 
termination of Galactic model parameters, at the level of 

4.3.4. Effects of Unrecognized Multiplicity 

Unrecognized multiplicity biases the density maps 
and the determination of Galactic model parameters 
by systematically making unresolved binary stars, when 
misidentified as a single star, appear closer then they 
truly are. It's effect is most strongly dependent on the 
fraction of observed "stars", fm, that are in fact unre- 
solved multiple systems. 

We model this effect by simulating a simplified case 
where all multiple systems are binaries. Because the frac- 
tion of binary systems is poorly known, we generate three 
mock catalogs with varying fractions fm of binary sys- 
tems misidentified as single stars, and observe the effects 
of fm on the determination of model parameters. Pho- 
tometric errors and photometric parallax dispersion (as 
discussed in Section [4. 3. 3p are also mixed in. 

The results are given in the last three rows of Table [51 
The effect of unresolved binary systems is a systematic 
reduction of all spatial scales of the model. Measured 
disk scale heights are underestimated by as much as 25% 
(/„ = 1), 20% (/,„ = 0.5) and 15% (/,„ = 0.25). Mea- 
sured scale lengths are similarily biased, with the thin 
disk scale length being underestimated by 25, 20, and 
13% and the thick disk scale length by 25, 15, and 9% 
for fm = 1, 0.5, 0.25, respectively. The thick disk density 
normalization is mildly overestimated (^ 10%) but not 
as strongly as the disk scales, and still within statistical 
uncertainty. 

4.3.5. Effects of Systematic Distance Determination Error 

We next measure the effect of systematically over- or 
underestimating the distances to stars due to absolute 
calibration errors of the photometric parallax relation. 
This can already be judged by comparing the values of 
model parameters determined from fits using the bright 
and faint photometric parallax relations (Tables [3] and S]), 
but here we test it on a clean simulated sample with a 
known underlying model. 

We generate a mock catalog by using the bright photo- 
metric parallax relation (eq. [2]) to convert from distances 



to magnitudes, and mix in SDSS photometric and par- 
allax dispersion errors (Section I4.3.3p . We process this 
catalog by assuming a parallax relation 0.5 magnitudes 
brighter, and 0.5 magnitudes fainter than Equation [21 
effectively changing the distance scale by ±23%. 

Fit results are shown in Table [HI including for compar- 
ison in the middle row the parameters recovered using 
the correct Mr{r — i) relation. The effect of system- 
atic distance errors is to comparably increase or decrease 
measured geometric scales. The thin and thick disk scale 
heights increase by 33% and 34%, and the scale lengths 
by 28% and 42%, respectively, if the distances are over- 
estimated by 23%. If they are underestimated by the 
same factor, the parameters are reduced by 14% and 19% 
(thin and thick disc scale height), 18% and 10% (thin and 
thick disk scale lengths). Interestingly, both increasing 
and decreasing the distance scale results in an increase 
of measured normalization, by a factor of 10 — 25%. 

4.3.6. Test of Cylindrical Symmetry 

In Sect ion [5TT] we argued based on the shapes of isoden- 
sity contours in Figures I12H14I and in particular in Fig- 
ure[Tl]that once the large overdensities are taken out, the 
Galactic density distribution is cylindrically symmetric. 
Therefore it was justifiable for the purpose of determin- 
ing the overall Galactic stellar number density distribu- 
tion, to measure the density along the same galactocen- 
tric annuli R and only consider and model the distribu- 
tion of stars in two-dimensional R — Z plane. 

Using Figure [25l we quantitatively verify this assump- 
tion. In the panels of the top row as solid black his- 
tograms we plot the distribution of 

Ap ^ p{R,<p,Z)~piR,Z) 

a ap{R,c^,Z) ^ > 

for four r — i color bins^^. This is the difference 
of the density measured in a pixel at (i?, (/>, Z) and 
the mean density 'p{R, Z) at annulus {R,Z), normal- 
ized by the expected Poisson fluctuation ap{R,Z) — 

^N{R,4>,Z)/V{R,4>,Z). 

The dotted red histogram in the panels shows a Poisson 
model of noise-normalized deviations expected to occur 
due to shot noise only. If all pixels were well sampled 
[N > 50 stars, which is not the case here), this distribu- 
tion would bea/i = 0,(T = l Gaussian. 

The data and the Poisson model show a high degree 
of agreement. However, in a strict statistical sense, for 

■^^ Analogous histograms of other r — i bins share the same fea- 
tures. 
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TABLE 9 

Effects of Mr{r — i) Calibration Errors 



Simulation 




Li 


Hi 


/ 


L2 


H2 


Mrir -i) - 0.5 


1.18 


3080 ± 55 


305 ±4 


0.123 ± 0.011 


4881 ± 78 


976 ± 19 


Mr {r - i) 


0.95 


2403 ± 40 


230 ±2 


0.111 ± 0.010 


3441 ± 57 


725 ± 13 


Mrir + 0.5 


1.23 


1981 ± 32 


198 ± 2 


0.138 ± 0.013 


3091 ± 52 


586 ± 11 



Note. — Effects of systematic error in the calibration of photometric paralax relation. Middle row lists the parameters 
recovered assuming the correct paralax relation (cq. [2j, from a Monte Carlo generated catalog with realistic photometric 
errors and dispersion ctmj. — O.Smag around the mean of Mj.{r — i). This is the same catalog as in row 3 of Table[8] The 
first and last row show parameters recovered when a paralax relation which systematically under-/overestimates the absolute 
magnitudes by 0.5 magnitudes (over-/underestimates the distances by ~ 23%) is assumed. 



0.10 < r-i < 0.15 0.70 < r-i < 0.80 1.00 < r-i < 1.10 1.00 < r-i < 1.10 




Fig. 25. — Distribution of noise-normalized deviations of density in pixels (X, Y, Z) from the mean density measured along their 
corresponding annuli {R = \/X'-^ + Y^, Z). Black solid histogram shows the data. Red dotted histogram shows a Poisson noise model. 
Histograms in the top row and bottom rows have been calculated assuming the bright (Equation [2J and faint (Equation [1} photometric 
paralax relation, respectively. The rightmost panel in the top row shows the same distributions derived from a Monte Carlo simulated 
catalog, with 25% unresolved binary fraction, crMr = 0-3 paralax dispersion and SDSS photometric errors. 



all but the 1.3 < r — i < 1.4 bin the data and the model 
are inconsistent with being drawn from the same distri- 
bution at a 5% confidence level. This is not very surpris- 
ing, as the effects of unresolved multiplicity and other 
observational errors may modify the residual distribu- 
tion. We verify this by examining the same statistic cal- 
culated from a Monte Carlo generated catalog with 25% 
unresolved binary fraction and mixed in photometric er- 
ror (Table [H fourth row) . The resulting "observed" and 
Poisson model histograms for the simulated catalog are 
shown in the top right-most panel of Figure [51 They 
show the same behavior as seen in the data. 

One may ask if these distributions could be used to 
further test (and/or constrain) the photometric parallax 
relations, as under- or overestimating the distances will 
break the cylindrical symmetry of Galaxy and distort the 
isodensity contours. The answer is, unfortunately, no. In 
the bottom row of Figure [55] we show the distributions 
analogous to those in the top row, but calculated from 
maps obtained using the faint parallax relation (Equa- 
tion [T]) . They are very similar to those in the top row, 
with no deviations that we can conclusively attribute to 
the change in photometric parallax relation, although 
there is an intriguing slightly better Data-Model agree- 
ment in 0.7 < r — i < 0.8 color bin for the faint, than 
for the bright relation. This is initially surprising, as one 



would intuitively expect the erroneous distance estimate 
to map the density from different real Galactocentric an- 
nuli {R, Z) to the same observed {Ro, Zo), and therefore 
widen the residual distribution. However, this effect (as 
we verified using the Monte Carlo generated catalogs) is 
indiscernible for the combination of density distribution 
seen in the Milky Way, and the portion in {R, Z) space 
where we have clean data. The region near R = where 
the assumption of cylindrical symmetry is most sensitive 
to errors in distance determination is contaminated by 
debris from the Virgo overdensity, making it unusable 
for this particular test. 

4.3.7. Resolved Substructure in the Dtsk 

The panels in Fig. [261 illustrate the fitting results and 
the revealed clumpy substructure. The columns, from 
left to right, show the data, the model and the model- 
normalized residuals. The bottom three rows are results 
of fitting a disk-only model, while the top row also in- 
cludes a fit for the halo. 

While the best-fit models are in good agreement with 
a large fraction of the data, the residual maps show 
some localized features. The most prominent feature 
is found at practically the same position (R ^6.5 kpc 
and Z ^ 1.5 kpc) in all color bins, and in figure [^Sj 
is the most prominent in the top right panel The fea- 
ture itself is not symmetric with respect to the Galactic 
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□ nta/MoclOl - 1 




Fig. 26. — Examples of model fits for four color bins, one per eacli row. Note the different scales. The left panel of each row shows the 
data, the middle panel the best-fit model and the right panel shows (data-model) residuals, normalized to the model. The residuals are 
shown on a linear stretch, from -40% to +40%. Note the excellent agreement of the data and the model for reddest color bins (bottom 
row), and an increasing number of overdensities as we move towards bluer bins. In the residuals map for the 0.35 < r — i < 0.40 bin (top 
row) the edges of the Virgo overdensity (top right) and the Monoceros stream (left), the overdensity at (i? ~ 6.5, Z ~ 1.5) kpc and a 
small overdensity at (_R ~ 9.5, Z ^ 0.8) kpc (a few red pixels) are easily discernible. The apparently large red overdensity in the south at 
(R ~ 12, Z ~ —7) kpc is an instrumental effect and not a real feature. 
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7 = +01500 pc {0.700 < r-I < O.ROO) 7 = +00600 pc (0 700 < r-i < O.ROO) 




4000 6000 8000 10000 12000 '4000 6000 8000 10000 12000 '4000 6000 8000 10000 12000 

R [pc] X [pc] X [pc] 



Fig. 27. — Ring-like deviations from the underlying density distribution detected after the best fit model subtraction. Left panel shows 
the data-model residuals in the R-Z plane, normalized to the model, for the 0.7 < r — i < 0.8 bin. Two overdensities detected on Figure [26l 
are clearly visible and marked by a dashed circle and rectangle. In the X — Y plane, shown in the middle panel, the R ~ 6.5 kpc feature 
reveals itself as a ring-like ^ 20% density enhancement over the smooth background at Z ~ 1.5 kpc. Similarly on the right panel, the 
R ~ 9.5 feature is detectable as a strong ~ 50% enhancement in the Z = 600 pc slice. 
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plane, though a weaker counterpart seems to exist at 
Z < 0. It may be connected to t he feature ob s erved by 
iLarsen fc Humphreys! (|1996D and lParker et al.l (|2003D in 
the POSS I survey at 20° < Z < 45°, 6 ~ 30° and which 
they interpreted at the time as a signature of thick disk 
asymmetry. We also show it in an X — y slice on the 
center panel of Figure d?! where it is revealed to have 
a ring-like structure, much like the Monoceros stream 
in Figure 1131 Another smaller overdensity is noticeable 
in all but the reddest Data-Model panel of figure US] at 
R ^ 9.5 kpc and Z ~ 0.8 kpc, apparently extending 
for ~1 kpc in the radial direction. When viewed in an 
X — Y slice, it is also consistent with a ring (Figure [23 
right panel); however, due to the smaller area covered in 
X — Y plane, an option of it being a localized clumpy 
overdensity is not possible to exclude. 

If this substructure is not removed from the disk as 
we have done in Section 14. 1[ it becomes a major source 
of bias in determination of model parameters. The ef- 
fect depends on the exact location, size and magnitude 
of each overdensity, and whether the overdensity is in- 
side the survey's flux limit for a particular color bin. For 
example, the effect of {R ^6.5, Z ^ 1.5) overdensity was 
to increase the scale of the thick disk, while reducing 
the normalization, to compensate for the excess number 
density at higher Z values. The effect of R ^ 9.5 over- 
density was similar, with an additional increase in scale 
length of both disks to compensate for larger than ex- 
pected density at i? > 9 kpc. Furthermore, these effects 
occur only in bins where the overdensities are visible, 
leading to inconsistent and varying best-fit values across 
bins. Removal of the clumps from the dataset prior to 
fitting the models (Section |4.1|) restored the consistency. 

4.3.8. Statistics of Best-fit Residuals 

The best-fit residuals of data where no apparent sub- 
structure was detected may hold statistical evidence of 
unrecognized low-contrast overdensities and dumpiness. 
If there are no such features in the data, the distribu- 
tion of residuals will be consistent with a Poisson noise 
model. Conversely, if a substantial degree of unresolved 
dumpiness exists, the distribution of residuals will be 
wider and may be skewed compared to the distribution 
expected from Poisson noise only. 

We begin by inspecting the statistics of residuals in 
R, Z plane, shown in Figure [551 for four color bins rep- 
resentative of the general behavior. The solid black his- 
togram shows the distribution of Data-Model deviation 
normalized by the shot noise. As pixels in (i?, Z) plane 
are well sampled (typically, N^tars > 50 in more than 
95% of pixels), the shot noise induced errors are Gaus- 
sian, and the residuals are expected to be normally dis- 
tributed {N{fj, = 0, cr = 1), dotted red curve) for a perfect 
model fit. The distribution seen in the data is wider than 
the expected Gaussian, by 10% at the red end and 85% 
at the blue end. 

Two-dimensional (i?, Z) plane pixels contain stars from 
all (potentially distant) X ~Y positions in the observed 
volume which map to the same R, thus making the 
the residuals difficult to interpret. We therefore con- 
struct analogous distributions of residuals for pixels in 
3D {X, Y, Z) space. They are shown in the panels of left 
column of Figure (solid black histograms). As this 
time not all (A", F, Z) pixels are well sampled, a full Pois- 



sonian noise model is necessary to accurately model the 
distribution of residuals. We overplot it on panels of Fig- 
ure [55| as a dotted red histograms. In the left column of 
the same figure, we also plot the measured distribution 
of model normalized residuals (solid black histogram), 
and the Poisson model prediction for residuals due to 
shot-noise only (dotted red histogram). To judge the ef- 
fects of observational errors and unresolved multiplicity, 
the bottom two panels show the distributions measured 
from a Monte Carlo generated catalog with 25% unre- 
solved binary fraction and photometric error (Table [H 
fourth row). Comparison of data and Poisson models, 
and the observed and simulated distributions, leads us 
to conclude that across all examined color bins, the dis- 
tribution of deviations is consistent with being caused by 
shot noise only. 

This is in apparent conflict with the analysis of residu- 
als in 2D (i?, Z) plane. The key in reconciling the two is 
to notice that different spatial scales are sampled in 3D 
and 2D case. The 3D analysis samples scales compara- 
ble to the pixel size. The effective sampling scale is made 
variable and somewhat larger by the smearing in line-of- 
sight direction due to unrecognized stellar multiplicity, 
but is still on order of not more than a few pixel sizes. 
On the other hand, the effective scale in 2D is the length 
of the arc over which 3D pixels were averaged to obtain 
the R—Z maps. This is on order of few tens of percent of 
the faint volume-limiting distance (Table [H column Di) 
for each bin. The deviations seen in 2D maps are there- 
fore indicative of data-model mismatch on large scales, 
such as those due to large scale overdensities or simply 
due to the mismatch of the overall shape of the analytic 
model and the observed density distribution. 

In support of this explanation in Figure [30] we plot a 
rainbow-coded shot noise normalized map of residuals in 
pixels at Z = 10 kpc slice, 0.1 < r — i < 0.15 color bin. 
On large scales a small but noticeable radial trend in the 
residuals is visible, going from slightly underestimating 
the data (larger median of residuals, more red pixels) at 
smaller R towards overestimating the data near the edge 
of the volume at higher R (smaller median of residuals, 
more blue pixels). This trend manifests itself as widening 
of residual distribution (and increase in x^) in Figure (281 

The small scale fluctuations are visible as the "nois- 
iness" of the data. They are locally unaffected by the 
large-scale trend, and consistent with just Poisson noise 
superimposed on the local density background. If exam- 
ined in bins along the R direction, the large scale trend 
does leave a trace: the median of residuals is slightly 
higher than expected from the Poisson model at low 
R and lower than expected at high R. But when the 
residuals of all pixels are examined together, this signal 
disappears as the opposite shifts from lower and higher 
radii compensate for each other. This leaves the residual 
distribution in Figure [29| consistent with being entirely 
caused by shot-noise. 

We conclude from this admittedly crude but neverthe- 
less informative analysis that i) it rules out significant 
dumpiness on scales comparable to the pixel size of each 
color bin ii) demonstrates there are deviations on scales 
comparable to radial averaging size, indicating the func- 
tional forms of the model do not perfectly capture the 
large-scale distribution, and iii) shows that these devia- 
tions are negligible for the disk and pronounced for the 
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Fig. 28. — Distribution of residuals in {R, Z) plane pixels. The solid black histogram shows the data, and the overplotted solid black 
curve is a Gaussian distribution with dispersion a determined from the interquartile range of the data. For comparison, the dotted red line 
shows a CT = 1 Gaussian, the expected distribution if the residuals were due to shot noise only. 



halo, pointing towards a need for halo profiles more com- 
plicated than a single power law. 

4.3.9. Wide Survey Area and Model Fit Degeneracies 

In a model with as many as ten free parameters, it is 
not easy to assess the uniqueness of a best-fit solution, 
nor to fully understand interplay between the fitted pa- 
rameters. We show two illuminating examples of fitting 
degeneracies. 

In Fig. [31] we plot the density distributions for two sig- 
nificantly different models: a thin plus thick disk model 
without a halo, and a single disk plus halo model. De- 
spite this fundamental intrinsic difference, it is possi- 
ble to fine-tune the model parameters to produce nearly 
identical Z dependence of the density profiles at i? = 
8 kpc. As shown in the bottom panel, significant differ- 
ences between these two models are only discernible at 
|Z| > 3 kpc and R significantly different from 8 kpc. 

Secondly, in the left column of Fig. [32] we reproduce 
the top two panels of Fig. [15] The density profile is 



well described by two exponential disks of scale heights 
Hi = 260 and H2 = 1000 and normalization of 4%. In 
the right column of the figure we plot the same data, 
but overplotted with best fit models from Table [H] The 
scales in this model are Hi = 245 and H2 — 750, with 
thick-to-thin normalization of 13%, and the bottom right 
panel also includes a contribution of the halo. Although 
significantly different, the two models are here virtually 
indistinguishable . 

This is a general problem of pencil beam surveys with 
a limited sky coverage. A single pencil beam and even a 
few pencil beams (depending on the quality of the data 
and positioning of the beams) cannot break such model 
degeneracies. We speculate that this in fact is likely the 
origin of some of the dispersi on in disk parame ter values 
found in the lite rature (e.g.. ISiegel et al.ll2002l . Table 1; 
iBilir et al.ll200i ). 

In our case, while we have not done a systematic search 
for degenerate models leading to similar given our sur- 
vey area, we have explored the possibility by attempting 
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Fig. 29. — Left column: the Poisson-noise normalized distri- 
bution of residuals in three-dimensional {X, Y, Z) pixels for three 
representative color bins. Right column: model-normalized distri- 
bution of residuals in each pixel, (Data - Model) / Model. The 
solid black histograms show the data, while the dotted red his- 
tograms show the expectation from residuals due to Poisson noise 
only. The bottom row shows the same distributions derived from a 
Monte Carlo simulated catalog, with 25% unresolved binary frac- 
tion, cTMr = 0-3 paralax dispersion and SDSS photometric errors. 

a 100 refits of the data starting witfi random initial pa- 
rameter values. In case of fits to individual bins, we 
find local minima, higher by ~ 20 — 30% than the 
global minimum, with parameter values noticeably dif- 
ferent from the best fit solution. However, when jointly 
fitting all r — z > 1.0 color bins, in all cases the fit either 
fails to converge, converges to a local minimum that is 
a factor of few higher than the true minimum (and pro- 
duces obviously spurious features in maps of residuals), 
or converges to the same best-fit values given in Tables [S) 
SDSS therefore seems largely successful in breaking the 
degeneracies caused by the limited survey area and pho- 
tometric precision, leaving local departures from expo- 
nential profiles as the main remaining source of uncer- 
tainty in best-fit model parameters. 

4.3.10. Physical Basis for the Density profile 
Decomposition into Disks and the Halo 

Although the density profiles shown in bottom right 
panel of Fig. [32] and the bottom panel of Fig. [15] ap- 




-10000 10000 20000 30000 

X [pc] 

Fig. 30.— A Z = 10 kpc X -Y slice of data from the 0.10 < 
r — i < 0.15 color bin. Only pixels with less than 5% of the density 
coming from the Galactic disk component are shown. The colors 
encode noise-normalized residuals in each pixel, saturating at —3cr 
(purple) and +3cr (red). The large red spot in the top part of the 
figure is due to the Virgo overdensity (this region was excluded 
during model fitting; it is shown here for completeness only). 

pear with high signal-to-noise ratios, it may be some- 
what troubling that as our range of observed distances 
expands, we need to keep introducing additional compo- 
nents to explain the data. Are these components truly 
physically distinct systems, or largely phenomenological 
descriptions with little physical basis? 

The question is impossible to answer from number den- 
sity data alone, and two companions papers use metallic- 
ity estimates (Paper II) and kinematic information (Pa- 
per III) to address it. Here we only look at a subset of this 
question, namely the differentiation between the disk and 
halo components. Disk stars (Population I and interme- 
diate Population II) have metallicities on average higher 
by about 1-2 dex than that of the halo. Such a large 
difference in metallicity affe cts the u — g color of turn-off 
stars (e.g.. lChen et al.ll2001h . An analysis of SDSS colors 
for Kurucz model atmospheres suggests that stars at the 
tip of the stellar locus with 0.7 < u — g < 1 necessarily 
have metallicities lower than about —1.0. These stars 
also have markedly different kinematics further support- 
ing the claim that they are halo stars (Paper II and HI) . 

We select two subsamples of stars from the 0.10 < 
r — i < 0.15 color bin: low metallicity halo stars with 
0.60 < u — g < 0.95, and high metallicity disk stars with 
0.95 < u — g < 1.15. This separation is of course only ap- 
proximate and significant mixing is expected both at the 
faint end (disk stars contaminated by the more numer- 
ous halo stars) and at the bright end (halo stars contam- 
inated by the more numerous disk stars). Nevertheless, 
the density profiles for these two subsamples, shown in 
Fig. [33] are clearly different. In particular, the disk pro- 
file is much steeper, and dominates for Z < 3 kpc, while 
the halo profile takes over at larger distances from the 
Galactic plane. This behavior suggests that the multiple 
components visible in the bottom panel in Fig. [15] are 
not an over-interpretation of the data. 

In addition to supporting a separate low-metallicity 
halo component, this test shows that a single exponential 
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Fig. 31. — An illustration of the degeneracies in fitting models for 
stellar distribution. The top panel shows a thin disk plus thick disk 
model, without any contribution from the halo (volume density on 
a logarithmic stretch, from blue to red, shown only for the regions 
with SDSS data), and the middle panel shows a single disk plus an 
oblate halo model. Both models are fine-tuned to produce nearly 
identical counts for R = 8 kpc and \Z\ < 8 kpc. The bottom panel 
shows the difference between the two models (logarithmic stretch 
for it a factor of 3, from blue to red, the zero level corresponds to 
green color). The models are distinguishable only at \Z\ > 3 kpc 
and R significantly different from 8 kpc. 



disk model is insufficient to explain the density profile 
of high-metallicity stars. This is ordinarily remedied by 
introducing the thick disk. However, with only the data 
presented here, we cannot deduce if the division into thin 
and thick disk has a physical basis or is a consequence 
of our insistence on exponential functions to describe the 
density profile. 

4.3.11. The Corrected Best Fit Parameters 



In Section |4~2| we have used two samples of stars to fit 
the parameters of the disk: the 1.0 < r — i < 1.4 sample 
of M dwarfs, and 0.65 < r — i < 1.0 sample of late K / 
early M dwarfs. Best fit results obtained from the two 
samples are very similar, and consistent with the stars 
being distributed in two exponential disks with constant 
scales across the spectral types under consideration. 

The fit to 0.65 < r — i < 1.0 sample required an ad- 
dition of a third component, the Galactic halo. This, 
combined with the photometric parallax relations that 
are inappropriate for low metallicity stars in this color 
range, may bias the determination of thick disk parame- 
ters. For example, while the measured scale height of the 
thick disk in 0.65 < r — i < 1.0 range is ^ 10% lower than 
in 1.0 < r — i < 1.4 range, it is difficult to say whether 
this is a real effect, or interplay of the disk and the halo. 

Furthermore, we detected two localized overdensities 
in the thick disk region (Section I4.3.7p . While every 
effort was made to remove them from the data before 
fitting the model, any residual overdensity that was not 
removed may still affect the fits. If this is the case, the 
0.65 < r — i < 1.0 bins are likely to be more affected 
than their redder counterparts, being that they cover a 
larger volume of space (including the regions where the 
overdensities were found). 

For these reasons, we prefer the values of disk param- 
eters as determined from 1.0 < r — i < 1.4 sample, as 
these are a) unaffected by the halo and b) least affected 
by local overdensities. 

Other dominant sources of errors are (in order of de- 
creasing importance) i) uncertainties in absolute calibra- 
tion of the photometric parallax relation, ii) the misiden- 
tification of unresolved multiple systems as single stars, 
and iii) Malmquist bias introduced by the finite width of 
Mr{r — i) relation. Given the currently limited knowl- 
edge of the true photometric parallax relation (Figure [2]), 
there is little one can do but try to pick the best one 
consistent with the existing data, and understand how 
its uncertainties limit the accuracy of derived parame- 
ters. Out of the two relations we use (bright, eq. [H and 
faint, eq. [T]), we prefer the bright normalization as it is 
consistent with the kinematic data (Paper III) and the 
analysis done with wide binary candidates (Section l2.2.2p 
shows its shape to be correct to better than O.lmag for 
r — i > 0.5. If we are mistaken, as discussed in Sec- 
tion l4.3.5l errors in Af,. of AMr = ±0.5 will lead to errors 
of 20 — 30% in parameter estimation. Given Figure[2]and 
the analysis of wide binary candidates in Section r2.2.2l we 
believe this to be the worst case scenario, and estimate 
that the error of each scale parameter is unlikely to be 
larger than ±20%. 

The dependence of best-fit parameters derived from 
mock catalogs on multiplicity (binarity) is shown in Fig- 
ure 1341 The challenge in correcting for multiplicity is 
knowing the exact fraction of observed "stars" which are 
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Fig. 32. — An illu stration of degeneracies present in fitting of Galactic models. The two panels in the left column are the same top 
two panels of Fig. 1151 The panels to the right show the same data, but are overplotted with best fit models from Table |6] In spite of 
substantially different best fit values, the two models are virtually indistinguishable when fitting the R = 8kpc, ±Z direction of the data. 



TABLE 10 
The Galactic Model 



Parameter 


Measured 


Bias-corrected Value 


Error estimate 


Zo 


25 




20% 


Li 


2150 


2600 


20% 


Hi 


245 


300 


20% 


/ 


0.13 


0.12 


10% 


L2 


3261 


3600 


20% 


H2 


743 


900 


20% 


fh 


0.0051 




25% 


1 


0.64 




<o.i 


n 


2.77 




< 0.2 



Note. — Best-fit Galactic model parameters (sec CQS. I21H24t , as directly measured from the apparent 
number density distribution maps (2"^ column) and after correcting for a 35% assumed binary fraction 
and Malmquist bias due to photometric errors and dispersion around the mean of the photometric 
paralax relation (3^^ column). 



unresolved multiple systems. While it is understood that 
a substantial fraction of Galactic field stars are in bi- 
nary or multiple systems, its exact value, dependence 
on spectral type, population, and other factors is still 
poorly known. Measure ments range from 20% for late 
type fL. M. K dwarfs - iReid et al.l l200d iReid fc Gizii 
119971: iFischer fc Marcvlll992D to upward of 6 0% for early 
tvpes fG dwarfs: 'Duqu ennov fc MavoH[l99ll ). The com- 
position (mass ratios) of binaries also poorly constrained, 
but appears to show a preference tow ards more equal- 
mass companions in late spectral types (jReid et al.ll2006l 
Fig. 8). Given our least biased disk parameters were de- 
rived from from the M-dwarf sample {r — i > 1.0), we 
choose to follow Reid & Gizis ( 1997) and adopt a binary 
fraction of 35%. 
We accordingly revise Li, Hi and H2 upwards by 15% 



and L2 by 10% to correct for multiplicity (see Fig- 
ure [34|l . We further include additional 5% correction 
due to Malmquist bias (Section I4.3.3p . and for the same 
reason correct the density normalization by -10%. The 
final values of measured and corrected parameters are 
hsted in Table [lOl 

5. THE VIRGO OVERDENSITY 

The X — Y projections of the number density maps at 
the heights above 6 kpc from the Galactic plane show a 
strong deviation from expected cylindrical symmetry. In 
this Section we explore this remarkable feature in more 
detail. We refer to this feature as "the Virgo overden- 
sity" because the highest detected overdensity is in the 
direction of constellation Virgo, but note that the feature 
is detectable over a thousand square degrees of sky. 
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Fig. 33. — The vertical (Z) distribution of SDSS stellar counts 
for R = 8 kpc, and 0.10 < r — i < 0.15 color bin. Stars are 
separated by their u — g color, which is a proxy for metallicity, into 
a sample representative of the halo stars (low metallicity, 0.60 < 
u — g < 0.95, circles) and a sample representative of the disk stars 
(high metallicity, 0.95 < u — g < 1.15, triangles). The line in 
the top panel shows the sum of the counts for both subsamples. 
The counts for each subsample are shown separately in the middle 
and bottom panels, and compared to the best fit models, shown 
as lines. Note that the disk stars are more concentrated towards 
the Galactic plane. Due to a simple u — g cut, both samples are 
expected to suffer from contamination: close to the Galactic plane 
(|Z| < 2 kpc) the halo sample is contaminated by the disk stars, 
while further away from the plane {\Z\ > 5 kpc) the disk sample is 
contaminated by halo stars. 

5.1. The Extent and Profile of the Virgo overdensity 

To quantily the extent and profile of the Virgo over- 
density, we consider the data in an X' — Z plane, per- 
pendicular to the Galactic symmetry axis, and is rotated 
from the X — Z plane by (/> = 30° clockwise around the 
Z axis. In Figures [T^ this plane would be seen edge on, 
as a straight line at a 30° angle from the X axis, passing 
through the Galactic center and the approximate center 
of the Virgo overdensity. Note that in this plane the dis- 



tance measured along the X' axis is just the cylindrical 
galactocentric radius R. 

In the top left panel of Figure [35] we show the corre- 
sponding number density map for the bluest color bin. 
Isodensity contours show a significant deviation from the 
expected monotonic decrease with X'{= R). Instead, 
they reveal the existence of an overdense region around 
X' ~ 7-8 kpc and Z ^ 10 kpc. This overdensity is also 
visible in the density profiles at Z = 10 kpc above the 
plane, shown in Y' > —3 kpc panels of Fig. [M] As dis- 
cernible from these figures, the Virgo overdensity is re- 
sponsible for at least a factor of 2 number density excess 
at Z = 10 kpc. 

To analyze this feature in more detail, we subtract a 
best-fit Galactic model from the data shown in the top 
right panel of Figure [35l We first fit a model described 
by Equations I^TI - f^ to the observations having Y < (or 
equivalents, 180° < / < 360°). As evidenced by Fig.[T2l 
this region does not seem significantly affected by the 
overdensity. We show the difference of the data from top 
right panel of Figure [551 and the so obtained model in the 
top middle panel of the same figure. The top right panel 
shows the same difference but normalized to the model. 

The model-normalized map reveals much more clearly 
the extent and location of the overdensity. A significant 
density excess (up to a factor of 2) exists over the entire 
sampled range of Z (6 < Z/kpc < 20). Importance 
of the overdensity, relative to the smooth Milky Way 
halo background, increases as we move away from the 
Galactic plane. This increase is however mainly due to 
a fast power-law decrease of the number density of the 
halo, which causes the increase in Virgo-to-MW ratio. 
The number density of stars belonging to the overdensity 
actually increases towards the Galactic plane, as seen in 
the top middle panel. 

For comparison, analogous density and residual plots 
from a parallel plane at y = — 9 kpc is shown in the 
bottom row of Figure [35] These show no large scale 
deviations from the model. The density contours rise 
smoothly and peak near X' — 0, the point closest to the 
Galactic center. The same is seen in Y' < —5 kpc slices 
of Figure [36l 

Because no local maximum of the overdensity is de- 
tected as Z approaches the observation boundary at 
Z — 6 kpc, with the data currently available we are un- 
able to quantify its true vertical (Z) extent. It is pos- 
sible that it extends all the way into the Galactic plane 
and, if it is a merging galaxy or a stream, perhaps even 
to the southern Galactic hemisphere. In the direction 
of Galactic radius, the Virgo overdensity is detected in 
the 2.5 < X'/kpc < 12.5 region. The X' position^^ of 
maximum density appears to shifts slightly from X' ~6 
kpc at Z — 6 kpc to X' ~7 kpc at Z = 15 kpc. The 
width ("full- width at half density") decreases by a fac- 
tor of ~ 2 as Z increases from 6 to 20 kpc. While not 
a definitive proof, these properties are consistent with a 
merging galaxy or stream. 

The thickness of the overdensity in the direction per- 
pendicular to the plane of the image in Figure [351 f the Y' 
direction) is no less than ~ 5 kpc. As in the case of the 
Z direction, the true extent remains unknown because of 



33 Note that in this Y' 
cylindrical radius. 



plane X' = R, the galactocentric 
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Fig. 34. — Effect of unrecognized binarity on fits of model parameters, derived from simulations listed in Table |8] Each of the panels 
shows the change of a particular model parameter when a fraction /(, of observed "stars" are unrecognized binary systems. 



the current data availability. Note that the size of the 
overdensity seen in the maps in the direction of the line 
of sight towards the Sun is a combination of true size 
and the smearing induced by the photometric measure- 
ment and parallax errors (Fig. [7|) and (most significantly) 
the effects of unrecognized stellar multiplicity. The true 
line of sight extent is therefore likely smaller, by at least 
30-35%. 

5.2. Direct Verification of the Virgo Overdensity 

Significant data processing was required to produce 
maps such as the one revealing the Virgo overdensity 
(e.g. the top panels in Fig. [55]) . In order to test its ex- 
istence in a more direct, and presumably more robust, 
way, we examine the Hess diagrams constructed for the 
region of the sky that includes the maximum overden- 
sity, and for a control region that appears unaffected by 
the Virgo feature. The boundaries of these two regions, 
which are symmetric with respect to the / = line, the 
corresponding Hess diagrams, and their difference, are 
shown in Fig. 1371 



The top left panel of Fig. [37] shows the northern (in 
Galactic coordinates) sky density of stars with 0.2 < g — 
r < 0.3 and 20 < r < 21 in the Lambert equal area 
projection of Galactic coordinates (the north Galactic 
pole is in the center, l—O is towards the left, and the 
outermost circle is 5 = 0°). This map proj ection does 
not p reserve shapes (it is not conformal, e.g. iGott et al] 
l2005f ). but it preserves areas - the area of each pixel in 
the map is proportional to solid angle on the sky, which 
makes it particularly suitable for study and comparison 
of counts and densities on the celestial sphere. The color 
and magnitude constraints select stars in a D ~ 18 kpc 
heliocentric shell, and can be easily reproduced using the 
publicly available SDSS database. The Virgo overdensity 
is clearly visible even with these most basic color and 
magnitude cuts, and extends over a large patch of the sky, 
roughly in the I = 300°, 6 = 65° direction. The overall 
number density distribution is clearly not symmetric with 
respect to the horizontal I = 0, 180 line. For example, in 
a thin shell at r ~ 21mag there are 1.85 ± 0.03 times 
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Fig. 35. — The top left panel shows the distribution of stellar number density similar to that in Fig. 1101 except that here we only show the 
data from a narrow Y' = slice in a X' , Y' , Z' coordinate system defined by rotating the X, Y, Z galactocentric system counterclockwise 
by (/) = 30° around the Z axis. In these coordinates, the Y' = plane cuts vertically through the center of the Virgo overdensity. The top 
middle panel shows the difference of the observed density and a best-fit model constrained using the data from the Y < region. The right 
panel show the same difference but normalized to the model. The bottom panels display analogous slices taken at Y' = —9 kpc. Compared 
to the top row, they show a lack of any discernable substructure. 



more stars in the I = 300°, b = 65° direction, than in the 
corresponding symmetric {I — 60°, b — 65°) direction, a 
~ 28(7 deviation from a cylindrically symmetric density 
profile. When the color range is sufficiently red (e.g. 
0.9 < g — r < 1.0), and in the same magnitude range, the 
asymmetry disappears (not shown). These stars have a 
smaller absolute magnitude, are therefore much closer, 
and do not go far enough to detect the overdensity. 

The two right panels in Fig. [571 show the Hess diagrams 
for two 540 deg^ large regions towards {I — 300°, b — 60°) 
and {I = 60°, b = 60°), and the bottom left panel shows 
their difference. The difference map reveals a strong over- 
density at g — r ~ 0.3 and r > 20. A more quantita- 
tive analysis of the Hess diagram difference is shown in 
Fig. 1381 For red stars the counts in two regions are indis- 
tinguishable, while for blue stars the counts difference is 
highly statistically significant. There is no indication for 
a turnover in blue star number count difference, which 
strongly suggests that the Virgo overdensity extends be- 
yond the SDSS faint limit. We conclude that the Hess 
diagram analysis robustly proves the existence of a sig- 
nificant star count overdensity towards I = 300°, b = 65°, 
from approximately r ~ 18 mag to r 21.5 mag. 

From the diagram in bottom left panel of Fig. [371 
a crude estimate of the surface brightness of the over- 
density can be made by summing up the fluxes of all 
stars in the CMD and dividing the total flux by the 
area observed. To isolate the overdensity, we only count 
the fluxes of stars satisfying 0.2 < g — r < 0.8 and 
18 < r < 21.5. This will effectively be a lower limit, 
because we will miss stars dimmer than the limiting mag- 



nitude (r — 21.5), and bright giants (r < 18). We obtain 
a value of: 

— 32.5 mag arcsec"^ (28) 

This is about a magnitude and a half fainter than the 
surface brightness of Sagi ttarius dwarf northern stream 
(Ey ^ 31 magarcsec"^; iMartmez-Delgado et al.ll20011 . 
iMartfnez-Delgado et al.l 120041 ). 

Assuming the entire overdensity covers ~ 1000 deg^ 
of the sky (consistent with what is seen in the top left 
panel of Fig. l37)) . and is centered at a distance of D ~ 10 
kpc, from the surface brightness we obtain an estimate 
of the integrated absolute r band magnitude, Mr = 
—7.7 mag. This corresponds to a total luminosity of 
Lr = O.O9-1O^L0, where we calculate the absolute r band 
magni tude of the Sun to be MrQ = 4.6, using eqs. 2 and 
3 from lJuric eFaLl ([2002D. and adopting (B-V )p, = 0.65 
and Vq = 4.83 from'Bi nnev fc Merrifieldl (|1998f) . This lu- 
minosity estimate is likely uncertain by at least a factor 
of a few. Most of the uncertainty comes from the un- 
known exact distance and area covered by the overden- 
sity. Uncertainty due to the flux limit depends on the 
exact shape of the luminosity function of stars making 
up the overdensity, but is likely less severe. For example, 
assuming that the luminosity function of the o verdensity 
is sim ilar to that of the Solar neighborhood (Reid et al.l 
[20021 Table 4), and that our sample of overdensity stars 
is incomplete for g — r > 0.5 (see bottom left panel of 
Fig. [57|) . the corrected luminosity and absolute magni- 
tude are Lr = 0.10 • lO^i© and Mr = —7.8 (note that 
only the red end of the luminosity function is relevant 
here). Taking a more conservative incompleteness bound 
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Fig. 36. — The distribution of stellar number density for 0.10 < r — i < 0.15 color bin at Z = 10 kpc a bove the Galactic plane. Each 
panel shows a narrow Y' crossection in coordinate system defined by = 30° (see the caption of Figure l35t . Note a clear and strong 
density excess around X' ~ 8 kpc in V > panels, coming from the Virgo overdensity. 



oi g — r > 0.3, the luminosity grows to Lr = 0.11 • IO^Lq 
(22% difference), or Mr — —8, in terms of absolute mag- 
nitude. Again, these are all lower limits. 

5.3. Metallicity of the Virgo Overdensity 

The SDSS u band measurements can be used to gauge 
metallicity of the Virgo overdensity. As already discussed 
in Section [3.3.1l stars at the tip of the stellar locus (0.7 < 
u — g < 1) typically have metallicities lower than about 
— 1.5. This u — g separation can be improved by using 
instead t he principal axes in the g — rvs.u — g color-color 
diagram (jivezic et al.ll2004b( l 

- 1.28 (29) 
0.234 (30) 



Pu = 0.415(5 -r) + 0.910(u - g) 
P2s = 0.545(5 -r)~ 0.249(u - g) ^ 



The main sequence stars can be isolated by requiring 

- 0.06 < P2s < 0.06 (31) 
and the low-metallicity turn-off stars using 

- 0.6 < Pis < -0.3, (32) 

with Pis = —0.3 approximately corresponding to 
[Fe/H] = -1.0. 

In Fig. [39] we show Hess diagrams of Pis color vs. r 
magnitude for the Virgo overdensity field and the control 
field, and their difference. A significant excess of stars 
with Pis < —0.3 exists in the Virgo overdensity field, 



while there is no statistically significant difference in star 
counts for stars having Pis > —0.3. The observed Pis 
distribution implies metallicities lower than those of thick 
disk stars, and similar to those of the halo stars (see also 
Paper II). 

5.4. Detections of Related Clumps and Overdensities 

There are a number of stellar overdensities reported 
in the litera ture that are p robably related to the Virgo 
overdensity. Newb erg et al.l (2002) searched for halo sub- 
structure in SDSS equatorial strips (d ~ 0) and reported 
a density peak at {l,b) ~ (297,63). They tentatively 
concluded that this feature is "a stream or other diffuse 
concentration of stars in the halo " and pointed out that 
follow-up radial velocity measurements arc required to 
ascertain that the grouping is not a product of chance 
and statistics of small numbers. 

Detections of RR Lyrae stars are particularly use- 
ful because they are excellent standard can dles. Using 
RR L yrae detect ed by the QUES T survey, IVivas et"aLl 
(l200l . see also IZinn et al] feoOSf ) discovered an over- 
density at ~ 20 kpc from the Galactic center at 
{l,b) ~ (314,62) (and named it the "12^4 clump"). 
The same clump is discernible in t he SDSS candidate 
RR Lyrae sample (jiv ezic ct al. 2000; ^Ivezic et al.ll2003"at 
llvezic et aLll2003Hl^zic et al . 2003bl. M ore rec ently, the 
NSVS RR Lvrae survev (fWozniak et^l2004l ) detected 
an overdensity in the same direction, and at distances ex- 
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Fig. 37. — The top left panel shows the sky density of stars with b > 0°, 0.2 < g — r < 0.3 and 20 < r < 21 in the Lambert projection 
(concentric circles correspond to constant Galactic latitude; equal area corresponds to equal solid angle on the sky) of Galactic coordinates 
(the north Galactic pole is in the center, (=0 is towards the left, and the outermost circle is 6 = 0°). The number density is encoded with 
a rainbow color map and increases from blue to red. Note that the sky density distribution is not symmetric with respect to the horizontal 
I = 0, 180 line. When the stellar color range is sufficiently red (e.g. 0.9 < g — r < 1.0), this asymmetry disappears (not shown). The 
two right panels show the Hess diagrams for two 540 deg^ large regions towards (I = 300°, fe = 60°, top) and (Z = 60°, fe = 60°, bottom), 
marked as polygons in the top left panel. The bottom left panel shows the difference of these Hess diagrams - note the strong statistically 
significant overdensity at g — r ~ 0.3 and r > 20. The pixel size in each of the three Hess diagrams is {d{g — r), dr) = (0.033, 0.1). 

tending to the sample faint limit, corresponding to about 
~6 kpc (P. Wozniak, private communication). 

2MASS survey offers an important advantage of an all- 
sky view of the Milky Wa y. We h ave followed a proce- 
dure developed by Maiewsk i et al.l ((2003) to select M gi- 
ant candidates from the public 2MASS database. We 
use M giant candidates that belong to the Sgr dwarf 
stream to fine-tune selection criteria. We also estimate 
the mean K band absolute magnitude by ty ing it to the 
stream distance implied by RR Lyrae stars (jivezic et al.l 
I2003bllah We adopt 1.0< J-K <1.3 and 9.2< K <10.2 
as the color-magnitude selection of M giant candidates 
at a mean distance of 10 kpc. 

Using a sample of 75,735 candidates selected over the 
whole sky (dominated by stars in the Galactic plane) , we 
study their spatial distribution in the high galactic lati- 



tude regions (see Fig. [40|) . We find a significant excess of 
candidate giants in the Virgo overdensity area, compared 
to a region symmetric with respect to the I = line, 
with the number ratio consistent with the properties of 
the Virgo overdensity inferred from SDSS data. For ex- 
ample, in a subsample restricted to 55° < b <80°, there 
are 66 stars with 240< I <360, and only 21 stars with 
0< I <120, with the former clustered around / ~300. 
There is no analogous counts asymmetry in the southern 
Galactic hemisphere. 

5.5. A Merger, Tri-axial Halo, Polar Ring, or? 

The Virgo overdensity is a major new feature in the 
Galactic halo: even within the limited sky coverage of the 
available SDSS data, it extends over a thousand square 
degrees of sky. Given the well defined overdensity out- 
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Fig. 38. — Quantitative analysis of t he H ess diagram difference 
shown in the bottom left panel in Fig. 1371 The left column cor- 
responds to the color bin 0.2<g — r<0.3 that is representative 
of the Virgo overdensity, and the right column is a control sam- 
ple with stars satisfying 1.2 < g — r < 1.3. The top panels show 
the counts difference as a function of apparent magnitude, and the 
middle panels shows the counts ratio. The inset in the middle 
right panel shows a histogram of the counts ratio for r < 21.5. 
The bottom panels show the counts difference normalized by the 
expected Poisson fluctuations. Note that for red stars the counts 
are indistinguishable, while for blue stars the difference is highly 
statistically significant. 

line, low surface brightness and luminosity, its most plau- 
sible interpretation is a tidally disrupted remnant of a 
merger event involving the Milky Way and a smaller, 
lower-metallicity dwarf galaxy. However, there are other 
possibilities. 

An attempt may be made to explain the detected 
asymmetry by postulating a non-axisymmetric compo- 
nent such as a triaxial halo. T his alternative is pa rticu- 
larly interesting because iNewberg fc Yannvl (|2005D , who 
essentially used the same data as analyzed here, have 
suggested that evidence for such a halo exists in SDSS 
starcounts. A different data analysis method employed 
here - the three-dimensional number density maps - sug- 
gests that the excess of stars associated with the Virgo 
overdensity is not due to a triaxial halo. The main argu- 
ment against such a halo is that, despite its triaxiality, it 
still predicts that the density decreases with the distance 
from the Galactic center. But, as shown in Figs. [T^ and 
I36i the observed density profile has a local maximum 
that is not aligned with the Galactic center. This can 
still be explained by requiring the axis of the halo not to 
be coincident with the Galactic axis of rotation. How- 
ever, even this model requires the halo density to attain 
maximal value in the Galactic center, and as seen from 
figure [551 a modest linear extrapolation of Virgo overden- 
sity to Z — still keeps it at i? ~ 6 kpc away from the 



Galactic center. Unless one is willing to resort to mod- 
els where the center of the stellar halo and the center of 
the Milky Way disk do not coincide, tri-axial halo cannot 
explain the geometry of Virgo overdensity. 

Although this makes the explanation of Virgo as a sig- 
nature of triaxial halo unlikely, it does not preclude the 
existence of such a halo. Unfortunately, it would be very 
difficult to obtain a reliable measurement of the halo tri- 
axiality with the currently available SDSS data because 
of contamination by the Virgo overdensity and uncertain- 
ties about its true extent. As more SDSS and other data 
become available in other parts of the sky, it may be- 
come possible to mask out the overdensity and attempt 
a detailed halo fit to reveal the exact details of its shape 
and structure. 

Another possible explanation of the overdensity is a 
"polar ring" around the Galaxy. This possibility seems 
much less likely than the merger scenario because there is 
no visible curvature towards the Galactic center at high 
Z in Fig.[35l Indeed, there seems to be a curvature in the 
opposite sense, where the bottom {Z ^ 6 kpc) part of the 
overdense region appears to be about 0.5 — 1 kpc closer 
to the Galactic center than its high-Z part. In addition, 
there is no excess of 2MASS M giant candidates in the 
southern sky that could be easily associated with the 
northern Virgo overdensity'^^. 

Finally, the coincidence of this overde nsity and the 
Virgo galaxy supercluster (|Binggelil 119991 ) could raise a 
question whether the overdensity could be due to faint 
galaxies that are misclassified as stars. While plausible 
in principle, this is most likely not the case because the 
star/galaxy classifier is kn own to be robust at the 5% 
level to at least r — 21.5 (jivezic et all 120021 ) ■ the over- 
density is detected over a much larger sky area (1000 deg^ 
vs. 90 deg^), and the overdensity is confirmed by ap- 
parently much brighter RR Lyrae stars and M giants. 

6. DISCUSSION 

6.1. A Paradigm Shift 

Photometric parallax methods have a long his- 
tory of use in studies of the Milky Way structure 
(e.g., iG ilmorc & Rcid 1983, K uiiken fc Gilmord Il989l . 
IChen et al.,,2001. .Sieeel et al.,,2002 ). An excellent recent 
example of the applic ation of this metho d to pre-SDSS 
data is the study by ISiegel et all (|2002f ). While their 
and SDSS data are of similar photometric quality, the 
sky area analyzed here is over 400 times larger than that 
analyzed by Siegel et al. This large increase in sample 
size enables a shift in emphasis from modelling to direct 
model-free mapping of the complex and clumpy Galac- 
tic density distribution. Such mapping and analysis of 
the maps allows for identification and removal of clumpy 
substructure, which is a necessary precondition for a re- 
liable determination of the functional form and best-fit 
parameters of the Galactic model. 

This qualitative paradigm shift was made possible by 
the availability of SDSS data. SDSS is superior to previ- 
ous optical sky surveys because of its high catalog com- 
pleteness and precise multi-band CCD photometry to 

■^■^ Note that the polar rings explanation is also unlikely for the- 
oretical reasons as these are thought to originate in large galactic 
collisions which would leave its imprint on other components of the 
Milky Way as well. We discuss it as an option here from purely 
observational standpoint, and mainly for completeness. 
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Virgo overdensity Control field Counts difference 




Fig. 39. — Hess diagrams of Pis color vs. r magnitude for the Virgo overdensity field (left), the control field as defined on figure l37l 
(middle), and their difference (right). The colors encode star counts within the fields. A significant excess of stars with Pis < —0.2 is 
present in the Virgo overdensity field. There is no statistically significant difference in star counts for stars having Pis > —0.2, implying 
that the stars that constitute of the Virgo overdensity have metallicities lower than disk stars, and closer to metallicities characteristic of 
the halo. 



faint flux limits over a large sky area. In particular, the 
results presented here were enabled by several distinctive 
SDSS characteristics: 

• A large majority of stars detected by the SDSS 
are main-sequence stars, which have a fairly well- 
defined color-luminosity relation. Thus, accurate 
SDSS colors can be used to estimate luminosity, 
and hence, distance, for each individual star. Ac- 
curate photometry (~ 0.02 mag) allows us to reach 
the intrinsic accuracy of photometric parallax rela- 
tion, and estimate distances to single stars within 
15-20%, and estimate the relative distances of stars 
in clumps of similar age and metallicity to better 
than 5%. 

• Thanks to faint flux limits (r ~ 22), distances 
as large as 15-20 kpc are probed using numerous 
main sequence stars 48 million). At the same 
time, the large photometric dynamic range and the 
strong dependence of stellar luminosities on color 
allow constraints ranging from the Sun's offset from 
the Galactic plane 25 pc) to a detection of over- 
densities at distances beyond 10 kpc. 

• Large sky area observed by the SDSS (as opposed 
to pencil beam surveys), spanning a range of galac- 
tic longitudes and latitudes, enables not only a 
good coverage of the (i?, Z) plane, but also of a 
large fraction of the Galaxy's volume. The full 
three-dimensional analysis, such as slices of the 
maps in X—Y planes, reveals a great level of detail. 

• The SDSS u band photometric observations can be 
used to identify stars with sub-solar metallicities, 
and to study the differences between their distri- 
bution and that of more metal-rich stars. 

6.2. The Best-Fit Galactic Model 

When the exponential disk models are used to describe 
the gross behavior of the stellar number density distribu- 



tion, we derive the best-flt parameter values summarized 
in Table [lOl 

Before proceeding to compare these results to the lit- 
erature, we note that a proper comparison with previous 
work is sometimes difficult due to the lack of clarity of 
some authors regarding which effects were (or were not) 
taken into account when deriving the model parameters. 
Of particular concern is the problem of unrecognized 
multiplicity: uncorrected for, or if using a significantly 
different binary fraction, it will affect the disk scales by 
up to 30% (Section l4.3.4p . In the discussion to follow we 
assumed, if not explicitly mentioned otherwise, that all 
appropriate corrections were taken into account by the 
authors of the studies against which compare our results. 

The derived 300 pc vertical scale of the thin disk 
(corrected for an assumed 35% binary fraction) is 
about 10% lower than the cannonical 325pc value, 
and near the middle of the range of values found 
in th e recent literature (240 — 350pc. iRobiii et ed] 
('19961: 'Larson fc Humphrevs ('1996'):'Bus er et al.l (|19990: 
Chen et al. (2001); Sicg el et al. ( 2002,) ). Similarly, 
the scale height of the th ic k disk is in the range 
found bv ISiegel et all (l2002f) : iBuser et all (|1999D and 
iLarsen fc Humphreys! (| 19961 ). and about 20% higher than 
the 580 — 790pc rang e spanned by measurements of 
iRobin et all ([1996); Oi ha et all (|1999D and IChen et all 
( 2001h . We note that uncorrected for unrecognized mul- 
tiplicity, our thin and the thick disk scale estimates (245 
and 740, respectively) would be closer to the lower end 
of the range found in the literature. 

We find the local thick disk normalization of ^ 12%, 
larger than most pr evious estim a tes bu t simi lar to recent 
determinations by IChen et al.l (|200lD and ISiegel et al.l 
(j2002) (> 10%). Models with normalizations lower than 
10% show increasingly large and, in particular, the 
combinations of parameters characteristic of early "low 
normalization/high thick disk scale height" models (e.g., 
Gilmore fc Rfiidl[l983t IRobin fc Crez d[l986l: lYoshii et all 
1987t lYamagata fc Yoshiil ll992: Rcid fc Maiewskilll993f) 
are strongly disfavored by the SDSS data. The root cause 
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Fig. 40.— The sky distribution of 189 2MASS M giant candidates with b > 45°, selected by 9.2< K <10.2 and 1.0< J - K <1.3. The 
symbols in the top panel are color-coded using their K band magnitude and J — K color, according to the scheme shown in the bottom 
panel (which shows all 75,735 candidates from the whole sky). The symbols in bottom right panel show the same sample as in the top 
panel in Lambert projection, with the SDSS density map from Fig.X shown as the gray scale background. At sin(b) > 0.8, there are 2.5 
times as many stars with I < than with I > 0. This asymmetry provides an independent confirmation of the Virgo overdensity revealed 
by the SDSS data. 
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Fig. 41. — The mass contribution to thin (solid) and thick (dot- 
ted) disks from different radii and heights in the Galaxy. Center 
panel shows the isodensity contours of thin (solid) and thick (dot- 
ted) disk having the bias-corrected parameters of Table 1101 Two 
rectangles centered around R = 8 kpc enclose the range in radii 
and vertical distances Z from the disk from which the model pa- 
rameters were measured. The bottom panel shows the cumulative 
fraction of disk mass enclosed outside a given radius R. Similarly, 
the side panel shows the fraction of disk mass enclosed at heights 
\Z\ > |.Zgivcnl- Note that while our coverage in Z direction is ade- 
quate, in the R direction we sample a region that contains less than 
20% of the total disk mass, and extrapolate the obtained density 
laws towards the Galactic center where most of the mass lies. 

of the apparent discrepancy may be sought for in the fact 
that all of these studies were pencil-beam surveys oper- 
ating on a single or at most a few lines of sight, usually 
towards the NGP. However, a single or even a few pencil 
beams are insufficient to break the degeneracies inher- 
ent in the multiparameter Galactic model (Section l4.3.9p . 
While adequately describing the observed lines of sight, 
these pencil-beam best-fit parameters are local minima, 
unrepresentative of the entire Galaxy. Only by using a 
wider and deeper sample, such as the one presented here, 
were we able to break the degeneracy and derive a glob- 
ally representative model. 
The value of thin disk scale length i s in ag r eement with 



cale length I S m ag r eement with 
Oiha et all (|1999h . IChen et all 



20021 ). and lower than the tra- 
4kpc. The scale length of the 



the recent estimates by 
l|200ll ) and iSiegel etldl 
ditionally assumed 3 — 
thick disk is longer than that of the thin disk. The qual- 
itative nature of this result is robust: variations of the 
assumed photometric parallax relation, binary fraction 
or the exact choice of and size of the color bins, leave 
it unchanged. Quantitatively, the ratio of best-fit length 
scales is close to 1.4, similar (within uncertainties) to typ- 
ical scale leng th ratios of ^ 1.25 seen in edg e-on late- type 
disk galaxies (jYoachim &: Dalcanton|[2"006r ). 

Assuming that exponential density laws correctly de- 
scribe the density distribution all the way to the Galactic 
center, our model implies that ~ 23% of the total lumi- 
nosity (and stellar mass) in K and M-dwarfs is contained 
in the thick disk. Despite being an extrapolation from a 
small region containing only a few tens of percent of the 
total mass of the disk (Figure [JT]), this is in good agree- 
ment wit h observatio ns of thick disks in external edge-on 
galaxies (jYoachim fc Dalcantonll2006L Figure 24). 



6.3. Detection of Resolved Substructure in the Disk 

Although important for understanding the current 
structure of the Milky Way, the disk mass ratios, details 
of density laws, or the exact values of Galactic model 
parameters are insufficient by themselves to attack the 
question of mechanism of Galactic disk formation (both 
thin and thick) . It is the departures from these laws that 
actually hold more information about the formation than 
the laws themselves. 

Thick disk formation scenarios can be broadly be di- 
vided into three classes: (1) slow kinematic heating of 
stars from the thin disk (Spitzer fc Schwarzschild 19531; 



iBarbanis fc Woltied Il967t iSellwood fc Carlberel a984). 



(2) pressure-supported slow collapse i mmedia.tely a fter 
an ELS-like monolithic collapse (e.g., iLarsonl [l976| ) or 

(3) merger- induced stirring of thin disk material and/or 
direct accretion of stellar content of the progenito r 
()Quinn et a"nil993HAbadi et al.ll2003HBrook et al.ll2004D . 
Scenarios (1) and (2) are usually disfavored due to the 
inability of either the giant molecular clouds or the spi- 
ral structure t o excite the stars to orbit s observed in the 
thin disk (e.g. iQuillen fc Garnettll200d[ l. and the appar- 
ent lack of vertical m etallicity gradient in the thick disk 
(|Gilmore et al.|[T995t see however Paper II for evidence 
to the contrary). 

The third scenario recently garnered increased 
attention, with detailed theoretical simulations of 
the formation of realistic galaxie s in ACDM hier- 
archical me r ger p icture context (jAbadi et al.l l2003t 
iBrook et al.l 120041 ) . an d the observation of prop- 
erties of thick di sks (^D alcanton fc BernsteinI 120021 : 
lYoachim fc Dalcant on 2006) and even a counter-rota ting 
thick disk in FGC 227 (Yoachim fc Dalcanton] [20051). A 
simulation reported by Abadi et al. (2003), while not 
directly comparable to the Milky Way (their galaxy is 
spheroid-dominated) is especially illuminating regarding 
the qualitative mechanisms that may build up the thick 
disk. Three of their conclusions are of particular conse- 
quence to our work: i) the thick disk is formed by direct 
accretion of stellar content from satellites on low incli- 
nation orbits, ii) the stars from a single disrupted satel- 
lite are not uniformly radially mixed, but rather form 
a torus-like structure at radii where the final disruption 
occurs, and iii) if formed through the same process, the 
disk of the Milky Way disk may still hold signatures of 
such early accretion events. 

Our finding that the thin and thick disk structure, 
similarly to that of the halo, is complicated by local- 
ized overdensities and permeated by ring-like departures 
from exponential pr ofiles may lend cred ence to the mech- 
anism described bv lAbadi et al.l (|2003f ). In addition to 
already known Monoceros stream, we found evidence for 
two more overdensities in the thick disk region (Fig- 
ure l27p . both consistent with rings or streams in the den- 
sity maps. While unlikely to be the relics from the age 
of thick disk formation (they would need to survive for 
~ 8 — 10 Gyr) it is plausible they, like the Monoceros 
stream, are remnants of smaller and more recent accre- 
tion events analogous to those that formed the thick disk. 

In case of the Monoceros stream, the three-dimensional 
maps offer an effective method to study its properties. 
The maps demonstrate this feature is well localized in the 
radial direction, which rules out the hypothesis that this 
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overdensity is due to disk flaring. The maps also show 
that the Monoceros stream is not a homogeneously dense 
ring that surrounds the Galaxy, pro viding support for the 
claim by iRocha-Pinto et al.1 ()2003D that this structure is 
a merging dwarf galaxy (see also Pcnarr ubia et al]l2005l 
for a comprehensive theoretical model). In Paper II, we 
demonstrate that stars in the Monoceros stream have 
metallicity distribution more metal-poor than thick disk 
stars, but more metal-rich than halo stars. 

Discoveries of this type of substructure point to a pic- 
ture of the thick disk filled with streams and remnants 
much in the same way like the halo. A crude extrapola- 
tion of three disk overdensities seen in our survey volume 
{\Z\ < 3 kpc, i? < 15 kpc) to the full Galactic disk leads 
to a conclusion that there may be up to ^--^15 - 30 clumpy 
substructures of this type in the Galaxy. These "disk 
streams" are still likely to carry, both in their physical 
(metallicity) and kinematic properties, some information 
on their progenitors and history. 

6.4. Stellar Halo 

We find it possible to describe the halo of the Milky 
Way by an oblate r""" power-law ellipsoid, with the 
axis ratio c/a = qn ^ 0.5 — 0.8 and the power-law in- 
dex of uh = 2.5 — 3 (with the formal best fit parameters 
qh = 0.64 and hh — 2.8 for galactocentric radii <20 
kpc). These values are consistent with previous studies: 
specifically, they are in excellent agreement with Besan- 
con prograin valu es {qh = 0.6 — 0.85, nn = 2.44 — 2.75; 
iRobin et al.ll2000h . with a more recent measurement of 
QH = 0.6 and ur = 2.75 bv lSiegel et all (120021 and with 
the p revious SDSS estimate of qh ~ 0.55 (jChen et al.l 
[2001]) . The convergence of best-fit values is encourag- 
ing, especially considering the differences in methodolo- 
gies (direct fitting vs. population synthesis modelling) 
and the data (photometric systems, limiting magnitudes, 
types of tracers, and lines of sight) used in each of these 
studies. 

The goodness of halo model fit is poorer than that 
of the disk fits (reduced ~ 2 — 3). Similar prob- 
lems with halo fits were previously noticed in studies of 
smaller sa mples of kinematically and metallicity-selected 
field stars llHartwick "1987':'S ommer-Larsen fc Zhe nlll990l; 
Allen et a lj 119911: Prestou 'eTaLl 119911: Kinma n et all 
1994c iCarnev et al. 1996; Chiba & Beers 200^, globu- 
lar clusters (Zinn 11993 ; Dincscu et al. 1999 ) and mam 
sequence stars ([Gilmore et all Il985l : ISiegel et al.ll2002[ ). 
and are unlikely to be explaine d away by instrumental 
or methodological reasons alone (jSiegel et al. 2002). Our 
own analysis of why this is so (Section 14.3.81 and Fig- 
ure I30|) points towards a need for a more complex den- 
sity distribution profile. For example, instead of a single 
power law, a two-component "dual-halo" , in which the 
stars are divided into a spherical and a flattened sub- 
comp onent, may be inv oked to explain the observations 
(e.g.. ISommer-Larsen fc Zhcn 1990). 

Such models, when applied to starc ounts, do show im - 
provements over a single power law (jSiegel et al.ll2002f ). 
Furthermore, this division may be theo retically moti- 
vated by an attempt to unify the ELS and lSearle fc Ziiml 
1197^1) pictures of Galaxy formation: the flattened sub- 
component being a result of the initial monolithic col- 
lapse, and the spherical compo nent originating from sub- 
sequent accretion of satellites (|Sandagelll990f . iMaiewskil 



ll993tlNorrisll994f) . While this explanation is circumstan- 
tially supported by the detection of on going accret ion in 
the halo today (e . g. lYannv et al."2000': Ivezic et al. ll2000l: 
Vivas et al.]|2001l : iMaiewski e t al. 2003; Bclokurov et al I 



20061 and references therein) , we would prefer a more di- 



rect line of evidence for it, derived from observations of 
halo stars themselves. 

For example, one may hope the component coming 
from accretion is visibly irregular, streamlike, and/or 
clumpy, thus lending credence to the hypothesis of its 
origin. However, our examination of the distribution of 
residuals in Section l4.3.8l revealed no signal of unresolved 
clumpy substructure in the halo on 1 — 2 ~ kpc scales. 
Instead, we found the large reduced is best explained 
by a poor choice of density law profile (a single power 
law) . A double power-law, or a more complicat ed pro- 
file such as the one used bv lPreston et al.l (|1991h . would 
likely better fit the data. 

The dumpiness may still be prevalent, but on a differ- 
ent spatial scale, or smalled in amplitude and harder to 
detect with a simple analysis employed here. We leave a 
detailed study of scale-dependent dumpiness in the halo 
and its possible two-component nature for a subsequent 
study. 

6.5. The Virgo Overdensity 

We report the discovery of Virgo overdensity. De- 
spite its large angular size and proximity, its low surface 
brightness kept it from being recognized by smaller sur- 
veys. Given the low surface brightness, its well defined 
outline, and low metallicity, the most plausible explana- 
tion of Virgo overdensity is that it is a result of a merger 
event involving the Milky Way and a smaller, lower- 
metallicity dwarf galaxy. For now, based on existing 
maps, we are unable to differentiate whether the observed 
overdensity is a tidal stream, a merger remnant, or both. 
However, it is evident that the Virgo overdensity is sur- 
prisingly large, extending in vertical (Z) direction to the 
boundaries of our survey (6 < Z < 15 kpc), and ^ 10 kpc 
in R direction. It is also exceedingly faint, with a lower 
limit on surface brightness of = 32.5 mag arcsec"^. 

A potential connection of Virgo overdensity and the 
Sagittarius stream is di s cussed in a foUowup paper by 
iMartfnez-Delgado et"aI1 (|2007f ). Their N-body simula- 
tions of the Sagittarius stream show that the Virgo over- 
density resides in the region of space where the leading 
arm of the Sagittarius stream is predicted to cross the 
Milky Way plane in the Solar neighborhood. This tenta- 
tive Virgo-Sagittarius association needs to be confirmed 
by measurement of highly negative radial velocities for 
the stars of the Virgo overdensity. 

A similar diffuse structure, the Triangulum- 
Andromeda feature (hereafter, TriAn d) , wa s re- 
cently identified by iRocha-Pinto et al"] (|2004D and 
iMajewski et al.l ()2004f ) in the southern Galactic hemi- 
sphere, as an overdensity of M giants observed with 
2MASS. They find an excess in M giant number den- 
sity over a large area of the sky (100° < / < 150°, 
-40° < 6 < -20°). TriAnd, just as the Virgo structure 
presen ted here, is very diffuse and shows no clear 
core. iRocha-Pinto et all (|2004[ ) estimate the distance 
to TriAnd of at least ~ 10 kpc. Recently, additional 
tenuous str uctures were discovere d in the same regio n 
of the sky (jMaiewski et all 120041 : iMartin et a"Lll2007f) . 
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pointing to the possibility that dilFuse clouds such as 
Virgo and TriAnd are quite common in the Galactic 
halo. 

Assuming that the Virgo overdensity is a part of a 
larger previously unidentified stream, it would be of in- 
terest to look for a possible continuation in the south- 
ern Galactic hemisphere. Our preliminary analysis of 
2MASS M-giants data did not reveal a similarly large 
density enhancement in the south. It would also be in- 
teresting to follow the stream towards the Galactic north, 
beyond the Z ~ 20 kpc limit of our survey, where a sig- 
nature o£_overdeiisit;Y_Ji^s been revealed by RR Lyrae 
stars (|Duffau et al.ll2006h . Above all, the understand- 
ing of the Virgo overdensity would greatly benefit from 
measurements of proper motion and radial velocity of its 
constituent stars. 

6.6. Mapping the Future 

This study is only a first step towards a better under- 
standing of the Milky Way enabled by modern large-scale 
surveys. Star counting, whether interpreted with tradi- 
tional modeling methods, or using number density maps, 
is limited by the number of observed stars, the flux limit 
and sky coverage of a survey, and the ability to differen- 
tiate stellar populations. All these data aspects will soon 
be significantly improved. 

First, the SDSS has entered its second phase, with 
a significant fraction of observing time allocated 
for the Milky Way studies (SEGUE, the Sloan Ex- 
tension for Galaxy Understanding and Explorat ion, 
Newberg fc Sloan Digital Sky Survev Collaboration! 
2003( 1. In particular, the imaging of low galactic lati- 
tudes and a large number of stellar spectra optimized for 
Galactic structure studies will add valuable new data to 
complement this work. In addition, the SDSS kinematic 
data, both from radial velocities and from proper 
motions (determined from astrometric comparison of 
the SD SS and the Palomar Observatory Sky Survey 
catalog. iMunn et~al . 2004) is already yielding significant 
advances in our understanding of the thin and thick 
disk, and halo kinematic structure (Paper III, in prep.). 

Another improvement to the analysis presented here 
will come from th e GAIA satellite mission (e.g. 
[Wilkinson et al.ll2005f) . GAIA will provide geometric dis- 
tance estimates and spectrophotometric measurements 
for a large number of stars brighter than ~ 20. Despite 
the relatively bright flux limit, these data will be invalu- 
able for calibrating photometric parallax relation, and 
for studying the effects of metallicity, binarity and con- 
tamination by giants. At the moment, the uncertainties 
of the photometric parallax relation are the single largest 
contributor to uncertainties in the derived parameters of 
Galactic models, and improvements in its calibration are 
of great interest to all practitioners in this field. 

A further major leap forward will be enabled by up- 
coming deep synopt ic sky surveys, such as Pan-STARRS 
(iKaiser et al. l Hool) and LSST (|Tvsonl [200l . Pan- 
STARRS has already achieved first light with its first 



telescope, and the four-telescope version may become 
operational around 2010. If approved for construction 
in 2009, the LSST may obtain its first Hght in 2014. 
These surveys will provide multi-band optical photom- 
etry of better quality than SDSS over practically the en- 
tire sky (LSST wiU be sited in Chile, and Pan-STARRS 
is at Hawaii; note that Pan-STARRS will not use the 
u band filter). One of their advantages will be signifi- 
cantly deeper data - for example, the LSST will enable 
studies such as this one to a 5 magnitudes fainter limit, 
corresponding to a distance limit of 150 kpc for the turn- 
off stars. LSST proper motion measurements will con- 
strain tangential velocity to within 10 km/s at distances 
as large as that of the Virgo overdensity reported here 
(~10 kpc). These next-generation maps will be based on 
samples including several billion stars and will facilitate 
not only the accurate tomography of the Milky Way, but 
of the whole Local Group. 



Acknowledgments 

We thank Princeton University, the University of 
Washington and the Institute for Advanced Study for 
generous financial support of this research. M. Juric 
gratefully acknowledges support from the Taplin Fellow- 
ship and from NSF grant PHY-0503584. Z. Ivezic and B. 
Sesar acknowledge support by NSF grant AST-0551161 
to LSST for design and development activity. Donald 
P. Schneider acknowledges support by NSF grant AST- 
0607634. We especially thank the anonymous referee for 
numerous helpful comments and suggestions which have 
significantly improved this manuscript. 

Funding for the creation and distribution of the SDSS 
Archive has been provided by the Alfred P. Sloan Foun- 
dation, the Participating Institutions, the National Aero- 
nautics and Space Administration, the National Sci- 
ence Foundation, the U.S. Department of Energy, the 
Japanese Monbukagakusho, and the Max Planck Soci- 
ety. The SDSS Web site is http://www.sdss.org/. 

The SDSS is managed by the Astrophysical Research 
Consortium (ARC) for the Participating Institutions. 
The Participating Institutions are The University of 
Chicago, Fermilab, the Institute for Advanced Study, 
the Japan Participation Group, The Johns Hopkins Uni- 
versity, the Korean Scientist Group, Los Alamos Na- 
tional Laboratory, the Max-Planck-Institute for Astron- 
omy (MPIA) , the Max-Planck-Institutc for Astrophysics 
(MPA), New Mexico State University, University of 
Pittsburgh, University of Portsmouth, Princeton Uni- 
versity, the United States Naval Observatory, and the 
University of Washington. 

This publication makes use of data products from the 
Two Micron All Sky Survey, which is a joint project of 
the University of Massachusetts and the Infrared Pro- 
cessing and Analysis Center/California Institute of Tech- 
nology, funded by the National Aeronautics and Space 
Administration and the National Science Foundation. 



REFERENCES 

Abadi, M. G., Navarro, J. F., Steinmetz, M., & Eke, V. R. 2003, Abazajian, K., Adelman-McCarthy, J. K., Agiieros, M. A., AUam, 

ApJ, 597, 21 S. S., Anderson, K. S. J., Anderson, S. F., Annis, J., Bahcall, 

N. A., Baldry, I. K., Bastian, S., Berlind, A., Bernardi, 
M., Blanton, M. R., Bochanski, Jr., J. J., Boroski, W. N., 
Brewington, H. J., Briggs, J. W., Brinkmann, J., Brunner, 
R. J., Budavari, T., Carey, L. N., Castander, F. J., Connolly, 
A. J., Covey, K. R., Csabai, I., Dalcanton, J. J., Doi, M., 
Dong, F., Eisenstein, D. J., Evans, M. L., Fan, X., Finkbeiner, 
D. P.. Friedman. S. D.. Frieman. J. A.. Fukueita. M.. GillesDie. 



52 



Juric et al. 



Abazajian, K., Adclman-McCarthy, J. K., Agiieros, M. A., AUam, 
S. S., Anderson, S. F., Annis, J., Bahcall, N. A., Baldry, I. K., 
Bastian, S., Berlind, A., Bernard!, M., Blanton, M. R., Blythe, 
N., Bochanski, Jr., J. J., Boroski, W. N., Brcwington, H., 
Briggs, J. W., Brinkmann, J., Brunner, R. J., Budavari, T., 
Carey, L. N., Carr, M. A., Castander, F. J., Chiu, K., CoUinge, 
M. J., Connolly, A. J., Covey, K. R., Csabai, I., Dalcanton, 
J. J., Dodelson, S., Doi, M., Dong, F., Eisenstein, D. J., Evans, 
M. L., Fan, X., Feldman, P. D., Finkbeiner, D. P., Friedman, 
S. D., Frieman, J. A., Fukugita, M., Gal, R. R., Gillespie, 
B., Glazebrook, K., Gonzalez, C. F., Gray, J., Grebel, E. K., 
Grodnicki, L., Gunn, J. E., Gurbani, V. K., Hall, P. B., Hao, 
L., Harbeck, D., Harris, F. H., Harris, H. C, Harvanek, M., 
Hawley, S. L., Hcckman, T. M., Hclmboldt, J. F., Hendry, J. S., 
Hennessy, G. S., Hindslcy, R. B., Hogg, D. W., Holmgren, D. J., 
Holtzman, J. A., Homer, L., Hui, L., Ichikawa, S.-i., Ichikawa, T., 
Inkmann, J. P., Ivezic, Z., Jester, S., Johnston, D. E., Jordan, 
B., Jordan, W. P., Jorgenscn, A. M., Juric, M., KaufTmanii, G., 
Kent, S. M., Kleinman, S. J., Knapp, G. R., Kniazcv, A. Y., 
Kron, R. G., Krzcsinski, J., Kunszt, P. Z., Kuropatkin, N., Lamb, 

D. Q., Lampcitl, H., Laubsclicr, B. E., Lcc, B. C, Lcgcr, R. F., 
Li, N., Lidz, A., Lin, H., Loli, Y.-S., Long, D. C., Lovcday, J., 
Lupton, R. H., Malik, T., Margon, B., McGchcc, P. M., McKay, 
T. A., Meiksin, A., Miknaitis, G. A., Moorthy, B. K., Munn, 
J. A., Murphy, T., Nakajima, R., Narayanan, V. K., Nash, T., 
Ncilscn, Jr., E. H., Ncwbcrg, H. J., Newman, P. R., Nichol, R. C, 
Nicinski, T., Nicto-Santistcban, M., Nitta, A., Odcnkirchen, M., 
Okamura, S., Ostrikcr, J. P., Owen, R., Padmanabhan, N., 
Peoples, J., Pier, J. R., Pindor, B., Pope, A. C., Quinn, T. R., 
Rafikov, R. R., Raymond, S. N., Richards, G. T., Richmond, 
M. W., Rix, H.-W., Rockosi. C. AI., Schayc, J., Schlcgcl, D. J., 
Schneider, D. P., Sclirocdci, .)., Scranton, R., Sckiguchi, M., 
SeljaJf, U., Sergey, G., Sesar, B., Sheldon, E., Sliimasaku, K., 
Siegmund, W. A., Silvestri, N. M., Sinisgalli, A. J., Sirko, 

E. , Smith, J. A., Smolcic, V., Snedden, S. A., Stebbins, A., 
Steinhardt, C, Stinson, G., Stoughton, C, Strateva, L V., 
Strauss, M. A., SubbaRao, M., Szalay, A. S., Szapudi, L, Szkody, 
P., Tasca, L., Tegmark, M., Thakar, A. R., Tremonti, C, Tucker, 
D. L., Uomoto, A., Vanden Berk, D. E., Vandenberg, J., Vogeley, 
M. S., Voges, W., Vogt, N. P., Walkowicz, L. M., Weinberg, D. H., 
West, A. A., White, S. D. M., Wilhite, B. C, Willman, B., Xu, 
Y., Yanny, B., Yarger, J., Yasuda, N., Yip, C.-W., Yocum, D. R., 
York, D. G., Zakamska, N. L., Zehavi, L, Zheng, W., Zibetti, S., 
& Zucker, D. B. 2003, AJ, 126, 2081 

Adelman-McCarthy, J. K., Agiieros, M. A., AUam, S. S., Anderson, 
K. S. J., Anderson, S. F., Annis, J., Bahcall, N. A., Baldry, 
I. K., Barcntinc, J. C, Berlind, A., Bcrnardi, M., Blanton, M. R., 
Boroski, W. N., Brcwington, H. J., Brinchmann, J., Brinkmann, 
J., Brunner, R. J., Budavari, T., Carey, L. N., Carr, M. A., 
Castander, F. J., Connolly, A. J., Csabai, I., Czarapata, P. C, 
Dalcanton, J. J., Doi, M., Dong, F., Eisenstein, D. J., Evans, 
M. L., Fan, X., Finkbeiner, D. P., Friedman, S. D., Fricman, 
J. A., Fukugita, M., Gillespie, B., Glazebrook, K., Gray, J., 
Grebel, E. K., Gimu, .J. E., Gurbani, V. K., dc Haas, E., Hall, 
P. B., Harris, F. H., Harvanek, M., Hawley, S. L., Hayes, J., 
Hendry, J. S., Hennessy, G. S., Hindslcy, R. B., Hirata, C. M., 
Hogan, C. J., Hogg, D. W., Holmgren, D. J., Holtzman, J. A., 
Ichikawa, S.-i., Ivezic, Z., Jester, S., Johnston, D. E., Jorgenscn, 
A. M., Juric, M., Kent, S. M., Kleinman, S. J., Knapp, G. R., 
Kniazev, A. Y., Kron, R. G., Krzcsinski, J., Kuropatkin, N., 
Lamb, D. Q., Lampcitl, H., Lee, B. C, Lcgcr, R. F., Lin, H., 
Long, D. C, Loveday, J., Lupton, R,. H., Margon, B., Martfnez- 
Delgado, D., Mandelbaum, R., Matsubara, T., McGehee, P. M., 
McKay, T. A., Meiksin, A., Munn, J. A., Nakajima, R., Nash, T., 
Neilsen, Jr., E. H., Newberg, H. J., Newman, P. R., Nichol, R. C, 
Nicinski, T., Nieto-Santisteban, M., Nitta, A., O'MuUanc, W., 
Okamura, S., Owen, R., Padmanabhan, N., Pauls, G., Peoples, 
J. J., Pier, J. R., Pope, A. C, Pourbaix, D., Quinn, T. R., 
Richards, G. T., Richmond, M. W., Rockosi, C. M., Schlegel, 
D. J., Schneider, D. P., Schroeder, J., Scranton, R., Seljak, U., 
Sheldon, E., Shimasaku, K., Smith, J. A., Smolcic, V., Snedden, 
S. A., Stoughton, C, Strauss, M. A., SubbaRao, M., Szalay, 

A. S., Szapudi, I., Szkody, P., Tegmark, M., Thakar, A. R., 
Tucker, D. L., Uomoto, A., Vanden Berk, D. E., Vandenberg, 
J., Vogeley, M. S., Voges, W., Vogt, N. P., Walkowicz, L. M., 
Weinberg, D. H., West, A. A., White, S. D. M., Xu, Y., Yanny, 

B. , Yocum, D. R., York, D. G., Zehavi, I., Zibetti, S., & Zucker, 
D. B. 2006, ApJS, 162, 38 



Allen, C, Poveda, A., & Schuster, W. J. 1991, A&A, 244, 280 
AUende Pricto, C, Beers, T. C, Wilhclm, R., Newberg, H. J., 

Rockosi, C. M., Yanny, B., & Lee, Y. S. 2006, ApJ, 636, 804 
Baade, W. 1944, ApJ, 100, 137 
Bahcall, J. N., & Soneira, R. M. 1980, ApJS, 44, 73 
Barbanis, B., & Woltjer, L. 1967, ApJ, 150, 461 
Belokurov, V., Zucker, D. B., Evans, N. W., Gilmore, G., Vidrih, 

S., Bramich, D. M., Newberg, H. J., Wyse, R. F. G., Irwin, M. J., 

Fellhauer, M., Hewett, P. C, Walton, N. A., Wilkinson, M. I., 

Cole, N., Yanny, B., Rockosi, C. M., Beers, T. C, Bell, E. F., 

Brinkmann, J., Ivezic, Z., & Lupton, R. 2006, ApJ, 642, L137 
Bilir, S., KaraaH, S., Ak, S., Yaz, E., & Hamzaoglu, E. 2006, New 

Astronomy, 12, 234 
Binggcli, B. 1999, in Lecture Notes in Physics, Berlin Springer 

Vcrlag, Vol. 530, The Radio Galaxy Messier 87, ed. H.-J. Roser 

& K. Meisenheimcr, 9 — h 
Binney, J., & Merrifield, M. 1998, Galactic astronomy (Galactic 

astronomy / James Binney and Michael Merrifield. Princeton, 

NJ : Princeton University Press, 1998. (Princeton series in 

astrophysics) QB857 .B522 1998 ($35.00)) 
Binney, J., & Trcmaine, S. 1987, Galactic dynamics (Princeton, 

NJ, Princeton University Press, 1987, 747 p.) 
Brook, C. B., Kawata, D., Gibson, B. K., & Freeman, K. C. 2004, 

ApJ, 612, 894 
Buscr, R., Rong, J., & Karaali, S. 1999, A&A, 348, 98 
Carney, B. W. 1979, ApJ, 233, 211 

Carney, B. W., Laird, J. B., Latham, D. W., & Aguilar, L. A. 1996, 

AJ, 112, 668 

Chen, B., Figucras, F., Torra, J., Jordi, C, Luri, X., & Galadi- 

Enrfqucz, D. 1999, A&A, 352, 459 
Chen, B., Stoughton, C, Smith, J. A., Uomoto, A., Pier, J. R., 
Yanny, B., Ivezic, Z., York, D. G., Anderson, J. E., Annis, J., 
Brinkmann, J., Csabai, I., Fukugita, M., Hindsley, R., Lupton, 
R., Munn, J. A., & the SDSS Collaboration. 2001, ApJ, 553, 184 
Chiba, M., & Beers, T. C. 2000, AJ, 119, 2843 
Dalcanton, J. J., & Bernstein, R. A. 2002, AJ, 124, 1328 
Dinescu, D. I., Girard, T. M., & van Altena, W. F. 1999, A J, 117, 
1792 

Duflfau, S., Zinn, R., Vivas, A. K., Carraro, G., Mendez, R. A., 

Winnick, R., & Gallart, C. 2006, ApJ, 636, L97 
Duquennoy, A., & Mayor, M. 1991, A&A, 248, 485 
Eggen, O. J., Lynden-Bell, D., & Sandage, A. R. 1962, ApJ, 136, 
748 

Eisenstein, D. J., Liebert, J., Harris, H. C, Kleinman, S. J., Nitta, 

A., Silvestri, N., Anderson, S. A., Barentine, J. C, Brcwington, 
H. J., Brinkmann, J., Harvanek, M., Krzcsiiiski, J., Ncilscn, Jr., 
E. H., Long, D., Schneider, D. P., & Snedden, S. A. 2006, ApJS, 
167, 40 
Fan, X. 1999, AJ, 117, 2528 

Finkbeiner, D. P., Padmanabhan, N., Schlcgcl, D. J., Carr, M. A., 
Gunn, J. E., Rockosi, C. M., Sekiguchi, M., Lupton, R. H., 
Knapp, G. R., Ivezic, Z., Blanton, M. R., Hogg, D. W., Adelman- 
McCarthy, J. K., Annis, J., Hayes, J., Kinney, E., Long, D. C, 
Scljak, U., Strauss, M. A., Yanny, B., Agiieros, M. A., AUam, 
S. S., Anderson, S. F., Bahcall, N. A., Baldry, I. K., Bcrnardi, 
M., Boroski, W. N., Briggs, J. W., Brinkmann, J., Brunner, 
R. J., Budavari, T., Castander, F. J., Covey, K. R., Csabai, I., 
Doi, M., Dong, F., Eisenstein, D. J., Fan, X., Friedman, S. D., 
Fukugita, M., Gillespie, B., Grebel, E. K., Gurbani, V. K., de 
Haas, E., Harris, F. H., Hendry, J. S., Hennessy, G. S., Jester, 
S., Johnston, D. E., Jorgenscn, A. M., Juric, M., Kent, S. M., 
Kniazev, A. Y., Krzosiriski, J., Leger, R. F., Lin, H., Lovcday, J., 
Mannery, E., Martfnez-Delgado, D., McGehee, P. M., Meiksin, 
A., Munn, J. A., Neilsen, Jr., E. H., Newman, P. R., Nitta, 
A., Pauls, G., Quinn, T. R., Rafikov, R. R., Richards, G. T., 
Richmond, M. W., Schneider, D. P., Schroeder, J., Shimasaku, 
K., Siegmund, W. A., Smith, J. A., Snedden, S. A., Stebbins, 
A., Szalay, A. S., Szokoly, G. P., Tegmark, M., Tucker, D. L., 
Uomoto, A., Vanden Berk, D. E., Weinberg, D. H., West, A. A., 
Yasuda, N., Yocum, D. R., York, D. G., & Zehavi, I. 2004, A J, 
128, 2577 

Finlator, K., Ivezic, Z., Fan, X., Strauss, M. A., Knapp, G. R., 
Lupton, R. H., Gunn, J. E., Rockosi, C. M., Anderson, J. E., 
Csabai, I., Hennessy, G. S., Hindsley, R. B., McKay, T. A., 
Nichol, R. C, Schneider, D. P., Smith, J. A., York, D. G., & 
the SDSS Collaboration. 2000, AJ, 120, 2615 

Fischer, D. A., & Marcy, G. W. 1992, ApJ, 396, 178 



The Milky Way Stellar Number Density Distribution 



53 



Fukugita, M., Ichikawa, T., Gunn, J. E., Doi, M., Shimasaku, K., 

& Schneider, D. P. 1996, AJ, 111, 1748 
Gilmore, G., & Reid, N. 1983, MNRAS, 202, 1025 
Gilmorc, G., Rcid, N., & Hewett, P. 1985, MNRAS, 213, 257 
Gilmore, G., Wyse, R. F. G., & Jones, J. B. 1995, AJ, 109, 1095 
Gilmore, G., Wyse, R. F. G., & Kuijken, K. 1989, ARA&A, 27, 

555 

Gott, J. R. I., Juric, M., Schlegel, D., Hoyle, F., Vogeley, M., 
Tegmark, M., Bahcall, N., & Brinkmann, J. 2005, ApJ, 624, 463 
Gunn, J. E., Carr, M., Rockosi, C., Sekiguchi, M., Berry, K., Elms, 

B. , de Haas, E., Ivezic, Z., Knapp, G., Lupton, R., Pauls, G., 
Simcoe, R., Hirsch, R., Sanford, D., Wang, S., York, D., Harris, 

F. , Annis, J., Bartozek, L., Boroski, W., Bakken, J., Haldeman, 
M., Kent, S., Holm, S., Holmgren, D., Petravick, D., Prosapio, 

A. , Rcchcnmachcr, R., Doi, M., Fukugita, M., Shimasaku, K., 
Okada, N., Hull, C., Siegmund, W., Mannery, E., Blouke, M., 
Heidtman, D., Schneider, D., Lucinio, R., & Brinkman, J. 1998, 
AJ, 116, 3040 

Gunn, J. E., Siegmund, W. A., Mannery, E. J., Owen, R. E., 
Hull, C. L., Leger, R. F., Carey, L. N., Knapp, G. R., York, 
D. G., Boroski, W. N., Kent, S. M., Lupton, R. H., Rockosi, 

C. M., Evans, M. L., WaddcU, P., Anderson, J. E., Annis, J., 
Barcntinc, J. C., Bartoszck, L. M., Bastian, S., Brackcr, S. B., 
Brcwington, H. J., Bricgcl, C. I., Brinkmann, J., Brown, Y. J., 
Carr, M. A., Czarapata, P. C, Drcnnan, C. C, Dombcck, T., 
Fcdcrwitz, G. R., Gillespie, B. A., Gonzales, C, Hansen, S. U., 
Harvanek, M., Hayes, J., Jordan, W., Kinney, E., Klacnc, M., 
Kleinman, S. J., Kron, R. G., Kresinski, J., Lee, G., Limmongkol, 
S., Lindcnmeyer, C. W., Long, D. C, Loomis, C. L., McGehee, 
P. M., Mantsch, P. M., Neilsen, Jr., E. H., Neswold, R. M., 
Newman, P. R., Nitta, A., Peoples, J. J., Pier, J. R., Prieto, 
P. S., Prosapio, A., Rivetta, C, Schneider, D. P., Snedden, S., & 
Wang, S.-i. 2006, AJ, 131, 2332 

Harris, W. E. 1996, AJ, 112, 1487 

Hartwick, F. D. A. 1987, in NATO ASIC Proc. 207: The Galaxy, 

ed. G. Gilmore & B. Carswell, 281-290 
Hogg, D. W., Finkbeiner, D. P., Schlegel, D. J., & Gunn, J. E. 

2001, AJ, 122, 2129 
Ibata, R. A., Irwin, M. J., Lewis, G. F., Ferguson, A. M. N., & 

Tanvir, N. 2003, MNRAS, 340, L21 
Ivezic, Z., AUyn Smith, J., Miknaitis, G., Lin, H., & Tucker, D. 

2007a^ ArXiv Astrophysics e-prints 
Ivezic, Z., Goldston, J., Finlator, K., Knapp, G. R., Yanny, 

B. , McKay, T. A., Amrose, S., Krisciunas, K., Willman, B., 
Anderson, S., Schaber, C, Erb, D., Logan, C, Stubbs, C, Chen, 
B., Neilsen, E., Uomoto, A., Pier, J. R., Fan, X., Gunn, J. E., 
Lupton, R. H., Rockosi, C. M., Schlegel, D., Strauss, M. A., 
Annis, J., Brinkmann, J., Csabai, I., Doi, M., Fukugita, M., 
Hennessy, G. S., Hindsley, R. B., Margon, B., Munn, J. A., 
Newberg, H. J., Schneider, D. P., Smith, J. A., Szokoly, G. P., 
Thakar, A. R., Vogclcy M. S., WaddcU, P., Yasiida, N., & York, 

D. G. 2000, AJ, 120, 963 

Ivezic, Z., Lupton, R., Johnston, D., Richards, G., Hall, P., 
Schlegel, D., Fan, X., Munn, J., Yanny, B., Strauss, M., Knapp, 

G. , Gunn, J., & Schneider, D. 2004a, in Astronomical Society of 
the Pacific Conference Scries, Vol. 311, AGN Physics with the 
Sloan Digital Sky Survey, cd. G. T. Richards & P. B. Hall, 437-+ 

Ivezic, Z., Lupton, R. H., Anderson, S., Eyer, L., Gunn, J. E., Juric, 
M., Knapp, G. R., Miknaitis, G., Gunn, J. E., Rockosi, C. M., 
Schlegel, D., Strauss, M. A., Stubbs, C, & Vanden Berk, D. E. 
2003, Memorie della Socicta Astronomica Italiana, 74, 978 

Ivezic, Z., Lupton, R,. H., Schlegel, D., Boroski, B., Adclman- 
McCarthy, J., Yaimy, B., Kent, S., Stoughton, C, Finkbeiner, 
D., Padmanabhan, N., Rockosi, C. M., Gunn, J. E., Knapp, 
G. R., Strauss, M. A., Richards, G. T., Eisenstein, D., Nicinski, 
T., Kleinman, S. J., Krzesinski, J., Newman, P. R., Snedden, S., 
Thakar, A. R., Szalay, A., Munn, J. A., Smith, J. A., Tucker, D., 
& Lee, B. C. 2004b, Astronomische Nachrichten, 325, 583 

Ivezic, Z., Lupton, R. H., Schlegel, D., Smolcic, V., Johnston, D., 
Gunn, J. E., Knapp, G. R., Strauss, M. A., & Rockosi, C. M. 
2003a, ArXiv Astrophysics e-prints 

— . 2003b, ArXiv Astrophysics e-prints 



Ivezic, Z., Menou, K., Knapp, G. R., Strauss, M. A., Lupton, R. H., 
Vanden Berk, D. E., Richards, G. T., Tremonti, C, Weinstein, 
M. A., Anderson, S., Bahcall, N. A., Becker, R. H., Bernardi, M., 
Blanton, M., Eisenstein, D., Fan, X., Finkbeiner, D., Finlator, 
K., Prieman, J., Gunn, J. E., Hall, P. B., Kim, R. S. J., 
Kinkhabwala, A., Narayanan, V. K., Rockosi, C. M., Schlegel, 
D., Schneider, D. P., Strateva, I., SubbaRao, M., Thakar, A. R., 
Voges, W., White, R. L., Yanny, B., Brinkmann, J., Doi, M., 
Fukugita, M., Hennessy, G. S., Munn, J. A., Nichol, R. C., & 
York, D. G. 2002, AJ, 124, 2364 

Ivezic, Z., Schlegel, D., Uomoto, A., Bond, N., Beers, T., Allende 
Prieto, C, Wilhelm, R., Lee, Y. S., Sivarani, T., Juric, M., 
Lupton, R., Rockosi, C, Knapp, G., Gunn, J., Yanny, B., Jester, 
S., Kent, S., Pier, J., Munn, J., Richards, G., Newberg, H., 
Blanton, M., Eisenstein, D., Hawlcy, S., Anderson, S., Harris, 
H., Kiuchi, F., Chen, A., Bushong, J., Sohi, H., Haggard, 
D., Kimball, A., Barentine, J., Brcwington, H., Harvanek, M., 
Kleinman, S., Krzesinski, J., Long, D., Nitta, A., Snedden, S., & 
for the SDSS Collaboration. 2007b, ArXiv Astrophysics c-prints 

Ivezic, Z., Smith, J. A., Miknaitis, G., Lin, H., Tucker, D., Lupton, 
R., Knapp, G., Gunn, J., Strauss, M., Holtzman, J., Kent, 
S., Yanny, B., Schlegel, D., Finkbeiner, D., Padmanabhan, N., 
Rockosi, C, Juric, M., Bond, N., Lee, B., Jester, S., Harris, 
H., Harding, P., Brinkmann, J., York, D., & for the SDSS 
Collaboration. 2007c, ArXiv Astrophysics c-prints 

Ivezic, Z., Tabachnik, S., Rafikov, R., Lupton, R. H., Quinn, T., 
Hammcrgren, M., Eyer, L., Chu, J., Armstrong, J. C, Fan, 
X., Finlator, K., Geballc, T. R., Gunn, J. E., Hennessy, G. S., 
Knapp, G. R., Leggett, S. K., Munn, J. A., Pier, J. R., Rockosi, 
C. M., Schneider, D. P., Strauss, M. A., Yanny, B., Brinkmann, 
J., Csabai, I., Hindsley, R. B., Kent, S., Lamb, D. Q., Margon, 
B., McKay, T. A., Smith, J. A., Waddel, P., York, D. G., & the 
SDSS Collaboration. 2001, AJ, 122, 2749 

Ivezic, Z., Vivas, A. K., Lupton, R. H., & Zinn, R. 2005, AJ, 129, 
1096 

Jackson, T., Ivezic, Z., & Knapp, G. R. 2002, MNRAS, 337, 749 
Juric, M., Ivezic, Z., Lupton, R. H., Quinn, T., Tabachnik, S., Fan, 
X., Gunn, J. E., Hennessy, G. S., Knapp, G. R., Munn, J. A., 
Pier, J. R., Rockosi, C. M., Schneider, D. P., Brinkmann, J., 
Csabai, I., & Fukugita, M. 2002, A J, 124, 1776 
Kaiser, N., Aussel, H., Burke, B. E., Boesgaard, H., Chambers, K., 
Chun, M. R., Heasley, J. N., Hodapp, K.-W., Hunt, B., Jedicke, 
R., Jewitt, D., Kudritzki, R., Luppino, G. A., Maberry, M., 
Magnier, E., Monet, D. G., Onaka, P. M., Pickles, A. J., Rhoads, 
P. H. H., Simon, T., Szalay, A., Szapudi, I., Tholen, D. J., 
Tonry, J. L., Waterson, M., & Wick, J. 2002, in Presented at 
the Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference, Vol. 4836, Survey and Other Telescope Technologies 
and Discoveries. Edited by Tyson, J. Anthony; Wolff, Sidney. 
Proceedings of the SPIE, Volume 4836, pp. 154-164 (2002)., ed. 
J. A. Tyson & S. Wolff, 154-164 
Karaali, S., BiUr, S., & Tunccl, S. 2005, Publications of the 

Astronomical Society of Australia, 22, 24 
Karaali, S., Karatas, Y., BiUr, S., Ak, S. G., & Hamzaoglu, E. 2003, 

Publications of the Astronomical Society of Australia, 20, 270 
Kinman, T. D., Suntzeff, N. B., & Kraft, R. P. 1994, AJ, 108, 1722 
Kleinman, S. J., Harris, H. C, Eisenstein, D. J., Liebert, J., Nitta, 
A., Krzesinski, J., Munn, J. A., Dahn, C. C, Hawley, S. L., 
Pier, J. R., Schmidt, G., Silvcstri, N. M., Smith, J. A., Szkody, 
P., Strauss, M. A., Knapp, G. R., CoUingc, M. J., Mukadam, 
A. S., Koester, D., Uomoto, A., Schlegel, D. J., Anderson, S. F., 
Brinkmann, J., Lamb, D. Q., Schneider, D. P., & York, D. G. 
2004, ApJ, 607, 426 
Kroupa, P., Tout, C. A., & Gilmore, G. 1993, MNRAS, 262, 545 
Kuijken, K., & Gilmore, G. 1989, MNRAS, 239, 605 
Laird, J. B., Carney, B. W., & Latham, D. W. 1988, AJ, 95, 1843 
Larsen, J. A. 1996, PhD thesis, PhD Thesis, University of 

Minnesota, (1996) 
Larsen, J. A., & Humphreys, R. M. 1996, ApJ, 468, L99-I- 
Larson, R. B. 1976, MNRAS, 176, 31 

Lupton, R. H., Ivezic, Z., Gunn, J. E., Knapp, G., Strauss, M. A., 
& Yasuda, N. 2002, in Presented at the Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference, Vol. 4836, Survey 
and Other Telescope Technologies and Discoveries. Edited by 
Tyson, J. Anthony; WolfF, Sidney. Proceedings of the SPIE, 
Volume 4836. pp. 350-356 (2002)., ed. J. A. Tyson & S. Wolff, 
350-356 

Luyten, W. J. 1968, MNRAS, 139, 221 



54 



Juric et al. 



Mai'z-Apellaniz, J. 2001, AJ, 121, 2737 
Majewski, S. R. 1993, ARA&A, 31, 575 

Majewski, S. R., Ostheimer, J. C, Rocha-Pinto, H. J., Patterson, 
R. J., Guhathakurta, P., & Reitzcl, D. 2004, ApJ, 615, 738 

Majewski, S. R., Skrutskie, M. F., Weinberg, M. D., &; Ostheimer, 
J. C. 2003, ApJ, 599, 1082 

Martin, N. F., Ibata, R. A., & Irwin, M. 2007, ArXiv Astrophysics 
e-prints 

Maxtmez-Delgado, D., Alonso-Garcia, J., Aparicio, A., & Gomez- 

Flechoso, M. A. 2001, ApJ, 549, L63 
Maxtmez-Delgado, D., Gomez-Flechoso, M. A., Aparicio, A., & 

Carrera, R. 2004, ApJ, 601, 242 
Maxtmez-Delgado, D., Peiiaxrubia, J., Juric, M., Alfaxo, E. J., & 

Ivezic, Z. 2007, ApJ, 660, 1264 
Munn, J. A., Monet, D. G., Levine, S. E., Canzian, B., Pier, 

J. R., Haxris, H. C., Lupton, R. H., Ivezic, Z., Hindsley, R. B., 

Hennessy, G. S., Schneider, D. P., & Brinkmann, J. 2004, A J, 

127, 3034 

Newberg, H. J., & Sloan Digital Sky Survey Collaboration. 2003, in 
Bulletin of the American Astronomical Society, Vol. 35, Bulletin 
of the American Astronomical Society, 1385 — h 

Newberg, H. J., & Yanny, B. 2005, in Astronomical Society of the 
Pacific Conference Series, Vol. 338, Astrometry in the Age of 
the Next Generation of Large Telescopes, cd. P. K. Seidclmann 
& A. K. B. Monet, 210-+ 

Newberg, H. J., Yanny, B., Rockosi, C, Grcbel, E. K., Rix, H.-W., 
Brinkmann, J., Csabai, I., Hennessy, G., Hindsley, R. B., Ibata, 
R., Ivezic, Z., Lamb, D., Nash, E. T., Odenkirchen, M., Rave, 
H. A., Schneider, D. P., Smith, J. A., Stolte, A., & York, D. G. 
2002, ApJ, 569, 245 

Norris, J. E. 1994, ApJ, 431, 645 

O'Connell, D. J. K. 1958, Ricerche Astronomiche, 5 

Ojha, D. K., Bienayme, O., Mohan, V., & Robin, A. C. 1999, A&A, 
351, 945 

Oke, J. B., & Gunn, J. E. 1983, ApJ, 266, 713 

Parker, J. E., Humphreys, R. M., & Larscn, J. A. 2003, AJ, 126, 
1346 

Peiiaxrubia, J., Martfnez-Delgado, D., Rix, H. W., Gomez- 
Flechoso, M. A., Munn, J., Newberg, H., Bell, E. F., Yanny, 
B., Zucker, D., & Grebel, E. K. 2005, ApJ, 626, 128 

Pier, J. R., Munn, J. A., Hindsley, R. B., Hennessy, G. S., Kent, 
S. M., Lupton, R. H., & Ivezic, Z. 2003, AJ, 125,' 1559 

Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 
1992, Numerical recipes in C. The axt of scientific computing 
(Cambridge: University Press, — cl992, 2nd ed.) 

Preston, G. W., Shectman, S. A., & Beers, T. C. 1991, ApJ, 375, 
121 

Quillen, A. C, & Garnett, D. R. 2000, ArXiv Astrophysics e-prints 

Quinn, P. J., Hcrnquist, L., & FuUagar, D. P. 1993, ApJ, 403, 74 

Reid, I. N., & Gizis, ,1. E. 1997, AJ, 113, 2246 

Reid, I. N., Gizis, J. E., & Hawley, S. L. 2002, AJ, 124, 2721 

Reid, I. N., Lcwitus, E., Allen, P. R., Cruz, K. L., & Burgasser, 
A. J. 2006, AJ, 132, 891 

Reid, I. N., van Wyk, F., Marang, F., Roberts, G., Kilkenny, D., 
& Mahoney, S. 2001, MNRAS, 325, 931 

Reid, I. N., Yan, L., Majewski, S., Thompson, I., & Small, I. 1996, 
AJ, 112, 1472 

Reid, M. J. 1993, ARA&A, 31, 345 

Reid, N., & Majewski, S. R. 1993, ApJ, 409, 635 

Richards, G. T., Fan, X., Newberg, H. J., Strauss, M. A., Vanden 
Berk, D. E., Schneider, D. P., Yanny, B., Boucher, A., Buries, 
S., Frieman, J. A., Gunn, J. E., Hall, P. B., Ivezic, Z., Kent, 
S., Loveday, J., Lupton, R. H., Rockosi, C. M., Schlegel, D. J., 
Stoughton, C, SubbaRao, M., & York, D. G. 2002, AJ, 123, 2945 

Robin, A., & Creze, M. 1986, A&A, 157, 71 

Robin, A. C, Haywood, M., Creze, M., Ojha, D. K., & Bienayme, 

O. 1996, A&A, 305, 125 
Robin, A. C, Rcylc, C, & Creze, M. 2000, A&A, 359, 103 
Rocha-Pinto, H. J., Majewski, S. R., Skrutskie, M. F., & Crane, 

J. D. 2003, ApJ, 594, L115 
Rocha-Pinto, H. J., Majewski, S. R., Skrutskie, M. F., Crane, J. D., 

& Patterson, R. J. 2004, ApJ, 615, 732 
Sandage, A. 1990, JRASC, 84, 70 

Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 

Seaxle, L., & Zinn, R. 1978, ApJ, 225, 357 

Sellwood, J. A., & Caxlberg, R. G. 1984, ApJ, 282, 61 



Sesar, B., Svilkovic, D., Ivezic, Z., Lupton, R. H., Munn, J. A., 
Finkbeiner, D., Steinhardt, W., Siverd, R., Johnston, D. E., 
Knapp, G. R., Gunn, J. E., Rockosi, C. M., Schlegel, D., Vanden 
Berk, D. E., Hall, P., Schneider, D. P., & Brunner, R. J. 2006, 
AJ, 131, 2801 

Siegel, M. H., Majewski, S. R., Reid, I. N., & Thompson, I. B. 

2002, ApJ, 578, 151 
Smith, J. A., Tucker, D. L., Kent, S., Richmond, M. W., Pukugita, 
M., Ichikawa, T., Ichikawa, S.-i., Jorgensen, A. M., Uomoto, A., 
Gunn, J. E., Hamabe, M., Watanabe, M., Tolea, A., Henden, A., 
Annis, J., Pier, J. R., McKay, T. A., Brinkmann, J., Chen, B., 
Holtzman, J., Shimasaku, K., & York, D. G. 2002, AJ, 123, 2121 
Smolcic, v., Ivezic, Z., Knapp, G. R., Lupton, R. H., Pavlovski, K., 
Ilijic, S., Schlegel, D., Smith, J. A., McGehee, P. M., Silvestri, 
N. M., Hawley, S. L., Rockosi, C, Gunn, J. E., Strauss, M. A., 
Fan, X., Eiscnstein, D., & Harris, H. 2004, ApJ, 615, L141 
Sommer-Larscn, ,1., & Zhcn, C. 1990, MNRAS, 242, 10 
Spitzer, L. J., & Schwarzschild, M. 1953, ApJ, 118, 106 
Stoughton, C, Lupton, R. H., Bernardi, M., Blanton, M. R., 
Buries, S., Castander, F. J., Connolly, A. J., Eisenstein, D. J., 
Frieman, J. A., Hemiessy, G. S., Hindsley, R. B., Ivezic, Z., Kent, 
S., Kunszt, P. Z., Lee, B. C, Meiksin, A., Munn, J. A., Newberg, 
H. J., Nichol, R. C, Nicinski, T., Pier, J. R., Richards, G. T., 
Richmond, M. W., Schlegel, D. J., Smith, J. A., Strauss, M. A., 
SubbaRao, M., Szalay A. S., Thakar, A. R., Tucker, D. L., 
Vanden Berk, D. E., Yanny, B., Adelman, J. K., Anderson, 
Jr., J. E., Anderson, S. F., Annis, J., Bahcall, N. A., Bakken, 
J. A., Bartelmann, M., Bastian, S., Bauer, A., Berman, E., 
Bohringer, H., Boroski, W. N., Bracker, S., Briegel, C, Briggs, 
J. W., Brinkmann, J., Brunner, R., Carey, L., Carr, M. A., Chen, 
B., Christian, D., Colestock, P. L., Crocker, J. H., Csabai, I., 
Czarapata, P. C, Dalcanton, J., Davidsen, A. F., Davis, J. E., 
Dehnen, W., Dodelson, S., Doi, M., Dombeck, T., Donahue, 
M., EUman, N., Elms, B. R., Evans, M. L., Eyer, L., Fan, X., 
Federwitz, G. R., Friedman, S., Fukugita, M., Gal, R., Gillespie, 

B. , Glaaebrook, K., Gray, J., Grebel, E. K., Greenawalt, B., 
Greene, G., Gunn, J. E., de Haas, E., Haiman, Z., Haldeman, 
M., Hall, P. B., Hamabe, M., Hansen, B., Harris, F. H., Harris, 
H., Haxvanek, M., Hawley, S. L., Hayes, J. J. E., Heckmaxi, T. M., 
Helmi, A., Henden, A., Hogan, C. J., Hogg, D. W., Holmgren, 

D. J., Holtzman, J., Huang, C.-H., Hull, C., Ichikawa, S.-L, 
Ichikawa, T., Johnston, D. E., KaufTmann, G., Kim, R. S. J., 
Kimball, T., Kinney, E., Klaene, M., Kleinman, S. J., Klypin, A., 
Knapp, G. R., Korienek, J., Krolik, J., Kron, R. G., Krzesinski, 
J., Lamb, D. Q., Lcger, R. F., Limmongkol, S., Lindenmeyer, C, 
Long, D. C, Looniis, C, Loveday, J., MacKinnon, B., Mannery, 

E. J., Mantsch, P. M., Margon, B., McGehee, P., McKay, 
T. A., McLean, B., Menou, K., Merelli, A., Mo, H. J., Monet, 

D. G., Nakamura, O., Narayanan, V. K., Nash, T., Ncilsen, Jr., 

E. H., Newman, P. R., Nitta, A., Odenkirchen, M., Okada, N., 
Okamura, S., Ostriker, J. P., Owen, R., Pauls, A. G., Peoples, J., 
Peterson, R. S., Petravick, D., Pope, A., Pordes, R., Postman, 
M., Prosapio, A., Quinn, T. R., Rechenmachcr, R., Rivetta, 

C. H., Rix, H.-W., Rockosi, C. M., Rosner, R., Ruthmansdorfer, 
K., Sandford, D., Schneider, D. P., Scranton, R., Sekiguchi, M., 
Sergey, G., Sheth, R., Shimasaku, K., Smee, S., Sneddon, S. A., 
Stebbins, A., Stubbs, C, Szapudi, I., Szkody, P., Szokoly, G. P., 
Tabachnik, S., Tsvetanov, Z., Uomoto, A., Vogeley, M. S., Voges, 
W., Waddell, P., Walterbos, R., Wang, S.-i., Watanabe, M., 
Weinberg, D. H., White, R. L., White, S. D. M., Wilhite, B., 
Wolfe, D., Yasuda, N., York, D. G., Zehavi, I., & Zheng, W. 
2002, AJ, 123, 485 

Tucker, D. L., Kent, S., Richmond, M. W., Annis, J., Smith, J. A., 
Allam, S. S., Rodgers, C. T., Stute, J. L., Adelman-McCarthy, 
J. K., Brinkmann, J., Doi, M., Finkbeiner, D., Fukugita, M., 
Goldston, J., Greenway, B., Gunn, J. E., Hendry, J. S., Hogg, 

D. W., Ichikawa, S.-L, Ivezic, Z., Knapp, G. R., Lampeitl, 
H., Lee, B. C, Lin, H., McKay, T. A., Merrelh, A., Munn, 
J. A., Neilsen, Jr., E. H., Newberg, H. J., Richards, G. T., 
Schlegel, D. J., Stoughton, C., Uomoto, A., & Yanny, B. 2006, 
Astronomische Nachrichten, 327, 821 

Tyson, J. A. 2002, in Presented at the Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference, Vol. 4836, Survey 
and Other Telescope Technologies and Discoveries. Edited by 
Tyson, J. Anthony; Wolff, Sidney. Proceedings of the SPIE, 
Volume 4836, pp. 10-20 (2002)., ed. J. A. Tyson & S. Wolff, 
10-20 



The Milky Way Stellar Number Density Distribution 



55 



Vivas, A. K., Zinn, R., Andrews, P., Bailyn, C, Baltay, C, Coppi, 
P., EUman, N., Girard, T., Rabinowitz, D., Schaefer, B., Shin, J., 
Snyder, J., Sofia, S., van Altena, W., Abad, C, Bongiovanni, A., 
Bricerio, C, Bruzual, G., Delia Prugna, P., Herrera, D., Magris, 

G. , Mateu, J., Pacheco, R., Sanchez, G., Sanchez, G., Schenner, 

H. , Stock, J., Vicente, B., Vieira, K., Ferrfn, I., Hernandez, 
J., Gebhard, M., Honeycutt, R., Mufson, S., Musser, J., & 
Rengstorf, A. 2001, ApJ, 554, L33 

Wilkinson, M. I., Vallenari, A., Turon, C., Munari, U., Katz, D., 
Bono, G., Cropper, M., Helmi, A., Robichon, N., Thevenin, 

F. , Vidrih, S., Zwitter, T., Arenou, P., Baylac, M.-O., BerteUi, 

G. , Bijaoui, A., Boschi, P., Castelli, P., Crifo, P., David, M., 
Gomboc, A., Gomez, A., Haywood, M., Jauregi, U., de Laverny, 
P., Lebreton, Y., Marrese, P., Marsh, T., Mignot, S., Morin, D., 
Pasetto, S., Perryman, M., Prsa, A., Recio-Blanco, A., Royer, 
P., Sellier, A., Siviero, A., Sordo, R,., Soubiran, C., Tomasella, 
L., & Viala, Y. 2005, MNRAS, 359, 1306 

Williams, C. C., Golimowski, D. A., Uomoto, A., Reid, I. N., 
Henry, T. J., Dieterich, S., Jue, S. L., Long, G. M., Neilsen, 
E. H., Spahn, E. Y., & Walkowicz, L. M. 2002, in Bulletin 
of the American Astronomical Society, Vol. 34, Bulletin of the 
American Astronomical Society, 1292 — h 

Wozniak, P. R., Vestrand, W. T., Akerlof, C. W., Balsano, R., 
Bloch, J., Casperson, D., Pletcher, S., Gisler, G., Kehoe, R., 
Kinemuchi, K., Lee, B. C., Marshall, S., McGowan, K. E., 
McKay, T. A., Rykoff, E. S., Smith, D. A., Szymanski, J., & 
Wren, J. 2004, AJ, 127, 2436 

Wyse, R. P. G., & Gilmore, G. 1989, Comments on Astrophysics, 
13, 135 

Yamagata, T., & Yoshii, Y. 1992, AJ, 103, 117 

Yanny, B., Newberg, H. J., Kent, S., Laurent-Muehleisen, S. A., 
Pier, J. R., Richards, G. T., Stoughton, C, Anderson, Jr., J. E., 
Annis, J., Brinkmann, J., Chen, B., Csabai, L, Doi, M., Pukugita, 
M., Hennessy, G. S., Ivezic, Z., Knapp, G. R., Lupton, R., Munn, 
J. A., Nash, T., Rockosi, C. M., Schneider, D. P., Smith, J. A., 
& York, D. G. 2000, ApJ, 540, 825 

Yoachim, P., & Dalcanton, J. J. 2005, ApJ, 624, 701 

— . 2006, AJ, 131, 226 

York, D. G., Adelman, J., Anderson, Jr., J. E., Anderson, S. F., 
Annis, J., Bahcall, N. A., Bakken, J. A., Barkhouser, R., Bastian, 
S., Berman, E., Boroski, W. N., Bracker, S., Briegel, C, Briggs, 
J. W., Brinkmann, J., Brunner, R., Buries, S., Carey, L., Carr, 
M. A., Castander, F. J., Chen, B., Colestock, P. L., Connolly, 
A. J., Crocker, J. H., Csabai, I., Czarapata, P. C, Davis, J. E., 
Doi, M., Dombeck, T., Eisenstein, D., Ellman, N., Elms, B. R., 
Evans, M. L., Pan, X., Federwitz, G. R., Piscelli, L., Friedman, 
S., Frieman, J. A., Pukugita, M., Gillespie, B., Gunn, J. E., 
Gurbani, V. K., de Haas, E., Haldeman, M., Harris, F. H., 
Hayes, J., Heckman, T. M., Hennessy, G. S., Hindsley, R. B., 
Holm, S., Holmgren, D. J., Huang, C.-h., Hull, C, Husby, D., 
Ichikawa, S.-L, Ichikawa, T., Ivezic, Z., Kent, S., Kim, R. S. J., 
Kinney, E., Klaene, M., Kleinman, A. N., Kleinman, S., Knapp, 
G. R., Korienek, J., Kron, R. G., Kunszt, P. Z., Lamb, D. Q., 
Lee, B., Leger, R. P., Limmongkol, S., Lindenmeyer, C, Long, 
D. C, Loomis, C, Loveday, J., Lucinio, R., Lupton, R. H., 
MacKinnon, B., Mannery, E. J., Mantsch, P. M., Margon, B., 
McGehee, P., McKay, T. A., Meiksin, A., MereUi, A., Monet, 
D. G., Munn, J. A., Narayanan, V. K., Nash, T., Neilsen, E., 
Neswold, R., Newberg, H. J., Nichol, R. C, Nicinski, T., Nonino, 
M., Okada, N., Okamura, S., Ostriker, J. P., Owen, R., Pauls, 
A. G., Peoples, J., Peterson, R. L., Petravick, D., Pier, J. R., 
Pope, A., Pordes, R., Prosapio, A., Rechenmacher, R., Quinn, 
T. R., Richards, G. T., Richmond, M. W., Rivetta, C. H., 
Rockosi, C. M., Ruthmansdorfer, K., Sandford, D., Schlegel, 
D. J., Schneider, D. P., Sekiguchi, M., Sergey, G., Shimasaku, K., 
Siegmund, W. A., Smee, S., Smith, J. A., Snedden, S., Stone, R., 
Stoughton, C, Strauss, M. A., Stubbs, C, SubbaRao, M., Szalay, 

A. S., Szapudi, L, Szokoly, G. P., Thakar, A. R., Tremonti, C, 
Tucker, D. L., Uomoto, A., Vanden Berk, D., Vogeley, M. S., 
Waddell, P., Wang, S.-i., Watanabe, M., Weinberg, D. H., Yanny, 

B. , & Yasuda, N. 2000, AJ, 120, 1579 



Yoshii, Y., Ishida, K., & Stobie, R. S. 1987, AJ, 93, 323 

Zinn, R. 1993, in Astronomical Society of the Pacific Conference 

Series, Vol. 48, The Globular Cluster-Galaxy Connection, ed. 

G. H. Smith & J. P. Brodie, 38-+ 
Zinn, R., Katherina Vivas, A., Gallart, C, & Winnick, R. 2003, 

ArXiv Astrophysics c-prints 



APPENDIX 

EFFECTS OP LOCUS PROJECTION 

The improvement in the estimate of r — i color resulting from the locus projection depends on the local slope of 
the locus. If the locus has a steep or almost vertical slope, as for stars with 5 — r ~ 1.4 (cf. figure [8]), the knowledge 
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oi g — r color does not further constrain r — i color. On the other hand, for shallow slopes the knowledge of 5 — r 
color determines the intrinsic r — i to a much better accuracy than the the r — i measurement alone. For most of the 
observed g — r color range we are closer to the second regime, with the locus having a slope of d{r — i)/d{g — r) ~ 0.3 
for < g — r < 1 . 

To illustrate the effects of locus projection, in Fig. |42] we simulate an ensemble of 10^ stars with the same color 
{g — r = 0.4, r — i = 0.143) subjected to photometric errors representative of SDSS observations (a^-i = fg-r = 
0.03 mag). Errors introduce a scatter in observed colors as shown in the inset. Using only the r — i information we 
obtain the expected ct^-z = 0.03 mag scatter in observed r — i color (dashed histogram). By projecting the colors to 
the locus, the scatter is reduced to fir-i — 0.01 mag. As we construct density maps for stars binned in color bins with 
A(r — i) = 0.05 mag width, this is a significant reduction in scatter. 

There are other benefits of locus projection. Figure l44l illustrates how photometric errors can bias the determination 
of the number of stars in regions with large gradients of the number density (in the g — r vs. r — i color-color diagram). 
The solid line shows the true density distribution of stars as a function oi r — i color. In this toy model, we used 
p{r — i) (X {r ~ z)^, with a sharp cutoff at r — i = 0.24 (solid histogram). Scattering the stars drawn from that 
distribution with the <Tr-i = f^g-r = 0.03 mag photometric errors, and binning in r — i color bins, produces the dashed 
histogram (effectively, a convolution of the original distribution with the SDSS photometric errors). It systematically 
underestimates the original density distribution by as much as 50% in the region of highest negative gradient (near 
the r — i = 0.24 cutoff), while overestimating the density in regions with positive density gradient. The exact under or 
overestimation depends on the scale at which the true density changes appreciably. If it is significantly smaller than 
Ur-i or (Tg-r, 8iS IS the case near the cutoff in figure [44l the biases become significant. Applying the locus projection 
gives a much better estimate of the original distribution (dotted histogram). 

This is not a purely hypothetical example - this type of strong density gradient is observed near g — r ^ 0.2 (cf. 
figure El) . Figure |44] compares the number of stars per r — i color bin for observed r — i colors (dashed histogram) and 
the colors estimated by locus projection (solid histogram) for the 48 million stars in our star catalog. In the regions 
of maximal gradients the locus projection has a ~ 20% effect on the total number of stars per r — i bin. Not including 
this correction would similarly bias the density normalization of Galactic models deduced from this sample. The bias 
would further propagate to luminosity function determination"^^. 

One may further be concerned about the dependence of the location of the r — z vs. g — r stellar locus on other 
parameters. We do observe a slight magnitude dependence of the locus, especially for stars with bluer g — r colors. 
Figure [45] shows the dependence of the measured r — i color on the r band magnitude for stars with 0.29 < g — r < 0.31 
(the "worst case g — r color", where we find the dependence of the locus location on the magnitude to be the largest). 
We find a weak, approximately linear dependence, with d(r — i)iocus/dr ~ 0.007 mag/mag (solid line). The horizontal 
dashed line shows the location of the locus for this bin as given by eq. El The two match for r = 18. We have chosen 
to disregard this magnitude dependence when performing the locus projection procedure. We reason that the real 
dependence is not on magnitude but on metallicity'^^ and probably other unknown factors for which we do not have a 
firm handle. By attempting to correct for those, we risk introducing additional unknown and more complicated biases. 
Secondly, the magnitude dependence is relatively small and only affects the bluest bins. 

Finally, the three panels in Fig. US] show the histograms of observed r — i colors for stars which have their locus 
projected colors in the 0.1 < (r — i)e < 0.15 color range. The top panel shows the histogram for all stars in the sample. 
The 0.03 mag scatter is comparable to SDSS photometric errors. To check for the effects discussed in the previous 
paragraph, the center and bottom panels show histograms for the brightest and faintest magnitude bin respectively. 
Slight magnitude dependence can be seen as a small shift of histogram median to the left (center panel) and right 
(top panel). Also, the worsening of photometric precision at the faint end is quite visible in the bottom panel, as the 
scatter increases to 0.08 mag. 



For this particular work, the discussed bias is not a problem as we make use of only the reddest bins {r — i > 1.0) to measure the 
Galactic model parameters 

For example, a part of the shift in locus is likely created by the distant metal-poor stars of the halo, which are preferentially found at 
fainter magnitudes in our sample. 
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Fig. 42. — An illustration of the reduction in the r — i color scatter after applying the locus projection. An ensemble of 10^ stars with 
the same color {g — r = 0.4, r — i = 0.143) is subjected to SDSS photometric errors (cTg = rTr = (7i = 0.02 mag), resulting in the color 
distribution shown in the inset. Its r — j histogram (dashed line) has a root-mean-square scatter of a, — i = 0.03 mag. After the colors are 
locus corrected (solid histogram), the scatter is reduced to cr, — i = 0.01 mag. The amount of reduction of the scatter depends on the slope 
of the locus - the more horizontal the locus gets, the better is the determination of the true r — i color. 
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Fig. 43. — An illustration of the bias in the determination of the number density introduced by photometric errors. The solid histogram 
shows a toy model distribution as a function of the r — i color (p(r — j) <x (r — i)^, with a sharp cutoff at r — i = 0.24). The dashed 
histogram is the convolution of the original distribution with SDSS photometric errors [ug = = Oi = 0.02 mag). The measured density 
of stars with r — i = 0.24 is underestimated by as much as 50%. Applying the locus projection gives an improved estimate of the original 
distribution (dotted histogram). 
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r— i 



Fig. 44. — A comparison of the number of stars per r — i bin for observed r — i (dashed histogram) and the locus corrected (r — i)e colors 
(solid histogram) for the 48 million stars in our star catalog. In the regions of maximal gradients the locus projection has a ~ 20% effect 
on the total number of stars per r ~ i bin. Not including this correction would similarly bias the density normalization of Galactic models. 
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Fig. 45. — The dependence of the measured r — i color on the r band magnitude for stars with 0.29 < g — r < 0.31. An approximately 
linear dependence, with d(r — i)/dr ~0.007 mag/mag is observed (solid line). The horizontal dashed line shows the location of the locus 
for this bin as given by eq. \5\ The two match for r = 18. 
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Fig. 46. — The effect of ignoring the magnitude dependence of the locus. The histograms show observed r — i colors of stars having 
locus-corrected colors in the 0.1 < (r — i)e < 0.15 color range. The top panel shows the histogram for all stars in the sample. The spread 
of (T = O.OS™''^ is comparable to SDSS photometric errors. The middle and bottom panels show histograms for the brightest and faintest 
magnitude bin, respectively. A weak magnitude dependence can be seen as a small shift of histogram median to the left (middle panel) 
and right (bottom panel). 



